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ABSTRACT 


This  report  is  a  description  of  a  computer  code  that  results  in  the  first-  and 
higher-order  statistics  of  RCS  of  a  complex  target  as  a  function  of  frequency, 
polarization,  aspect  angle,  and  constituent  parameters  of  the  target.  The 
premise  is  that  a  complex  target  can  be  represented  by  a  set  of  simple  scatter¬ 
ed.  There  is  a  general  discussion  of  RCS  data;  the  capabilities  of  the  code 
from  the  viewpoint  of  the  RCS  analyst  or  data  librarian;  the  theory  of  RCS, 
statistics,  lobe  structure,  and  glint  embodied  in  the  code;  a  user’s  guide  from 
the  viewpoint  of  the  computer  analyst;  and  a  program  listing. 
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1.0  BACKGROUND 


The  computer  code  described  in  this  report  -  The 
APL  RCS/Statistics  Code  -  represents  years  of  de¬ 
velopment  by  E.  Shotland,  J.  W.  Follin,  Jr.,  and  F. 
C.  Paddison  (APL)  and  J.  W.  Crispin,  Jr.,  A.  L. 
Maffett,  and  K.  M.  Siegel  (the  University  of 
Michigan).  The  results  are  in  the  form  of  first-  and 
higher-order  statistics  for  the  radar  cross  section 
(RCS)  of  a  complex  target  as  a  function  of  fre¬ 
quency,  polarization,  aspect  angle,  and  constituent 
parameters  of  the  target.  The  statistics  are  developed 
from  the  premise  that  a  complex  target  can  be  broken 
down  into  a  set  of  simple  scatterers  whose  RCS  can 
be  more  manageably  calculated  from  available  tech¬ 
niques  (thus,  the  program  can  be  conveniently  up¬ 
dated  to  include  the  most  accurate  techniques  cur¬ 
rently  at  hand)  and  whose  relative  phases  are  ran¬ 
domly  (i.e.,  uniformly)  and  independently  dis¬ 
tributed. 

Methods  for  determining  both  first-order  RCS  sta¬ 
tistics  for  a  complex  target  and  also  simple  com¬ 
ponent  scatterer  RCS  data  were  developed  by  mem¬ 
bers  of  Siegel’s  Radiation  Laboratory  at  the  Univer¬ 
sity  of  Michigan,  principally  Crispin  and  Maffett. 
The  extensions  of  these  methods  to  include  the 
(higher-order)  lobe  and  glint  statistics  were  developed 
by  Shotland  and  Follin  under  the  sponsorship  of  the 
Advanced  Research  Projects  Agency  in  an  extension 
of  the  Advanced  ALBIS  Program.  The  (higher- 
order)  beta  statistics  were  developed  by  Follin  during 
the  Cruise  Missiles  Observables  Program.  All  of 
these  statistical  methods  were  organized  into  a  coher¬ 
ent  program  by  Maffett,  and  the  result  was  coded  by 
H.  W.  Klimach.  Paddison  served  as  Program 
Manager  for  the  APL  efforts. 

The  techniques  upon  which  the  simple  scatterer 
RCS  calculations  are  based  are,  to  a  large  extent,  the 
methods  of  geometrical  and  physical  optics  with  their 
attendant  assumptions  and  restrictions  on  body 
smoothness  and  size  relative  to  illuminating  wave¬ 
length.  Included  in  the  code,  however,  are  other 
methods  that  are  appropriate,  for  example,  to  travel¬ 
ing-wave  phenomena  and  wedges  (with  straight 
edges). 

For  regions  where  the  ratio  of  body  length  to  wave¬ 
length  may  become  important  (particularly  if  radar 
absorbing  materials  have  been  used  to  reduce  large 
RCS  contributions),  diffraction  effects  have  not  been 


taken  into  account,  but  several  methods  are  being 
evaluated  to  determine  which  would  be  most  appro¬ 
priate  to  the  nature  of  the  APL  code.  The  methods 
include  the  geometrical  and  physical  theories  of 
diffraction,1-2  the  equivalent  current  technique,3  and 
the  numerical  electromagnetic  code.4  The 
geometrical  theory  is  based  on  the  tracing  of  rays  and 
can  be  described  by  currents  induced  on  the  illum¬ 
inated  portion  of  the  structure.  The  physical  theory 
adds  currents  induced  in  the  shadowed  area  by  dif¬ 
fraction  of  the  incident  field.  The  equivalent  current 
technique  adds  edge  currents  on  assumed  filamentary 
edges.  The  numerical  electromagnetic  code  is  an 
integral-equation  surface-current  determination  in 
which  the  structure  is  broken  up  into  small  cells,  and 
the  induced  charge  and  current  on  each  piece  are 
calculated  from  the  incident  field  and  from  the 
charge  and  current  of  each  other  cell.  It  is  planned  to 
include  one  or  more  of  these  methods  in  the  APL 
code  to  account  for  various  diffraction  effects  caused 
by  some  of  the  scattering  components  of  a  complex 
target. 

Bistatic  RCS  computational  capabilities  are  not  yet 
fully  incorporated  into  the  code.  The  complexity  of 
bistatic  RCS  estimation  over  monstatic  is  greatly  in¬ 
creased  by  its  additional  aspect  variables,  so  the  sta¬ 
tistical  approach  gains  still  more  importance  as  a  re¬ 
ducer  of  data  bulk.  The  bistatic  capability  will  be  in¬ 
cluded  in  the  APL  code  as  soon  as  possible. 

New  coordinate  system  arrangements  have  been  in¬ 
troduced  recently  to  simplify  the  execution  of 
conical-aspect  views  of  a  target  (see  the  Appendix). 


'.I.  B.  Keller,  "(icoim.-mi.ai  theory  of  Diffraction,''  J.  Opt. 
S ot.  Am.,  52.  I  16-140(1962). 

•|*.  fa.  Ufrtniises ,  "Method  ol  Edge  Waves  m  the  Physical 
theory  ol  Diffraction.-  /a/  lti.  Sot'.  Hudt o.  I -243  (l%2); 
translation  hy  IJ.S  Ail  force  foreign  technology  Division 
<IV7I|. 

'f.  I  Knon  and  T.  II.  A.  Senior,  "Comparison  .it  Three  Higli- 
I  resiliency  Diffraction  I  celling  tics,"  Pro*  ,  /fcfcfc  62.  1468- 
1474(19741. 

*Ci.  .1.  Burke  and  A.  J.  Poggio,  Sumeneat  Lltunmia/tnelti 
(  i  it  If  -  Method  ol  Moments,  Naval  Electronics  Systems  Com¬ 
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2.0  GENERAL  DISCUSSION  OF  RCS  DATA 


Experimental  and  predicted  RCS  data  usually  are 
in  the  form  of  a  received  signal  referenced  to  the 
transmitted  signal,  versus  aspect  angle.  Typically, 
there  are  separate  tabulations  for  each  frequency  of 
interest.  The  radars  used  for  aircraft,  missile,  and 
space-object  detection,  tracking,  and  identification 
are  usually  at  wavelengths  that  are  small  in  compari¬ 
son  to  the  dimensions  of  the  targets.  Usually  the 
targets  are  complex  shapes.  At  some  aspects  such  as 
broadside,  the  RCS  of  a  target  may  be  many  orders 
of  magnitude  greater  than  at  other  aspects,  may  be 
dominated  by  a  single  component,  and  hence  can  be 
characterized  simply.  However,  at  most  aspects  there 
is  a  complicated  interference  pattern  with  maxima 
spaced  roughly  \/L  apart  in  angle  and  c/L  apart  in 
frequency  (where  X  is  the  wavelength,  c  is  the  velocity 
of  light,  and  L  is  a  typical  length  parameter).  If,  for 
example,  L  is  10  meters  (i.e.  a  small  missile),  roughly 
109  data  points  are  needed  to  cover  X  band  for  ail 
polarizations.  A  statistical  description  is  needed  to 
reduce  the  data  storage  and  retrieval  problem. 


The  magnitude  of  RCS  data  is  not  the  only  reason 
for  an  interest  in  a  statistical  description.  The  main 
users  of  these  data  are  the  hardware  designer,  the 
radar  system  designer,  and  the  performance  analyst. 
The  latter  two  have  long  searched  for  a  statistical  des¬ 
cription  of  the  amplitude  scintillation  of  RCS.  The 
literature  contains  many  attempted  fits  of  experi¬ 
mental  data  with  single-parameter  statistical  des¬ 
criptors. 

The  parameters  of  interest  in  the  vicinity  of  any 
one  aspect  are  the  mean  RCS  and  its  probability  dis¬ 
tribution,  the  lobe  widths  in  angle  and  frequency, 
and  the  mean  centroid  and  its  variances. 

Over  the  years,  the  technique  of  breaking  down  a 
complex  target  into  a  finite  set  of  component  scatter¬ 
ed  has  been  evolved  together  with  statistical  descrip¬ 
tions  of  lobe  structure  and  the  appropriate  descriptor 
of  amplitude  scintillation,  the  two-parameter  beta 
distribution  function.  The  computer  code  described 
in  this  report  embodies  the  aforementioned  predic¬ 
tive  RCS  and  its  statistical  descriptors. 
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3.0  THE  CAPABILITIES  OF  THE  CODE  FROM  THE  VIEWPOINT 

OF  THE  RCS  ANALYST 


The  APL  RCS/Statistics  Code,  designed  for  the 
PDP  11/60  computer,  has  been  checked  out  with 
several  missile  configurations.  Its  computational 
capabilities  fall  into  two  categories,  either  of  which, 
with  various  options,  can  be  executed  from  a  master 
code.  They  are: 

1 .  First-order  statistics  for  RCS,  and 

2.  Lobe  and  higher-order  statistics  for  RCS. 


The  entire  development  is  based  on  the  premise 
that  a  complex  scattering  body  can  be  broken  down 
into  a  finite  set  of  component  scatterers  (usually 
chosen  to  be  elementary  scatterers,  some  of  which 
may  correspond  to  simple  geometric  shapes)/  If  the 
contributions  from  this  set  of  scatterers  are  combined 
with  proper  phases,  the  result  is  an  estimate  of  rela¬ 
tive-phase  RCS;  this  option  is  available  in  the  APL 
code.  Figure  1  is  an  example  of  relative-phase  RCS. 


V> 

a 
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At  this  stage,  the  first-  and  higher-order  statistics  are 
limited  to  monostatic  situations.  However,  the  lobe 
statistics  are  coded  in  such  a  way  as  to  be  applicable 
to  bistatic  situations  as  soon  as  that  capability  is  in¬ 
troduced. 


J.  W,  Crispin.  Jr.,  ami  A.  1..  Mai  ten .  "Esiiinaiiou  of  RCS  for 
Simple  Sha|ies."  Isiimanon  of  RCS  for  Complex 
Shapes,"  Prm.  UAL  53.  <  I%5|. 
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Frequency  17.96  GHz 
Horizontal  polarization 


20 


Pitch  0 
Roll  0 
“  Yaw  0 
Tilt  0 


[_  (port) 


Observation  angle  (deg) 


Figure  1  —  Example  of  calculated  monostatic  relative- 
phase  RCS  1/24th  scale  Convair  990. 


If  it  is  assumed  that  phases  among  scattering  com¬ 
ponents  are  distributed  independently  and  randomly 
(i.e.,  uniformly,  over  the  interval  0  to  2x),  a 
hierarchy  of  statistics  can  be  developed.  Under  first- 
order  statistics  fall  the  usual  mean  RCS  and  root- 
mean-square  (rms)  spread  about  the  mean  (see  Fig. 
2).  Under  higher-order  statistics  fall  higher  moments 
from  which  skewness  and  kurtosis  may  be  derived  to 
show  the  appropriate  beta  distribution  characteristics 
for  the  RCS  data  under  examination.'’*  (These 
higher-order  statistics  can  also  be  used  to  analyze 


6J.  W.  Follin,  Jr.,  Statistical  Pro/terties  of  RCS.  J  HU/  API  (JM- 

81-115(1981). 

J.  W.  Follin,  Jr.,  and  A.  L.  Maffell,  "RCS  Scintillations  and 
Their  Statistical  Description,"  in  Proc.  Second  Annual  lac- 
tical  Air  Surveillance  Control  Con/.,  Rome  Air  Development 

Center  (1981). 

*J.  W.  Follin,  Jr.,  F.  C.  Paddison,  and  A.  L.  Vtaffcn,  “  The  RCS 
of  Two  Cruise  Missiles  and  Their  Statistical  Description," 
DARPA  8th  Strategic  Space  Symposium  ( 1982). 
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64.1%  shadowed 
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Figure  2  —  Example  of  calculated  monostatic  mean  RCS 
with  rms  spread  (including  lobe  statistics);  1/25th  scale 
Convair  990. 


measured  RCS  data;  this  feature  is  not  part  of  the 
APL  code.)  Various  properties  of  the  autocorrelation 
of  RCS  can  be  used  to  describe  lobe  statistics  (in  both 
frequency  and  aspect).1'  Lobe  characteristics  are  illus¬ 
trated  in  Fig.  2  at  observation  angle  intervals  of  30°. 
The  terms  “major"  and  “minor”  refer  to  half  the 
long  and  short  axes  of  the  ellipse;  they  give  the  mean 
lobewidth  at  that  aspect  in  both  the  horizonal  pL..c 
(minor)  and  the  normal  plane  (major).  “Angle” 
refers  to  the  tilt  of  the  lobe  ellipse  with  the 
observation  axis.  N,  (GHz)  refers  to  the  number  of 
lobes  of  RCS  per  gigahertz  of  frequency. 

At  this  writing,  the  component  scatterer  analyses 
rest  largely  upon  physical  and  geometrical  optics 
methods.  However,  the  code  could  be  expanded 
readily  to  include  diffraction  treatments  and  even  the 
more  exact  integral-equation  analyses  exemplified  by 
the  numerical  electromagnetic  code.4 

The  program  as  written  and  described  here  was  in¬ 
tended  for  targets  of  modest  complexity  at  observing 
frequencies  in  the  centimeter  region.  We  plan  to  re¬ 
code  the  program  for  a  more  capable  computer  in  the 


T.  .Midi land.  Statistical  Analysts  »/  Radar  I uriici  Scintillation. 
JHU/APl.  rti 955  (19671. 
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near  future  to  allow  the  handling  of  more  complex 
targets  up  to  the  mid-millimeter  range  of 
frequencies. 

The  effects  of  a  surface  treatment  with  radar  ab¬ 
sorbing  material  can  be  taken  into  account  through 
the  uniform  application  of  a  reflection  coefficient  to 
the  component  on  which  the  material  is  placed.  No 
account  is  taken  of  the  possibility  of  volumetric  scat¬ 
tering  for  either  radar  absorbing  material  or  dielec¬ 
tric  components. 

The  major  emphasis  of  the  code  is  on  a  statistical 
examination  of  the  RCS  of  a  target.  Considering  the 
number  of  variables  upon  which  RCS  depends,  a 
very  voluminous  collection  of  data  can  be  amassed 
for  a  single  target,  even  when  domains  of  variables 
are  rather  severely  restricted.  A  statistical  examina¬ 
tion  of  RCS  offers  a  condensation  of  facts  and,  thus, 
an  economy  of  output  not  available  with  a  determin¬ 
istic  treatment  of  similar  scope. 

All  of  the  RCS  statistics  are  essentially  unaffected 
by  small  changes  in  target  configuration  or  frequency 
that  may  drastically  affect  the  lobe  structure  of  deter¬ 
ministic  treatments  (especially  for  frequencies  of  1 
GHz  and  above). 

Last  but  not  least,  statistical  RCS  results  mesh  well 
with  the  requirements  of  detection  and  tracking  sys¬ 
tems  for  statistical  formulations. 


3.1  Inputs 

The  inputs  usually  required  for  the  code  are  wave¬ 
length;  polarization  (linear,  but  extendable  to  ellip¬ 
tical);  plane  of  observation;  range  of  observation 
angle  in  the  observation  plane;  pitch,  roll,  and  yaw  of 
the  target;  and  tilt  of  the  target  rotation  plane  (to 
make  comparisons  with  RATSCAT  measurements). 

For  a  particular  complex  target,  a  suitable  mathe¬ 
matical  model  must  be  derived  in  the  form  of  a  set  of 
component  scatterers  matched  to  the  target  features 
with  the  help  of  drawings,  pictures,  and  artists’  con¬ 
ceptions.  The  set  of  component  scatterers  must,  of 
course,  be  limited  to  the  types  presently  available  to 
the  APL  code  (some  20  as  of  this  writing). 

Computations  of  RCS  are  performed  in  coordinate 
systems  generic  to  the  various  components  and  so 
must  be  transformed  to  the  coordinate  systems  in 
which  the  statistics  are  formulated. 


3.2  Outputs 

The  first-order  and  lobe  statistics  have  already 
been  illustrated  in  Fig.  2.  Note  that  for  the  latter 
case,  for  any  line  through  the  lobe  ellipse  center,  the 


distance  from  the  center  to  the  ellipse  intercept  gives 
the  lobewidth  for  that  direction  of  observation.  For 
example,  at  the  30°  observation  angle,  the  ellipse  has 
a  major  axis  oriented  82°  from  the  observation  axis, 
indicating  thin  lobes  in  the  horizontal  plane  of  the 
target  and  fat  lobes  in  the  perpendicular  plane. 

The  output  for  the  higher-order  statistics  may  be 
either  printed  or  graphical.  Let  us  consider  the  latter. 
One  has  available  skewness  versus  width  or  kurtosis 
versus  width  (width,  skewness,  and  kurtosis  are  de¬ 
fined  in  Section  4  of  this  guide).  Skewness  and 
kurtosis  are  plotted  in  Figs.  3  and  4,  respectively,  ver- 


Width 


Figure  3  —  Calculated  width-skewness  statistics  for  the 
RCS  scintillation  of  a  l/25th  scale  Convair  990. 


sus  width  for  a  typical  target  at  10°  observation  angle 
intervals.  It  should  be  noted  that  the  data  presented 
in  Figs.  3  and  4  may  be  based  either  on  experimental 
data  or  on  theoretical  calculations.  The  code  can  pro¬ 
vide  the  theoretical  values.  The  points  are  plotted 
within  boundaries  that  denote  the  limits  for  beta-dis¬ 
tributed  component-RCS  values.  Locations  denoting 
normal,  lognormal,  chi-square,  and  Rayleigh  distri- 
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Figure  4  —  Calculated  width-kurtosls  statistics  for  the  RCS 
scintillation  of  a  1/25th  scale  Convair990. 


butions  all  fall  outside  of  these  boundaries,  for  both 
width/skewness  space  and  width/kurtosis  space. 

For  example,  consider  Fig.  4  where  the  observation 
direction  100*  produces  in  width/kurtosis  space  the 
point  (0.7, 0.7).  These  statistics  yield  the  parameter 
values 


a  =  0.96,  b  =  3.82 


for  the  beta  distribution  given  by 


&  (x;a,b)  = 


T(a  +  b) 
V(a)T(b) 


jf-'  (1 


-  x 


where  T(a)  denotes  the  gamma  function  of  an  argu¬ 
ment  a.  Thus,  the  component  scatterers  that  con¬ 
tribute  to  the  RCS  at  100*  from  nose-on  produce, 
under  the  radom-phase  assumption,  an  RCS  that  is 
distributed  according  to  the  above  beta  distribution. 
[Note:  The  skewness  can  be  calculated  to  be  0.88, 
which  agrees  very  well  with  the  skewness  for  100* 
found  at  width  0.7  in  Fig.  3.  The  beta  distribution  for 
the  above  values  of  a  and  b  is  skewed  to  the  left,  with 
a  long  tail  to  the  right.  ] 

The  beta-distributed  .quantity  x  is  actually  o/omax, 
where  a  is  the  RCS  for  some  fixed  Aspect,  frequency, 
and  polarization.  So,  for  example,  since  amttX  (100°) 
=  1.14  m2  and  since  the  mean  of  &(x;a,b)  is 
a/(a+b)  =  0.2,  the  mean  RCS  at  100*  is  a  = 
0.2(1.14)  =  0.23  m2  or  -6.4  dBsm  (dBsm  is  decibels 
relative  to  a  square  meter).  Note  that  this  value 
agrees  with  the  mean  RCS  at  100*,  which  can  be  read 
on  Fig.  2.  The  beta  distribution  of  RCS  at  a  fixed 
observation  angle  arises  from,  and  depends  on,  the 
assumption  that  the  phases  among  components  are 
randomly  distributed  (i.e.,  uniformly,  between  0  and 
2w)  for  every  observation  angle  value. 


3.3  Glint 

Illustrations  of  the  output  of  the  program  predict¬ 
ing  glint  are  not  included  in  this  writing;  that  portion 
of  the  program  has  yet  to  be  validated.  However,  a 
discussion  of  the  theory  used  is  given  in  Section 
4.6. 


j 
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4.0  CALCULATIONS 


4.1  Relative-Phase  RCS 

Consider  a  class  of  complex  targets,  each  of  whose 
members  can  be  broken  down  into  a  finite  set  of 
component  scatterers.  Each  component  scatterer  is 
assumed  to  be  in  the  far  field  of  the  transmitter  so 
that  it  is  illuminated  by  a  plane  wave;  it  is  also  in  the 
far  field  of  the  receiver  so  that  the  observed  field  has 
a  plane  wavefront.  The  major  dimension  of  a  com¬ 
ponent  is  assumed  to  be  large  with  respect  to  the 
wavelength  X.  Polarizations  of  the  transmitter  and 
receiver  are  taken  to  be  horizontal  (H)  or  vertical  (V) 
in  this  discussion;  however,  they  could  be  circular  or 
elliptical.  The  possible  combinations  are  HH,  HV, 
VH,  and  W;  only  the  first  and  last  are  considered 
here. 

Let  the  signal  scattered  by  the  tth  component  be 
denoted  by 


where  the  asterisk  means  complex  conjugate  and 
<  •  •  •  >  indicates  an  average  over  all  phases.  The  de¬ 
viation  from  the  mean  is 

F=  AA‘  -  E[a 1 

=  2  £  a,ak  cos  (i,  -  ik ) 


=  2  Yi  cos (i,  -  ) 


whence  the  variance  is 


=  </*>  =  2S2  =  2  £  a,ak 


a,exp(ji,),  j  =  V- 1 ",  (1) 

where  the  amplitude  ai  (which  is  usually  taken  to  be 
where  a,  is  the  RCS  of  the  rth  component)  and 
the  phase  i,  may  depend  on  aspect  (for  both  trans¬ 
mitter  and  receiver,  or  on  their  bistatic  separation  0), 
frequency,  polarization,  and  the  constituent  param¬ 
eters  of  the  component.  The  scattered  amplitude  is 
then 

^  Vff ,  exp in,).  (2) 


-(£«.)-  E  «.’• 

/=  1  1=  I 

4.3  Higher-Order  Statistics 

If  the  triple  and  quadruple  sums  are  written  as 

s3  =  £ 


The  relative-phase  RCS,  a,  for  the  complex  target  is 


a  -  A  A*  =  I  £  Vff,  exp  CM) 


=  £  'fok  expOV, )  exp(  -jik )  . 


S4  —  ^  o,  ok  o,am  , 


the  third  and  fourth  central  moments  are,  respective¬ 
ly.6 


4.2  Random-Phase  RCS 

For  the  random-phase  approximation,  the  mean 
RCS,  E[a],  over  all  random  phases,  is 

v,  !-£[*])  -  <AA*>  =  E<7„  (4) 


=  <£*>  =  12S, 


Vk  =  (F*)  =  6 S22  +  12S,S3  +  13254  .  (10) 


v 


v 
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In  order  to  examine  the  nature  of  the  distributions 
governing  RCS  for  a  complex  target,  we  have 
modified  Pearson’s  method  somewhat.  Instead  of 
examining  target  RCS  behavior  in  skewness-kurtosis 
space,  we  have  set  up  two  spaces,  one  a  width-skew¬ 
ness  space,  the  other  a  width-kurtosis  space,  where 
width  w  is  defined  as 

w=/j,/»V.  (11) 

In  fact,  to  simplify  the  analysis  somewhat,  we  have 
defined  a  modified  skewness  as 

Yi  =  Pj'V’i'i  •  (12) 

Kurtosis  is  defined  in  the  standard  manner  as 

>2  =  -  3  .  (13) 

M  2 

From  the  nature  of  the  sums  S„  /  =  1,  2,  3,  4, 
above,  it  can  be  shown  that  the  allowable  regions  in 
width-skewness  (actually,  read  modified  skewness) 
and  width-kurtosis  spaces  are  almost  triangular  and 
are  restricted  as  follows  (see  Figs.  3  and  4): 


4.3.1  Width-Skewness 
Right  boundary: 


y;  =  4w  -  2,  -  <  w  <  1  ; 


where 


7,'  +  2  .  ,  . 

n  =  - - -  and>,  is  as  above, 

2w  -  7| 


<  w<  1  ,  lim  y:  =  6  ; 


Lower  boundary: 


>  =  -  1.5  ,  0<  w<  -  ; 

2 

Left  boundary: 

part  1:  w  =  0 ,  -  1.5  <7,  <0; 


part  2:  w  = 


/(/  +  2) 
(/+  D2 


(as  above) , 


6  t(t  +  4) 

*  =  •  (19b) 

The  location  of  points  within  the  restricted  regions 
will  vary  with  the  aspect  at  which  a  complex  target  is 
illuminated  and  viewed.  The  corresponding  beta  dis¬ 
tribution  parameters  may  be  calculated  from  those 
coordinates.  The  beta  distribution  is  defined  as 


Lower  boundary: 


0U)  = 


T(o  +  b ) 
r(o)f(b) 


jf“- 1  (1  -  x)' 


7.'  =  0,  o  <  w  <  -  ; 


Left  boundary: 

t(t  +  2) 
w  = - -  , 

( t  +  l)2 


2 /(/  +  3) 

( t  +  IK/  +  2)  ’ 


Ostsoo. 


4.3.2  Width-Knrtosis 
Right  boundary: 


6n2  13  1 

y*  ~  (n+  lKn  +  2)  W  +  n2  w  n  n2  ]  ’  0  ) 


0  <  x  <  1  , 

where  T(a)  is  the  gamma  function  of  argument  a, 
and  x  =  .  Then,  for  example,  the  parameters  a 

and  b  can  be  determined  from  width-kurtosis  space. 
The  result  can  be  used  to  find  skewness,  which  can  be 
checked  against  the  location  of  the  aspect  point  in 
width-skewness  space.  The  equations  are 


a(a  +  b  +  1) 
2 (b  -  a) 


a(a  +  b  +  2)  ’ 

6(o  +  b  +  I) 

”  (a  +  b  +  2K»  4-  *  +  3) 

/a  b  1  \ 


-WWW WWKK! 


(21) 
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The  values  of  a  and  b  may  be  determined  from  any 
pair  of  these  equations,  but  the  first  two  should  be 
used  for  experimental  data.  If  one  takes,  for 
instance,  the  aspect  of  100°  (Fig.  2),  one  can  compute 
a  =  0.96  and  b  =  3.82. 

In  order  to  check  the  beta  distribution  theory 
against  measured  data,  a  slightly  different  viewpoint 
must  be  taken  because  the  output  of  a  measurement 
exercise  is  usually  in  the  form  of  RCS  versus  aspect. 
If  the  RCS  is  sufficiently  oscillatory  in  the  chosen 
aspect  interval  (usually  10°,  but  ranging  from  5  to 
30°,  depending  on  frequency),  the  random  phases 
have  had  the  opportunity  to  range  over  most  of  the 
intervals  of  definition,  and  so  the  RCS  in  the  aspect 
interval  should  be  distributed  in  essentially  the  same 
way  that  the  RCS  at  a  fixed  aspect  was  owing  to  the 
random  phase  fluctuations. 


4.4  Modified  Statistical  Distribution 

An  extension  of  the  statistical  theory  just  presented 
has  been  developed  where  the  requirement  of  con¬ 
stant  component  RCS  is  relaxed  and  an  individual 
variance  and  third  and  fourth  moments  about  the 
mean  are  assigned.  The  physical  regions  of  the  statis¬ 
tical  space  resulting  in  this  modification  will  be  pub¬ 
lished  shortly. 


4.5  RCS  Lobe  Structure 

If,  instead  of  taking  the  relative  phases  to  be  ran¬ 
dom  variables,  one  examines  their  variation  (and 
hence  the  variation  of  relative-phase  RCS)  as  a  func¬ 
tion  of  receiver  aspect  (in  a  bistatic  situation,  either 
aspect  or  bistatic  angle),  it  is  possible  to  obtain 
scintillation  statistics.  The  results  are  derived  from 
the  RCS  correlation  function,  using  the  zero  crossing 
theorem  of  Rice.10  They  represent  an  extension  of 
results  for  the  monostatic  case  obtained  by 
Shetland.9 

Suppose  the  target  is  at  the  origin  of  a  £,  tj,  ("coor¬ 
dinate  system  whose  £f  plane  contains  the  transmitter 
and  receiver.  Then  the  phase  term  4>,  of  Eq.  1  is  given 
by 


-  *[f/(l  +  COS 0)  +  £,  sin  0)  ,  (22) 


l0S.  O.  Rice,  "MaihemaliL'al  Analysis  ol  Random  Noise,"  in 
Selected  Papers  on  ISoise  and  Stochastic  Processes,  Nelson 
Wax,  cd.,  IJover  Publieanons,  New  York  ( I V54). 


where  A:  is  the  wavenumber.  The  results  obtained  are 
the  number,  N,,  of  lobes  of  RCS  per  hertz  of  fre¬ 
quency: 


N,  = 


1 

cS 


<>,  ",  1  ( f.  -  U  )  ( 1  +  COS  ,i) 


+  (£,  -  £, )  sin  dl") 


(23) 


where  c  is  the  speed  of  light  and  S  =  \p:  is  from  Eq. 
6.  Also  obtained  are  the  number,  N  ,  of  lobes  of  RCS 
per  radian  in  traverse  relative  to  a  rotation  axis 
oriented  at  an  angle  6  with  the  £  axis  in  the  £f  plane: 

N  =^s(LAd,d)  .  r,s  =1,2,  (24) 


where 

d  i  =  sin6,  d,  =  -cosiS,  (25) 


and 


l<£,  -  £, )  <  1  +  cos  if) 

i  ■  k 


-  (C,  -  f, )  sin  d)4' ' 


x  l(t  -»/,)( 1  +  cosd)]"'  :  •  (26) 


The  numbers  /V,  and  N  ,  considered  as  a  function  of 
f>  for  a  fixed  aspect  angle,  form  ellipses.  They  can 
also  be  converted  into  local  lobewidths  of  frequency 
and  aspect,  respectively,  by  considering  their  recipro¬ 
cals.  For  any  line  through  the  ellipse  center,  the  dis¬ 
tance  from  the  center  to  the  ellipse  intercept  gives  the 
lobe  rate  for  that  direction  of  rotation  axis. 


4.6  Glint 

The  apparent  electrical  tocation  of  the  target  re¬ 
flection  depends  on  the  type  of  radar  seeker  used  to 
measure  it.  For  the  relative-phase  return  (Eq.  2)  with 
the  receiver  location  at  X,  Y,Z;  X,  Y<  < Z;  we  have 


d*i  =  1  {kr  k{X  ~  il  =  1 

dX  dx  r)  r,  r 


(27) 
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and  the  variance  of  the  glint  is 


W,  k(Z-t)  ^ 
dZ  r,  ' 


<A£2>  =  ^  E  W  +  «>  "  •  <32> 


The  apparent  source  location  is 


i  =  r tan'  ( - — —  ) 
^dA/dZ  ' 


So  far  we  have  discussed  only  one  component  of 
glint,  but  it  is  clear  that  expressions  similar  to  Eqs.  30 
and  32  are  needed  for  ??.  In  addition,  the  variance  is  a 
tensor  operator  so  that  an  off-diagonal  term  is 
needed. 

Defining 


S  ika,  £,  expQY', ) 
Lika,  expOY,-) 


B*.  =  E  +  ij  -  2|r 


E  «(«/{/ expo'll  -  *,  )) 


X  (,,  +  ijy  -  2)j)m+"-2 ,  m,/i  =1,2  (33) 


^  a,aj  exp(i('h  -  tj)) 


*  =  E  - 


^  E  «/«;*/  W.  “  */>. 

£  u 


where  a,  =  cos  a,  a2  =  sin  a,  the  standard  deviation 
of  glint  in  the  a  direction  is 


where  the  real  part  is  taken  in  the  last  step  since  the 
rf  signal  is  actually  proportional  to  the  real  part  of 
Eq.2. 

The  di  ision  is  normally  performed  by  an  auto¬ 
matic  gain  control  (AOQ  circuit  with  a  time  constant 
so  that  the  numerator  and  denominator  are  averaged 
independently.  The  result  in  the  random-phase  ap¬ 
proximation  is 


s--2£F5'®- 


Further  in  analogy  to  the  derivation  of  Eqs.  23  and 
26,  we  define 


cm„  =  e  wMi  +  tj  -  w-" 


X  (Vi  +  Vj  -  2ii)n 


With  a  fast  AGC,  the  fluctuations  of  Eq.  29  become 
larger  but  cannot  be  expressed  in  closed  form.  In 
addition,  the  result  depends  on  the  dynamic  range  of 
the  circuits.  Note  that  in  a  bistatic  geometry  the  glint 
is  measured  relative  to  the  receiver  and  is  indepen¬ 
dent  of  bistatic  angle  except  for  the  effects  of 
shadowing  of  components. 

Subtracting  Eq.  30  from  Eq.  29,  we  have  the  devi¬ 
ation  from  the  mean  for  the  slow  AGC  case. 


X  [(£  -  fy)(l  +  COS0)  +  ({,  -  £,)  sin0]2 , 


Dmnn  -  E  *i*y«i  +  -  2i)4 


x  (,,  +  ij,  -  2,>"+"-2 


x  [(£,  -  l,)(l  +  cos/3)  -  (f,  -  ft)  sin/3]4 


E  ~  Dcos(*,  “  tj) 


X  t(n/  -  %)( 1  +  cos/3)]' 


■  1  E  aiajUi  +  *y  “  2|)cosM  -  ij) .  (31) 

®  i<j 


c=  E  » 


v  (  •  A 


v.v.v 
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and 


D=  £  Dmnnama„drds.  (39) 


Then  the  number  of  glint  lobes  in  the  a  direction  per 
cycle  is 


N, 


°J  ~  2 c5„  ^ 


(40) 


and  the  number  of  glint  lobes  in  the  a  direction  per 
radian  of  rotation  about  the  axis  oriented  in  the  6 
direction  is 


(41) 


Note  that  Eqs.  40  and  41  are  correct  even  though 
the  assumption  of  a  slow  AGC  has  been  used;  only 
the  amount  of  the  displacement  is  in  error. 


5.0  USER’S  GUIDE  TO  THE  CODE  FROM  THE  VIEWPOINT 
OF  THE  COMPUTER  ANALYST 


5.1  Introduction 

The  APL  RCS/Statistics  Code  is  a  computer  soft¬ 
ware  package  capable  of  computing  and  displaying 
the  collective  cross  section  from  a  group  of  simple 
components  using  physical  optics  as  the  main  model¬ 
ing  basis.  In  addition  to  average  cross  section,  it  is 
also  able  to  compute  and  display  relative-phase  cross 
section,  RCS  lobe  and  glint  lobe  characteristics,  and 
higher-order  beta  statistics  of  the  collective  object  as 
a  function  of  observation  angle,  in  either  printed  or 
graphical  form.  Input  for  the  program  is  taken  from 
a  data  base  in  batch  mode  process,  and  output  results 
are  displayed  on  a  graphics  terminal,  plotter,  or  line 
printer.  The  RCS  analytical  basis  is  given  in  Refs.  1 1 
and  12. 

The  program  is  currently  running  on  a  PDP  11/60 
computer  using  the  RSX-1 1M  V3.2  operating  system. 
The  graphical  peripherals  available  on  this  system  are 
the  Tektronix  terminal  and  the  ZETA  plotter.  The 
program  is  executed  by  issuing  the  computer  com¬ 
mand  RUN  RCS.  It  responds  by  prompting  the  user 
for  the  name  of  the  input  data  base  file.  The  data 


"A.  L.  Maffett,  Bistatic  Lobe  Statistics,  JHU/APL  M80-7  (Nov 
19W). 

A.  L.  Maffett,  Linear  Transformations  Relating  Theoretical  and 
Experimental  Coordinate  Systems  Used  in  RCS  Determina¬ 
tion,  JHU/APL  QM-80-049  (M804)  (Feb  1980). 


base  name  must  then  be  entered,  followed  by  any 
optional  switches.  The  available  options  are: 

/G  enable  graphical  debugging 
/GC  enable  graphical  debugging 
(combined  components) 

/P  enable  printed  debugging 

/Z  direct  graphics  output  to  plotter 

(default  is  terminal) 

The  program  then  performs  the  functions  defmed  in 
the  data  base  and  returns  to  computer  command 
mode. 

This  section  of  the  report  provides  the  user  with  a 
basic  understanding  of  the  program  along  with  a 
comprehensive  reference  manual  for  actual  program 
operations.  The  internal  structure  and  program  capa¬ 
bilities  are  presented  in  Section  S.2.  Section  3.3  pre¬ 
sents  the  display  options,  including  the  available  de¬ 
bugging  output.  General  features  of  the  input  data 
base,  such  as  looping  and  command  priorities,  are 
described  in  Section  3.4.  The  actual  RCS  commands 
are  presented  in  detail  in  Section  3.3. 


5.2  RCS  Program  Operations 

The  program  is  designed  to  compute  the  combined 
RCS  (called  random-phase  RCS)  of  several  simple 
components  as  a  function  of  observation  angle, 
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wavelength,  and  polarization.  In  addition  to  average 
cross  section,  the  program  also  computes  relative- 
phase  cross  section,  RCS  lobes,  average  glint,  and 
glint  lobes.  Finally,  it  may  also  be  used  to  compute 
higher-order  beta  statistics. 

A  complex  object  to  be  analyzed  in  terms  of  RCS 
must  first  be  separated  into  its  constituent  com¬ 
ponents.  The  basic  components  available  to  the  pro¬ 
gram,  such  as  ogives,  cylinders,  corner  reflectors, 
etc.,  are  described  in  Section  5.5.  Each  component 
must,  at  a  minimum,  be  specified  by  its  name,  physi¬ 
cal  dimensions,  and  orientation  in  the  radar  coordi¬ 
nate  system.  The  location  of  each  component  (or 
scattering  center)  in  the  radar  coordinate  system  must 
also  be  specified  if  lobe  or  glint  output  is  required. 
The  geometry  describing  the  scattering  centers  of 
each  available  component  has  been  specified  in  Ref. 
13. 

The  program  may  be  divided  functionally  into  four 
areas:  input,  cross-section  computation,  statistics 
and  lobe/glint  computations,  and  output  display. 
The  input  function,  which  is  actually  in  continuous 
operation  throughout  the  computations,  is  described 
in  detail  in  Section  5.4  in  the  context  of  a  data  base. 
In  summary,  the  input  must  be  ordered  to  supply  the 
program  with  global  parameters,  followed  by  local 
component  parameters,  and  ending  with  display  pa¬ 
rameters.  Since  knowledge  of  the  input  function  or 
output  display  function  is  not  necessary  to  under¬ 
stand  the  program’s  computational  algorithms,  they 
will  not  be  discussed  further  in  this  section. 

The  program  first  zeros  an  array  of  accumulators 
and  then  adds  in  the  cross-section  contributions  of 
each  component  read  on  input.  This  process  con¬ 
tinues  until  the  advent  of  an  output  display  instruc¬ 
tion,  which  causes  the  accumulated  results  to  be  dis¬ 
played  as  specified.  The  RCS  program  operates  com¬ 
pletely  on  one  component  at  a  time  by  computing  the 
component  RCS  at  a  fixed  wavelength,  polarization, 
and  radar  phi  angle  and  then  varying  the  radar  theta 
angle  (phi  and  theta  are  the  radar  polar  coordinates). 
Each  cross-section  value  as  a  function  of  theta  is 
stored  in  a  separate  accumulator.  The  range  and 
number  of  theta  values  desired,  specified  on  input, 
must  not  require  more  than  181  accumulators  (a  pro¬ 
gram  limitation).  This  limitation  may  be  circum¬ 
vented  in  some  instances  as  described  in  Section  5.3 
and  defined  in  Section  5.5. 

In  addition  to  accumulating  average  cross-section 
responses,  the  program  also  saves  the  square  of  the 
values  of  the  responses  in  order  to  compute  and  dis- 
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play  the  standard  deviation  of  the  overall  response  at 
each  theta  increment.  If  relative- phase  rather  than 
average  cross-section  output  is  requested,  it  will  use 
the  two  accumulator  vectors  to  save  the  complex 
components  of  cross  section  rather  than  the  detected 
sum  and  sum  squared  values.  All  results  generated 
thus  far  may  be  stored  in  computer  memory  since 
they  are  relatively  few  in  number.  The  lobe  and/or 
glint  options,  however,  require  the  temporary  storage 
of  large  volumes  of  intermediate  data.  Those  data, 
which  consist  of  the  cross  section  and  three-dimens¬ 
ional  position  of  each  subreflector  in  the  object  for 
each  theta  angle  of  interest,  are  stored  on  a  disk  file 
for  later  processing. 

The  following  paragraphs  summarize  the  RCS  pro¬ 
gram  procedure  for  computing  cross  section  for  a 
single  component.  Cross-section  calculations  do  not 
actually  begin  until  a  component  description  has 
been  read  from  the  input  stream. 

The  program  begins  to  process  a  component  by 
computing  a  transformation  matrix  from  component 
space  through  object  space  to  radar  space.  The  com- 
ponent-to-object  transformation  matrix  (TRANS) 
will  have  been  read  previously  as  input.  The  object  - 
to-radar  transformation  matrix  will  have  been  com¬ 
puted  from  previously  read  values  of  radar  platform 
pitch,  roll,  yaw,  and  tilt  (PRYT). 

The  program  then  enters  the  THETA  loop  to  com¬ 
pute  the  cross  section  at  each  required  theta  angle.  In 
the  loop,  it  determines  if  the  current  radar  theta  and 
phi  angles  are  valid  for  that  component.  If  valid,  they 
are  transformed  into  the  component  coordinate  sys¬ 
tem.  At  that  point,  calculations  that  correspond  to 
the  specific  component  of  interest  are  performed. 

Although  they  differ  in  detail,  ail  groups  of  com¬ 
ponent  calculations  (i.e.,  subroutines)  have  several 
things  in  common.  They  must  first  collect  their  ap¬ 
propriate  inputs  and  then  determine  if  a  valid 
response  is  possible  from  the  current  parameter  con¬ 
figuration.  Each  component  subroutine  must  then 
compute  its  RCS  as  a  function  of  geometry,  radar 
characteristics,  and  observation  angle.  If  lobe  output 
is  also  required,  the  component  subroutines  must 
also  compute  the  location  of  the  reflecting  source  or 
sources  for  multiple  reflectors.  Component  position 
information,  after  being  transformed  to  the  radar  co¬ 
ordinate  system,  is  stored  in  a  file  for  future  lobe 
and/or  glint  processing. 

When  the  component  cross  section  has  been  com¬ 
puted,  it  is  added  to  the  accumulator  bin  appropriate 
to  the  current  theta  angle.  The  powers  of  these  results 
are  computed  and  stored  in  another  file  if  optional 
beta  statistics  are  desired.  Beta  statistics  are  com¬ 
puted  in  the  output  display  section  of  the  program 
using  the  relatively  simple  equations  from  Ref.  7.  In 
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any  case,  the  observation  angle  is  incremented  by  by  the  angle  delta;  in  contrast,  RCS  and  glint 

repeating  the  THETA  loop.  After  the  component  moments  are  defined  by  the  observation  angle  theta, 

cross  section  has  been  computed  for  all  required  When  computations  for  all  values  of  delta  have  been 

values  of  theta,  program  execution  returns  to  the  in-  completed,  the  results  are  displayed,  and  execution 

put  section  to  obtain  new  component  information  control  is  returned  to  the  main  program, 

from  the  data  base,  and  the  process  is  repeated. 

Basic  RCS  computations  pause  when  a  display  in¬ 
struction  is  read  on  input.  This  causes  the  program  to 

display  its  cross-section  results  (described  in  Section  5.3  Display  Options 


S.3)  and  then  begin  lobe  and/or  glint  calculations  if 
required.  The  remainder  of  this  section  highlights  the 
general  procedure  for  computing  lobe  and  glint  char¬ 
acteristics.  A  complete  analytical  description  is  given 
in  Refs.  14  and  15. 

The  RCS  lobe  calculations  begin  by  reading  tem¬ 
porary  results  stored  by  the  main  program.  Each 
component’s  contribution  to  RCS  at  the  observation 
angle  of  interest  (specified  by  DELTA,  a  subset  of 
THETA)  is  read  by  the  lobe  subprogram  along  with 
its  position  in  space.  The  first  parameters  computed 
from  these  data  are  a  simple  cross-section  sum  and 
standard  deviation.  Cross-section  products  of  com¬ 
ponent  cross  section  and  position  are  then  computed 
and  summed.  The  square  root  of  this  result  divided 
by  the  standard  deviation  and  the  speed  of  light 
yields  the  number  of  lobes  of  radar  cross  section  per 
hertz  of  frequency.  A  different  set  of  summed  cross 
products  is  used  to  compute  lobe  width. 

The  glint  calculations  begin  by  computing  average 
moments  for  the  components  read  in  the  lobe  sec¬ 
tion.  Sums  of  cross  products  using  component  cross 
section,  position,  average  moments,  and  an  intersect 
plane  defined  by  the  angle  alpha  are  used  to  compute 
the  glint  standard  deviation.  The  average  moment 
arm  length,  along  with  standard  deviation,  as  a  func¬ 
tion  of  observation  angle  (theta)  is  the  primary  out¬ 
put  from  the  glint  section. 

Glint  lobes  are  computed  from  the  same  set  of  data 
used  to  compute  glint  moments.  In  this  case,  an  ex¬ 
panded  set  of  summed  cross  products  is  computed 
and  used  in  conjunction  with  glint  standard  deviation 
to  yield  the  number  of  glint  lobes  per  hertz  of  fre¬ 
quency  with  respect  to  observation  angles  delta  and 
alpha.  An  even  greater  expanded  set  of  summed  cross 
products  is  then  computed  and  combined  with  glint 
standard  deviation  to  yield  glint  lobe  width. 

The  RCS  lobe  and  glint  lobe  computations  are  re¬ 
peated  for  each  required  observation  angle  defined 
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The  real  power  of  the  RCS  program  lies  in  its 
ability  to  display  cross-section  output  results  in  a 
complete  and  elegant  manner  in  both  printed  and 
graphical  form.  This  section  discusses  the  standard 
display  options  available  in  the  program  followed  by 
additional  debugging  output.  The  standard  display 
options  are  RCS,  RCS  lobes,  glint,  glint  lobes,  and 
higher-order  statistics. 

The  commands  that  initiate  the  display  function 
are  GRAPH  and  PRINT  (defined  in  Section  5.5).  As 
the  names  imply,  GRAPH  produces  graphical  output 
on  a  terminal  or  a  plotter,  and  PRINT  produces 
printed  output  on  a  terminal  or  a  line  printer.  Since 
the  results  displayed  by  these  commands  are  general¬ 
ly  the  same,  only  differences  in  the  data  presented 
will  be  described  in  conjunction  with  the  overall  data 
available  for  display.  Program  results  are  not  altered 
because  of  their  being  displayed.  This  allows  multiple 
display  commands  to  be  included  in  the  data  base  in 
order  to  view  either  accumulated  results  or  new 
results  if  program  accumulators  are  cleared  (see  the 
CLEAR  directive). 

The  basic  program  shows  monostatic  full-scale 
RCS  as  a  function  of  observation  angle.  Both  aver¬ 
age  cross  section  (shown  with  the  optional  plus  or 
minus  standard  deviation)  and  relative-phase  cross 
section  may  be  displayed  in  square  meters  on  a 
logarithmic  scale  (i.e.,  in  decibels  relative  to  a  square 
meter).  Relative-phase  cross  section  may  also  be 
shown  in  square  meters  on  a  linear  scale.  Graphical 
RCS  output  also  lists  certain  program  parameters, 
such  as  frequency,  polarization,  and  percent 
shadowed,  on  the  edge  of  the  graph.  The  output 
from  this  section  is  always  displayed  as  part  of  the 
GRAPH  or  PRINT  operation. 

The  RCS  lobe  output  is  displayed  on  top  of  the 
basic  cross-section  output  in  graph  mode  and  follow¬ 
ing  it  in  print  mode  if  LOBE  is  active.  The  lobe  out¬ 
put  shows  the  number  of  RCS  lobes  per  hertz  of  fre¬ 
quency  and  the  number  of  lobes  per  spatial  increment 
around  a  circle.  In  graph  mode,  the  spatial  lobes  are 
displayed  inversely  to  yield  lobe  width.  Samples  of 
lobe  width  are  computed  around  a  circle  and  plotted 
as  an  ellipse  centered  at  I  he  observation  angle  of  in¬ 
terest.  The  major  and  minor  axes  of  the  ellipse  along 
with  the  orientation  angle  may  then  be  found  and  dis- 
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played.  Since  these  values  are  searched  out  rather 
than  computed,  lobe  width  must  be  sampled  often 
enough  (from  4  to  180  times  as  defined  by 
DELTA(4))  to  converge  on  the  required  accuracy. 

If  GLINT  is  active,  monostatic  glint  is  displayed  as 
a  function  of  all  the  independent  parameters  men¬ 
tioned  thus  far  and  also  the  angle  alpha,  which  de¬ 
fines  an  observation  plane  that  cuts  through  the  lobe 
structure.  Since  a  complete  graph  or  print  list  show¬ 
ing  glint  results  is  produced  for  each  value  of  alpha, 
only  certain  values  ,  namely  0  and  90°,  are  used.  The 
varying  independent  parameter  on  a  glint  graph  or 
print  list  is  the  observation  angle  theta.  The  depend¬ 
ent  parameter  represents  the  glint  moment  arm  dis¬ 
placement  in  degrees;  values  of  standard  deviation 
are  included  as  an  option. 

Monostatic  glint  lobe  output  is  displayed  in  a  man¬ 
ner  completely  analogous  to  cross-section  lobe  out¬ 
put.  It,  therefore,  similarly  shows  the  number  of 
lobes  of  glint  per  hertz  of  frequency  and  the  number 
of  lobes  per  spatial  angle  at  samples  evenly  spaced 
around  a  circle.  The  amount  of  glint  lobe  data  output 
is,  as  with  the  RCS  lobes,  completely  defined  by  the 
parameter  array  DELTA. 

The  final  standard  display  option  presents  higher- 
order  beta  statistics  computed  from  the  basic  RCS 
results.  The  type  and  amount  of  output  produced  by 
this  option  are  very  different  between  graph  and 
print  modes.  The  two  plots  generated  in  graph  mode 
show  skewness  and  kurtosis  as  a  function  of  beta 
width.  Distribution  acceptance  boundaries  are  drawn 
on  each  plot,  followed  by  the  actual  computed  data 
values,  which  are  annotated  with  their  corresponding 
observation  angle.  Many  additional  statistics  are 
computed  during  print  mode.  The  results  are  dis¬ 
played  within  the  normal  RCS  output  listing  rather 
than  at  the  end. 

Three  mutually  exclusive  debugging  display 
options  are  available  in  the  RCS  program.  All  three 
are  intended  to  provide  partial  results  corresponding 
to  each  separate  component  rather  than  the  com¬ 
bined  object.  When  the  cross-section  contribution  of 
a  component  is  computed,  it  is  immediately  displayed 
on  either  a  printer  (option  /P)  or  a  graphics  output 
device  (option  /G).  After  it  is  displayed,  the 
accumulator  vectors  are  reset  to  zero,  and  the  next 
component  is  computed  and  displayed  on  a  new 
listing  or  graph.  This  process  continues  until  the 
entire  data  base  has  been  read.  Prints  or  graphs  for 
which  there  is  no  component  contribution  are  not 
displayed.  Furthermore,  component  display  may  be 
aborted  in  graph  mode  by  entering  a  control  -C 
from  the  terminal  keyboard  for  each  component  not 
desired.  The  final  debugging  option,  /GC,  instructs 


the  program  to  combine  all  component  curves  onto 
one  graph. 

Debugging  display  options  may  not  be  used  in  con¬ 
junction  with  lobe,  glint,  or  higher-order  statistics 
output.  Any  PRINT  or  GRAPH  commands  found  in 
the  data  base  are  simply  ignored  by  the  program  in 
this  mode. 


S.4  Data  Base  Description 

All  inputs  to  the  RCS  program  must  be  suplied 
from  within  a  data  base.  A  data  base  is  a  file  that 
may  be  created  or  modified  by  an  editor.  It  contains 
all  the  instructions  necessary  to  control  the  program. 
Each  line  in  the  data  base,  except  comment  and  title 
lines,  constitutes  a  single,  complete  instruction  to  the 
program.  Each  instruction  line  begins  with  either 
four  or  five  characters,  defining  the  directive  type, 
followed  by  a  space  and  up  to  nine  numeric  parame¬ 
ters.  All  directives  (i.e. ,  instructions)  recognized  by 
the  program  and  their  functions  are  defined  in  detail 
in  Section  5.5. 

In  this  section  we  will  discuss  directive  interactions 
and  usage  protocols,  the  concepts  of  global  and  local 
directives  followed  by  their  relative  ordering  within 
the  data  base,  and  then  looping  capabilities  within 
the  data  base. 

All  valid  directives  in  the  RCS  data  base  fall  into 
one  of  two  functional  categories:  component  defini¬ 
tion  or  program  control.  All  component  definition 
directives  are  considered  to  be  local  parameters.  The 
program  control  directives  may  be  either  global  or 
local.  A  global  parameter,  once  defined,  retains  its 
specified  value  until  explicitly  altered  later  in  the  data 
base.  A  local  parameter  only  retains  its  assigned 
value  until  the  occurrence  of  a  component  definition 
directive,  which  has  the  effect  of  resetting  all  local 
RCS  program  control  directives  and  then  canceling 
itself.  The  local  parameters  must  therefore  be  reas¬ 
signed  before  the  occurrence  of  the  next  component 
directive,  if  required.  The  local  RCS  program  control 
directives  and  their  defaults  or  reset  values  are: 
ORGN  (zero),  RAMDB  (zero  dB),  TRANS  (identity 
matrix),  and  VALID  (all  angles). 

As  previously  indicated,  local  parameters  must  be 
ordered  so  that  all  control  directives  related  to  a 
component  directive  appear  before  the  component  in 
the  .data  base.  Global  RCS  directives  that  define 
radar  parameters  must  also  appear  before  the 
component  directives.  The  radar  directives  are: 
DELTA,  GLINT,  LOBE,  PHIR,  POLAR,  PYRT, 
RPRCS,  STATS,  THETA,  and  WLEN.  The  display 
directives  GRAPH  and  PRINT  (either  or  both  may 
be  specified)  must  appear  in  the  data  base  after  the 
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required  component  directives.  A  single  title  line 
must  always  be  placed  immediately  after  the  display 
directives.  The  two  directives  NEWTH  and 
NEWWL,  which  control  data  base  looping,  must 
always  be  placed  at  the  end  of  the  data  base,  if  used. 
If  both  directives  are  used,  NEWTH  must  precede 
NEWWL.  The  remaining  directives  may  appear 
anywhere  in  the  data  base  within  the  constraints 
given. 

The  RCS  program  may  be  instructed  to  reexecute 
the  entire  data  base  by  using  the  directives  NEWTH 
for  new  THETA  and/or  NEWWL  for  new 
wavelength  (WLEN).  The  maximum  number  of 
THETA  increments  allowed  on  one  pass  of  the  data 
base  is  181.  That  limit  may  be  extended  by  using  the 
NEWTH  directive  in  conjunction  with  the  parameter 
THETA(4).  The  actual  number  of  increments  used 
(NOI)  is  computed  as  1  +  (THETA(2)  - 

THETA(1))/THETA(3).  For  the  first  pass  through 
the  data  base,  the  observation  angle  ranges  from 
THETA(l)  to  THETA(2)  in  NOI  steps.  If  the 
directive  NEWTH  is  detected,  the  RCS  program 
stops  current  processing  and  starts  reading  the  data 
base  from  the  beginning.  If,  when  the  THETA 
directive  is  read,  THETA(2)  is  larger  than  or  equal  to 
THETA(4),  the  RCS  program  immediately  jumps  to 
the  first  line  following  the  NEWTH  directive.  The 
only  directives  that  may  follow  NEWTH  and 
NEWWL  are  comments.  If,  on  the  other  hand, 
THETA(4)  is  larger  than  THETA(2),  the  starting 
value,  THETA(l),  is  replaced  with  THETA(2),  and 
the  ending  value,  THETA(2),  is  replaced  with  the 
minimum  of  THETA(2)  +  NOI  and  THETA(4). 
This  processing  loop  continues  as  described  while 
THETA(4)  remains  larger  than  THETA(2).  All 
grapical  output  resulting  from  this  loop  is  placed  on 
the  same  graph.  Debugging  options  and  the  directives 
LOBE,  GLINT,  and  STATS  are  not  valid  if  THETA 
looping  has  been  specified. 

The  set  new  wavelength  directive,  NEWWL, 
instructs  the  RCS  program  to  reexecute  the  data  base 
in  a  manner  similar  to  NEWTH.  However,  in  this 
case,  each  time  the  directive  WLEN  is  read,  the  next 
value  in  the  parameter  list  will  be  used  as  the  new 
wavelength.  For  example,  if  the  WLEN  directive  has 
three  associated  parameter  values  from  a  maximum 
of  nine,  the  wavelength  will  be  set  to  the  first  param¬ 
eter  value  on  pass  one,  the  second  value  on  pass  two, 
and  the  third  on  pass  three.  The  RCS  program  will 
not  loop  for  pass  four  because  a  fourth  parameter 
value  was  not  specified  on  the  WLEN  command  line. 
The  only  valid  lines  following  the  NEWWL  directive 
are  comment  lines. 


5.5  Data  Base  Directives 

This  section  presents  a  complete  list  of  valid  RCS 
directives  and  associated  parameters.  The  directives 
have  been  divided  into  two  groups:  RCS  control  and 
component  definition.  A  complete  description  of  the 
available  RCS  components  is  given  in  Refs.  1 1  and 
12.  All  angles  must  be  specified  in  degrees  and  dis¬ 
tances  in  meters. 

5.5. 1  RCS  Control  Directives 

CLEAR  Clear  RCS -clears  all  ac¬ 

cumulated  radar  cross-sec¬ 
tion  results. 

COMM  n  Comment  -  ignores  the  next 

n  lines  in  the  data  base.  This 
instruction  may  be  used  to 
insert  a  specified  number  of 
comment  lines  in  the  data 
base. 

COME  Comment  End  -  ends  com¬ 

ment  block  mode.  This  in¬ 
struction  will  cancel  the  ef¬ 
fect  of  COMS.  The  next  data 
base  line  will  be  processed 
normally. 

COMS  Comment  Start  -  starts  com¬ 

ment  block  mode.  All  data 
lines  following  this  one  in  the 
data  base  will  be  ignored 
until  the  COME  directive  is 
found.  These  two  directives 
(COMS  and  COME)  provide 
the  user  with  an  efficient  tool 
for  ignoring  blocks  of  data. 
DELTA  n I,  .  .  .,«4  Delta  Angles  -  defines  the 
lobe  display  angles.  RCS 
lobe  and  glint  lobe  results 
will  be  displayed  at  angles 
ranging  from  n I  to  w2  in  in¬ 
crements  of  n 3  degrees.  The 
number  of  lobe  samples,  be¬ 
tween  4  and  180,  is  defined 
by  n4.  Only  delta  angles  that 
lie  within  the  theta  range  will 
be  used. 

GLINT  n  Glint  Switch  -  activates  the 

glint  lobe  calculations  if  n  is 
set  to  1  or  2  and  LOBE  is  ac¬ 
tive.  Glint  data  will  not  be 
computed  otherwise.  If  n  is 
set  to  2,  the  standard  devia¬ 
tion  will  be  added  to  the  dis¬ 
placement  curve. 
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GRAPH  n  1 . n9  Graph  Display  -  collected  re¬ 

sults  will  be  output  in  graph¬ 
ical  form  (see  Section  5.3). 

For  average  cross-section 

output,  plus/minus  standard  NEWWL 

deviation  curves  will  be 

added  to  the  graph  if  the 

absolute  value  of  nl  is  equal 

to  3  instead  of  I.  If  n  1  is 

negative,  terminal  graphs 

will  be  automatically  sent  to 

a  hard-copy  device.  The  x- 

axis  end  points  are  specified 

by  n2  and  «3.  The  number  of 

x-axis  labels  and  tick  marks 

per  label  are  specified  by  n4  ORGN  x,y,z 

and  nS.  The  y-axis  end  points 

are  specified  by  n6  and  nl. 

The  number  of  y-axis  labels 
and  tick  marks  per  label  are 
specified  by  n 8  and  n9.  The 
next  line  following  this  in¬ 
struction  is  always  used  for 
the  title.  This  instruction  will 
not  affect  any  RCS  results.  PAUSE 

GSIZE  nl,  ...  ,nS  Graph  Size -locates  the 
graph  boundaries  on  the 
plotter  page.  The  lower  left- 
hand  corner  of  the  graph  will 
be  placed  at  nl,n2  inches 
from  the  paper  corner  and  PH1R  n 
the  upper  right-hand  corner 
at  «3,n4  inches.  The  column  POLAR  n 
of  data  displayed  at  the  right 
of  the  graph  may  be  further 
displaced  by  n5  inches.  The 
default  values  for  nl  through  PRINT  n 
n5  are  2.25,  2.25,  8.75,  6.25, 
and  0  inches,  respectively.  If 
used,  GSIZE  must  appear 
before  the  GRAPH  di¬ 
rective. 

LOBE  n  Lobe  Switch  -  activates  the 

cross-section  lobe  calcula¬ 
tions  if  n  is  set  to  1.  Lobe 
data  will  not  be  computed 
otherwise. 

NEWTH  New  Theta -reruns  data 

base  with  new  theta  angle 
range.  This  instruction  will  PRYT  p,r,y,t 
cause  the  RCS  program  to  re- 
execute  the  data  base  from 
the  beginning  using  the  next 
theta  range  computed  by  the 
THETA  directive.  The  in¬ 


struction  will  be  ignored  if  no 
new  theta  range  is  required 
(see  THETA  and  Section 
5.4). 

New  Wavelength  -  reruns 
data  base  with  new  wave¬ 
length.  This  instruction  will 
cause  the  RCS  program  to  re- 
execute  the  data  base  from 
the  beginning  using  the  next 
wavelength  on  the  WLEN  di¬ 
rective.  The  instruction  will 
be  ignored  if  no  new  wave¬ 
length  is  available  (see 
WLEN  and  Section  5.4). 
Component  Origin -the 
component  following  this  in¬ 
struction  will  have  its  origin 
located  at  x,  y,  and  z  in  radar 
coordinates.  This  instruc¬ 
tion,  useful  for  lobe  and  glint 
computations  only,  will  be 
reset  to  zero  after  the  next 
component. 

Pause  Program -the  RCS 
program  will  stop  and  wait 
until  a  carriage  return  is  re¬ 
ceived  from  the  user  termi¬ 
nal.  This  command  should 
not  be  used  in  batch  mode. 
PHI  Radar -sets  radar  phi 
angle  to  n  degrees. 
Polarization  -  sets  system 
polarization  to  either  vertical 
(n  =  1)  or  horizontal  ( n  = 
2). 

Print  Display  -  collected  re¬ 
sults  will  be  output  in  printed 
form  (see  Section  5.3). 
Printed  output  will  be  dis¬ 
played  on  the  user  terminal  if 
n  =  0.  It  will  be  sent  to  a 
nonspooled  line  printer  if  n 
=  1 ,  a  spooled  line  printer  if 
n  =  2,  and  a  file  if  n  =  3. 
The  next  line  following  this 
instruction  is  always  used  for 
the  title.  This  instruction  will 
not  affect  any  RCS  results. 
Pitch,  Roll,  Yaw,  Tilt -de¬ 
fines  radar  rotation  angles. 
These  angles  will  be  used  to 
construct  the  radar-to-object 
coordinate  system  transfor¬ 
mation  matrix.  PRYT  will 
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default  to  the  identity  matrix 
if  not  defined.  If  all  angles 
are  zero,  the  resultant  matrix 
will  have  its  x  and  y  com¬ 
ponents  switched.1" 

RAMDB  n  Radar  Absorption  Materi¬ 

al -modifies  the  effective 
cross  section  of  a  compo¬ 
nent.  The  RCS  of  the  next 
component  following  this  in¬ 
struction  will  have  n  dB 
added  to  it  before  it  is  used. 
Thereafter,  the  value  of  n 
will  be  reset  toOdB. 

RPRCS  n\,  n2  Relative  Phase  Switch -ac¬ 

tivates  the  relative  phase 
RCS  calculations  if  nl  is  set 
to  1.  Relative-phase  data  will 
not  be  computed  otherwise. 
Relative-phase  cross  section 
will  be  displayed  on  a  loga- 
rthmic  scale  unless  n2  is  spec¬ 
ified,  which  defines  the  maxi¬ 
mum  linear  value  to  be 
graphed  in  square  meters. 

STATS  n  Statistics  Switch -activates 

the  beta  higher-order  statis¬ 
tics  calculations  if  n  is  set  to 
1.  Higher-order  statistics  will 
not  be  computed  otherwise. 

THETA  n\,  .  .  /|4  Theta  Angle  Range -sets 

radar  theta  range  and  incre¬ 
ment.  The  radar  observation 
angle,  theta,  will  vary  from 
nl  to  n2  degrees  in  increments 
of  n3.  The  total  number  of 
discrete  values  may  not  ex¬ 
ceed  181.  If  the  total  maxi¬ 
mum  angle,  n4,  is  present, 
theta  will  vary  from  nl  to  n4 
degrees  in  groups  of  (n2  - 
nl)/n3  increments.  Since  this 
option  reruns  the  RCS  pro¬ 
gram  for  each  THETA 
group,  the  number  of  dis¬ 
crete  observation  angles  al¬ 
lowed  becomes  unlimited. 
The  command  NEWTH 
must  be  included  at  the  end 
of  the  input  stream  (but  be¬ 
fore  the  directive  NEWWL, 
if  present)  to  force  the 
program  to  loop  to  the  next 
THETA  group  (see  Section 
5.4).  If  the  parameter  n4  is 


specified,  the  following  com¬ 
mands  are  not  allowed: 
LOBE,  GLINT,  STATS,  and 
debugging  options. 

TRANS  nl,  .  .  .  ,n9  Transformation  Matrix  as¬ 
signs  component  orientation. 
The  nine  values  nl  through 
n9,  specified  in  degrees,  de¬ 
fine,  by  row,  the  direction 
cosine  matrix  of  the  compo¬ 
nent  following  this  instruc¬ 
tion.  The  transformation 
matrix  will  be  checked  for 
validity  before  being  used  by 
a  component  and  then  reset 
to  the  identity  matrix  after 
being  used. 

VALID  nl,  .  .  .  ,n8  Valid  Range -sets  range  of 
valid  angles  for  component. 
The  values  of  nl  through  the 
maximum  n 8,  specified  in  de¬ 
grees,  will  define  observation 
angles  for  which  the  follow¬ 
ing  component  is  visible.  The 
parameter  values  must  be 
specified  in  pairs  and  imply 
the  following:  nl,  n2  =  first 
valid  range  of  THETA;  n3, 
n4  -  associated  valid  range 
of  PHIR  (all  PHIR  if  not 
present);  n5,  n6  =  optional 
second  valid  range  of 
THETA;  n7,  n8  =  associ¬ 
ated  valid  range  of  PHIR. 
All  ranges  are  assumed  to  be 
counterclockwise  from  the 
first  element  in  a  pair  to  the 
second.  The  valid  range  will 
be  reset  to  all  angles  after  be¬ 
ing  used  for  a  component. 

WLEN  nl,  n2  .  .  .  Wavelength  -  assigns  pro¬ 
gram  wavelengths.  This  in¬ 
struction  sets  the  initial  pro¬ 
gram  wavelength  to  nl 
meters.  Each  following  pass 
through  the  data  base  in 
looping  mode  will  set  the 
program  wavelength  to  the 
next  n  on  the  command  line 
(pass  2  would  use  a  wave¬ 
length  of  n2,  etc.)  up  to  a 
maximum  of  n9.  The  instruc¬ 
tion  NEWWL  must  be  used 
at  the  end  of  the  data  base  if 
multiple  wavelength  passes 
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(i.e.,  looping)  are  desired  (see 

center  to  nearest  side  length  1®, 

Section  5.4). 

a.  If  ogive  is  truncated,  pa-  • 

rameter  b  and  angle  A 

$ 

5.S.2  Component  Definition  Directives 

specify  the  amount  of  trun- 

CAVA  a 

Cavity  -  defines  a  cavity  with 

PARAB  p 

cation  (b  is  zero  otherwise). 

Parabola  -  defines  a  parabo- 

m 

equivalent  circular  area  a. 

la  where  the  parabola  equa- 

CAVD  d 

Cavity  -  defines  a  cavity  with 

lion  is  p  =  (jr**2  + 

circular  diameter  d. 

y"2)/-4z.  .>■ 

CFLAT  r 

Circular  Flat  Plate  -  defines 

RCORN  a,b,c 

Rectangular  Corner  -  defines  \-j 

a  flat  plate  with  radius  r. 

a  rectangular  corner  with  V-\ 

%N 

CONE  A,n,a, 

Truncated  Elliptical  Cone- 

edges  a,  b,  and  c. 

LI,  L2 

defines  a  cone  with  lengths 

RFLAT  a,b 

Rectangular  Flat  Plate -de- 

LI  and  L2  from  the  im- 

fines  a  flat  plate  with  sides  a  Wf 

g$ 

aginary  cone  tip.  Cone  tip  to 

SOLID  a,b,c 

and  b. 

Li  is  truncated.  The  elliptical 

Solid  -  defines  an  ellipsoid  ^ 

cone  base’s  major  radius 

with  radii  a,  b,  and  c.  S\" 

length  is  a  and  the  major-to- 

TCORN  a.b.c 

Triangular  Corner  -  defines 

hi 

minor  ratio  is  n.  The  cone 

a  triangular  corner  with 

• 

half-angle  is  defined  by  A . 

edges  a,  b,  and  c. 

CYLIN  a.b.L 

Elliptical  Cylinder -defines 

TORUS  a,b 

Torus- defines  a  torus  with  v,-'. 

» 

a  cylinder  with  major  radius, 

inner  radius  a  -  b  and  outer 

a,  minor  radius  b,  and  length 

radius  a  +  b. 

L. 

TWAVE  p.L.g 

Traveling  Wave -defines  a 

D1HED  a,b,c 

Dihedral  -  defines  a  90°  di- 

traveling  wave  with  the  pa-  \V" 

|2 

hedral  with  sides  a  and  c  and 

rameters  p  (relat  i ve-phase  ve- 

■ 

height  b. 

locily),  L  (length),  and  g  (re-  5f1 

ECORN  a,b,c 

Elliptical  Corner -defines  an 

flection  coefficient). 

8 

fcd 

elliptical  corner  with  edges  a. 

WEDGE  A,L,B 

Wedge -defines  a  wedge 

b,  and  c. 

with  angle  A  and  height  L.  c'S 

LWIRE  a 

Loop  Wire -defines  a  cir- 

Parameter  B  is  used  to  limit 

cular  wire  loop  with  radius 

the  cross  section  to 

■ 

a. 

S/WLEN**2. 

5 

OGIVE  R,a,b,A 

Ogive  -  defines  an  ogive  with 

WIRE  L,a 

Wire -defines  a  wire  with 

o 

radius  of  curvature  R  and 

length  L  and  wire  radius  a. 

A 

C.v 

tv 

4-*  - 

K* 

.s'/ 

S'.- 

V.* 

If 

1 

i 

ac 

,vr 

S' 

VJ! 

m 

3 
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6.0  PROGRAM  LISTING 


»>  FILE  ■  8Y0:t111,12)RCS.FTN 


27-SEP-83  14:92:52  <<< 


RRRR  CCCC  SSSS 
ARC  S 

RRRR  C  SSS 

R  R  C  S 

R  R  CCCC  SSSS 


RADAR  CROSS  SECTION 


NOTE!  THIS  PROGRAM  STORES  MATRICES  BY  ROWS  (NOT  COLUMNS) 

MAXRCS  CONTAINS  THE  CURRENT  MAXIMUM  RCS  BUFFER  SIZE  1181) 

REALMS  T1 ,T2,T3,TA,T5,T8,SO.ZERO 

REALMS  CMOS (26) .OBJSI18) , RCS MSI  I ,RCS2(181 ) .TRANSIS) .VALIO(S) 
REALMS  ARRIS) .PRYTI9) , MAT  19) .C0MPI9) .PRYTV(S) , PRYT2I9) 

LOQICALM2  ATTN 

LOQICALMl  F I LE (SO) ,BUF(S0) .SLASH. AQ.AP, OBJECT, AZ. AC, AN, 00T.DRCSI4) 
REALMS  POLVHI2) , YYVI181 ,2) ,XV(1S1 ) , SL0BE(S, IS1 ,2) 

EQUIVALENCE  (FILE.BUF) , (YYV.SLOBE) , (XV.SLOBEI1 ,1 .2) ) 

REALMS  I  DENT  1 9 ) . ORARR 1 9 ) , S I DE  U , 4 ) , DEBOUT ( 0 1 , ORQN ( 3 1 . M60DB 
REALMS  DELTA(S) .ALPHA ( 4 ) , SCALE ( 6 ) ,THETA(4I .GSIZE(S) 

INTEQERM2  LSIDEI4) .ROUTE 

INTEQERM2  DEBUG , COMB  IN, 9L I  NT 

LOB  I CALM2  GOULD , LP , OPENE , LPF . OLDGRF .MQRAPH . QDAUTO 
COMMON/SWITCH/  I  Of , THETA , GOULD , LP , OPENE , LPF , OLDGRF .MQRAPH . 

1  GDAUTO.DE BUG, COMB  IN, GLINT, IBUQ.LBUG, I  STATS. IZETA .STATS.NRCS, I OR 

COMMON/RCSCOM/CUTOFF .WLEN.PI , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP . SPHP , 
S  C2PHP.S2PHP, RAD, POLAR, LOBE, 

S  ARR 

COMMON /RCS LOB/ I OL , ORQN, COMP , SLOBE .BETA . 

S  ANQ . CANS , SANQ .PHR.CPHR, SPHR , KLOBE 
COMMON  /RPRCS/  I  RPR, I RCS, RAM, RCS, RCS2 
COMMON  /GRCOM/  QRLIB 

DATA  POLVH/’ VER HOR ’/, CLEAR/ ’ CLEA ’/ , COME/ ’ COME ’ / , SLASH/ ’ / • / 

DATA  CMOS/ 'COMM' ,'COMS' ,'WLEN' , ’ POLA' , ' TRAN’ ,'VALI' ,’PHIR’ ,’THET' , 
S  'QRAP' , ’PRIN’ ,’PRYT’ , 'PAUS' , ’CLEA’ . 'DEBU' ,  NEWW’ , ’ RAMD ’ , 'LOBE' . 

S  ’BETA’ , 'ORQN' ,’DELT’ ,’ALPH’ ,’RPRC ,’STAT’ ,’QLIN' ,’NEWT’ ,’QSIZ’/ 

DATA  OBJS/’OGIV’ ,’CYLI' ,’CONE' , • PARA' , ’ SOL  I ’ , ’ TORU’ , ’ WE DO’ , 

S  ’WIRE’ ,’LWIR’ , ’CAVA’ ,’CAVD’ ,’CFLA’ ,’RFLA’ ,’TWAV’ ,'DIHE’ ,’TCOR’ , 

S  ’RCOR’ ,’ECOR’/ 

DATA  IDENT/t . ,0. ,0.,0.,1.,0.,0. ,0. , 1 ./.AQ/’Q’/.AP/’P’/.AC/’C’/ 
DATA  DEBOUT/O.,’  COM' , ’ PONE ’ . ’ NT  0’ , ’ UTPU’ , ’ T • / , ALL/’ ALL  •/ 

DATA  QRARR/-1 ,0,180, 7, 3, >90, 60, 7, 1/ , AZ/’ Z’ / , AN/’ N' / .DOT/’ . ’ / 

OATA  SIDE/’ STAR' ,’ BOAR’ ,’D’ , ’PORT’ ,0,0, ’ TOP’ ,0,0, ’BOTT’ , ’OM’ ,0/ 
DATA  LSIDE/7, 4, 3,8/, DRCS/' .’ ,’R’ ,’C’ ,’$’/, ZERO/O.DO/ 

DATA  Q8 I ZE/2. 29, 2. 25, 8. 75, 8.29,0./ 

DEBUG  OUTPUT 

IBUG»1«>RCS  OUTPUT 
I BUa«2«>TRANS , PRYT , COMP 
| BUa>4<> ANGLES ( RADARAPR IME ) 
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c 

C  SET  RCS  CONSTANTS 
C 

NCMDS26 
NOBJsIS 
MAXRCS=1 81 
P|s3. 1415926 
PI2SPI/2. 

RADsPI/180. 

CUTOFF *1 .E-10 
LINEsO 

QOULDs. FALSE. 

C 

C 

C  ASSIGN  LUNS  AND  ATTACH  INPUT  FILE 
C 

I  OR  si 
IOtfs2 
IOPS3 
IOQS4 
IOFS5 
lOLsfl 
IOE*7 
C 

CALL  ASNLUNI IOR, ' Tl * ,0) 

CALL  ASNLUNI  lOW.’TT  ,0) 

CALL  ASNLUNI IOF, ’ST* ,0) 

CALL  ASNLUNI I  OP,’ QD’ ,0) 

CALL  ASNLUNI lOL.'SY* .01 
CALL  ASNLUNI lOE.'SY' ,0) 

OPENEs. FALSE. 

C 

WRITE ( IOW, 1 ) 

1  FORMAT  I*  ENTER  RCS  DATA  FILE  NAME’) 

READIIOR.2.ENDSSS99)  FILE 

2  FORMAT  I S0A1) 

CALL  TRIMSOINF ,F ILE) 

F I LE (NF*1 ) »0 
C 

C  INPUT  OPTIONS 

C 

C  /Z  s  DRAW  GRAPHS  ON  ZETA  PLOTTER 

C  /P  s  PRINT  RESULTS  FOR  EACH  COMPONENT 

C  /G  s  GRAPH  RESULTS  FOR  EACH  COMPONENT 

C  /QC  s  DRAW  ALL  /Q  OUTPUT  ON  ONE  GRAPH  (COMBINE) 
C 

DEBUG *0 
COMBINsO 
I ZETA SO 
IF*1 
I  DOT SO 

DO  6  J*IF,NF 
JJsJ 

IFIFILEIJJ) .EQ.DOT)  IDOTsJJ 
IFIFILEIJJ) . EQ. SLASH)  GO  TO  7 

6  CONTINUE 
GO  TO  9 

7  NFsjJ-1 

IFIFILEIJJ+1) .EG.AZ)  IZETAs200 
IFIFILEIJJ*!) . EQ.AGI  DEBUGsl 
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1 200»  IF IFILEIJJ+1 ) .EQ.AP)  DEBUQ-2 

1  2 1  0 >  IF  I  DEBUG . EQ. 1 . AND. F I LE I J J  *2 1 . EQ. AC)  COMB  I  Ns -1 

1 220>  I F  =  J  J ♦ 2 

1 230>  •  IF  1 1  DOT. BT. 01  QO  TO  9 

1 240>  C  DEFAULT  FILE  EXTENT  ION  TO  .RCS 

1 250>  DO  13  J:l,4 

12B0>  NFsNF+1 

1 270*  FILE(NF)sORCSIJ) 

1 280>  13  CONTINUE 

1290>  9  F ILE INF+1 I =0 

1 300>  OPEN  I UN  I Ts I OF , NAME  a F I LE , TYPE s ’ OLD’ .READONLY) 

131 0>  C 

1 320>  TIME=SECHDS(0.0) 

1 330>  CALL  QRINITIO, IZETA, IOQ, IER) 

1 340>  CALL  GRHV('H’) 

1 3 5 0 >  CALL  QRP 1C  12 .25 ,2.25,8.73,6.25) 

130O>  IFI IZETA. EQ.O)  CALL  QRDEVI ’TEN* . ’HARD’ , IER) 

1 370>  C 

1 3B0>  C  INITIALIZE  RCS  VECTOR  A  OTHER  STUFF 
1 390>  C 

1400  CMD*0. 

1410  IWL*1 

1420  OLDQRFs.  FALSE. 

1430  QO  TO  70 

1440>  C 
1450  C 

1460  C  IHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHI 

1470  C  READ  NEXT  LINE  IN  RCS  DATA  FILE 

1 460>  C  IHHHHHHHMHMHMHHHHHMHHHHMMHMHHHMI 

1490  C 

1 500>  3  CONTINUE 

1510  IFIOBJECT. AND. DEBUG. E0.1)  GO  TO  45 

1S20>  IFIOBJECT. ANO. DEBUG. EQ. 2)  GO  TO  50 

1 530*  C 

1 540>  LINE*LINE*1 

1350  READ! IOF ,4,END=9990,ERR=9993)  CMD.ARR 

1 560>  4  FORMAT (A4 , IX ,9F 10 . 0) 

1 S70>  C 

1 5S0>  C  DETERMINE  IF  CMD  IS  A  DIRECTIVE 
1590  C 

1600>  DO  5  Jsl.NCMD 

1610  IF (CMD.EQ.CMDSI  J) ) 

1020  6  GO  TO  (10,12,15,20,25,30,39,40,45,50,55,60,70,60.66,65. 

1630>  9  500,510,520,540,550,560,670,560,590,600),  J 

1 040>  5  CONTINUE 

1 650>  C 

106O>  C  QO  TO  OBJECT  SEARCH  IF  NO  DIRECTIVE  MATCH  FOUND 
107O>  C 

10SO>  QO  TO  90 

1690>  C 
1700>  C 
1 71 0>  C 
1720>  C 

1 730>  C  EXECUTE  APPROPRIATE  DIRECTIVE 

1740>  C  «MN«N»««NN«KN««ltlt*NN«ll««NII««N»« 

1750>  C 
1760>  C 

1770>  C  "COMM”  -  SKIP  COMMENT  LINES 

1760>  C 

1 790>  10  NCOMMiARRM) 
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1S00> 

IF INCOMI. LE . 0)  Q0  TO  3 

1S10> 

DO  11  J  * 1 , NCOMM 

1I20> 

LINEsLINE+1 

r*' 

1930>  11 

READ! I0F |4, END-9990)  CMD 

1S40> 

Q0  TO  3 

1950>  c 

"COMS"  •  SEARCH  FOR  "COME1* 

1960>  12 

READ! IOF,4,END:9993)  CMD 

1970> 

LINEbLINE+1 

■ 

1880> 

1 F (CMD. EQ. COME  1  QO  TO  3 

1390> 
1900  C 

QO  TO  12 

1910  C 
1920>  C 

"WIEN"  -  SET  WAVELENGTH 

ft 

1930  15 

WLEN=ARR( IWL) 

tv 

1940> 

IFIWLEN. LE .0. AND. IWL . EO. 1 )  GO  TO  9991 

19S0> 

IF (WLEN.LE .0)  GO  TO  9990 

jn 

1 990> 

1 970>  C 

QO  TO  3 

'vt 

1 990>  C 

"POLAR"  -  SET  POLARIZATION  VECTOR 

1 990>  C 
2000>  20 

POLAR=AMAX1 11 . , AM  INI (2 . , AAA  11 ) 1 1 

2010> 

GO  TO  3 

2020>  C 

« 

2030>  C 
2040>  C 

"TRANS"  -  COMPUTE  COORDINATE  TRANSFORMATION  MATRIX  1 

2050>  25 

CONTINUE 

$ 

20S0> 

DO  29  Jal ,9 

f  >i 

2070>  29 

TRANS (J ) =COS 1 ARR 1 J) KRAD) 

Uu 

2050> 

CALL  TRANSP (TRANS .MAT) 

2090> 

CALL  MATMULI TRANS, MAT. MAT) 

rn 

2100>  C 

MAKE  SURE  TRANS  IS  VALID 

■ 

2 1 1 0> 

DO  24  J«1,9 

■ 

2120> 

IF (ABS (MAT  1 J ) - 1  DENT! J) ) .QT. 1 . E-3)  QO  TO  27 

21 30>  24 

CONTINUE 

2140> 

QO  TO  3 

ft 

21 50>  27 

WRITEI2.29) 

sl'l 

2160>  21 

FORMAT  1 ’  RCS>  —  BAD  TRANS  DIRECTION  COSINE  MATRIX’/’  RCS>  ) 

21 70» 

WRITE  1 2, 29)  ARR 

21 50>  29 

FORMAT!’  RCS>  \3F10.3) 

■ 

2190> 

WRITE (2. 29) 

2200> 

WRITE  1 2, 29)  MAT 

w 

2210> 
2220>  C 

QO  TO  9993 

2230>  C 
2240>  C 

"VALID"  -  SET  VALID  ANGLE  REGIONS 

B?fl 

22S0>  30 

CONTINUE 

2260> 

DO  31  Jal ,8 

2270>  31 

VALID! J)<ARR(J) 

m 

::”j 

2290> 
2290>  C 
2300>  C 
231 0>  C 
2320>  35 
2330> 
2340>  39 
2950> 
2990>  C 
2970>  C 

QO  TO  3 

"PHIR"  -  SET  RADAR  PHI  ANGLE 

PHRaARRM) 

IF(ABSIPHR) .GT. 190. )  WRITE! I0W, 39)  PHR 

FORMAT  1 ’  RCS>  —  PHI R*’ , F6 . 1 , ’  19  OUT  OF  BOUNDS’) 

GO  TO  3 

"THETA”  -  SET  THETA  RANGE  IMIN.MAX, INC)  A  TOTAL  MAX 

3» 

2990>  C 

3® 

2990>  40 

IFIOLDGRF)  QO  TO  43 

iC 
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2400> 

THETA  11  l-ARRM  1 

• 

241 0> 

IF IARRI3I . LE .0. )  GO  TO  42 

-•r* 

*  4 

2420> 

THETA  12 J*AMAX1 (THETA! 1 ) , ARRI2) ) 

2430> 

THETA (3) SARA (3) 

2440> 

THETA(4)=ARRI4) 

•  J 

2450> 

MQRAPHs. FALSE. 

2460> 

1 F ( THETA  141 .GT. THETA  121)  MGRAPHs . TRUE . 

“ 

247Q> 

C 

2490> 

NRCS= (ARR (2) -ARR  Mil  /ARR  (31*1  .001 

2490> 

1 F ( NRCS . LE .MAXRCS )  GO  TO  3 

2500> 

WRITE  1 10V ,41 ) 

2  5 1 0  > 

41 

FORMAT ( ’  RCS>  --  TOO  MANY  OUTPUT  POINTS’) 

2920> 

GO  TO  9993 

2930> 

C 

v 

2940> 

42 

THETA(2)sARR(1) 

2990> 

THETA (31=1  . 

■» 

2990> 

NRCS=t 

2970> 

GO  TO  3 

2980> 

C 

'/l 

2990> 

43 

TsTHETAMI 

2600> 

THETA (1 )sTHETAI2) 

291 0> 

THETA  12 1 sAMI N1 (2 .HTHETAU )-T,THETAI4) ) 

2920> 

GO  TO  3 

2930> 

C 

ta 

2940> 

C 

"GRAPH"  -  GRAPH  RESULTANT  RCS  VECTOR  s  CMDI9) 

w 

2990> 

c 

2990> 

49 

CONTINUE 

* 

2670> 

C 

IFIISTATS.EQ.il  GO  TO  4942 

*  ■ 

2690> 

IFINRCS.LE.1 .OR. DEBUG. EQ. 2)  GO  TO  3 

• 

2690> 

1 F ( 1 ZETA. EQ.O. AND. COMB IN.LE.O. AND. . NOT.OLDQRF) 

2700> 

1  CALL  GRMODEC  ERASE'  ,3) 

2710> 

ATTN*. FALSE. 

2720> 

1 F 1 CMD . EQ . CMOS  191)  GO  TO  4949 

Sf 

2730* 

c 

2740> 

C  GRAPH  DEBUG  MODE 

2790> 

M90DBS1./10MM9 

V 

2760> 

DO  4941  Jsl , NRCS 

v\ 

v  * 

2770> 

IFIRCSIJ) .GT.M90DBI  GO  TO  4940 

’.V 

2790> 

4941 

CONTINUE 

2790> 

GO  TO  4999 

2900> 

C 

3 

2910> 

4940 

CALL  TRAPCCI 1  OR, ATTN) 

> 

29  20> 

DO  4949  J*1 ,9 

2930> 

4949 

ARRI J)sQRARR(J) 

2940> 

4949 

CALL  QRSCLURRI2) , ARRIS) , ARRIS) ,ARR(7)I 

u- 

2990> 

IF  I IRPR.GT.O. AND.RPRCSG.GT.O. I CALL  GR8CLIARRI2) ,0. , ARR 13) .RPRCSG) 

2»90> 

SCALE  ID ‘ARR (2) 

2970> 

SCALE (2lsARR(6) 

2990* 

SCALE  13) >ARR(3) 

2»90> 

SCALE  1 4) SARR 1 7 ) 

2900> 

SCALE  19) SAAR (41 

\< 

2910> 

SCALE  19) >ARR 19) 

2920> 

IFIC0MBIN.EQ.1 .0R.0LD8RF)  GO  TO  4947 

2930* 

C 

V* 

2940> 

YsGSIZE(2) 

2990> 

IF IKLOBE. LE. 0)  GO  TO  44 

2990> 

Yst-,99 

2970> 

CALL  QRPLOTIG8IZEI3) .GSIZEI2) ,01 

29f0> 

CALL  QRPLOT (GSIZE 11 ) .GSIZE 12) , 1 ) 

■*  * 
s‘( 

2990> 

C 

m 
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3000 

3010 

3020 

3030 

3040 

3050 

3000 

3070 

3000 

3000 

3100 

3110 

3120 

3130 

3140 

3150 

3100 

3170 

3100 

3100 

3200 

3210 

3220 

3230 

3240 

3250 

3200 

3270 

3200 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3300 

3370 

3300 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

3400 

3470 

3475 

3490 

3490 

3500 

3310 

3520 

3530 

3540 

3550 

3500 

3570 

3500 


44  CONTINUE 

CALL  GAAXIS(0,GSIZEI1)  ,  Y ,QSI ZE (3) -Q8IZE ( 1 1  ,0.  ,  AHA  12 1 ,  AAR  (31 . 

1  IF  IX (AHA  14) ) , IP  I X (AAR  1 5 1 ) , 'OBSERVATION  ANGLE’, 99,1) 

DYsQSIZE (4) -GSIZE (2) 

I F II APR . EQ . 0 . OR . RPRCSG . EG . 0 . )  GO  TO  4500 
C 

CALL  GRAXmS.GSIZEM)  .QSIZE12)  ,DY,90.  ,0. . RPRCSG. LRPRC8, 2, 

1  ' RCS  ISM)’, 9,1) 

GO  TO  4547 
C 

4500  CALL  GRAXIS(3,Q$IZE(1 ) , GS I ZE 1 2 ) ,DY,90. .ARRIO) , AAA  17) , IF  I X I AAR  10 1 ) . 
9  IF IX (AAR (91 ) , 'RCS  (DBSM) ’ ,10,1) 

C 

4547  ANQsTHETA  1 1 ) 

00  40  J*1,NACS 
XV <J) SANG 
ANQ>ANG*THETAI3> 

YYYsRCSI J ) 

IF « IRPR.EQ.O)  GO  TO  4501 
YYY=RCS I J )N«2«RCS2 ( J ) mm2 
IF (RPRCSG. EQ.O. )  GO  TO  4501 
YYVIJ , 1 ) sYYY 
GO  TO  40 

4501  YYVIJ,  1)310. «ALOQ10(AMAX1  (CUTOFF, VYYI ) 

40  CONTINUE 

CALL  QRVI NO ( -NRCS , XV , Y YV , 18,0, .1 .0. .1) 

IF (COMBIN.EQ. 1 )  GO  TO  4553 

IFIABSIARRII ) I.LE.1. OR. DEBUG. EQ.1. OR.IRPR.EQ.il  GO  TO  49 
SIGN*-1. 

DO  47  K»1,2 
DO  49  Jsl.NRCS 

S>SQRT ( AMAX1 (0. , RCS( J )H*2-RCS2( J) ) ) 

49  YYVIJ, K)«10.NAL0G10(AMAX1 (CUTOFF, RCS(J)«SIGN«S) I 

CALL  QRWINDI-NRCS.XV.YYVd ,K) , IB.O, .1 .0. ,1 ) 

47  SIQNE1. 

49  IF( IZETA.EQ.O)  CALL  QRMODE ( ’ ALPHA’ ,2) 

IF (CMD .EQ.CMDS t 9) )  GO  TO  4540 
DEBOUT  1 1 ) *CHD 

IF (COMBIN.EQ. -1 )  DEBOUT II ) BALL 
CALL  GRTXT (DEBOUT, 21) 

GO  TO  4549 

4541  READIIOF.2.ENDS9993)  BUF 
LINE*LINE»1 
IF (OLDGRF )  GO  TO  4559 
CALL  TRIMOOINB.BUF) 

CALL  GRTXT ( BUF, NB) 

QDAUTO>. FALSE. 

IFIARRI1 1 .LT.O. )  QDAUTO« .TRUE . 

4549  XMID* IQSIZE (3) *6SIZE ( 1 ) )/2 . 

CALL  GRSX Y ( XM I D , QS I ZE 1 4 ) * , 75 , XX , YY ) 

CALL  QRPRNT(0,XX,-1 ,YY, .13,0.) 

CALL  GRTXT I ’MONOSTATIC  FULL-SCALE  ’,0) 

IFI IRPR.EQ.O)  CALL  GRTXT ( 'RCS' ,0) 

IFI IRPR.EQ. 1 )  CALL  GRTXTI ’ APACS’ ,01 
CALL  QRSXY IXMID.QSIZE 1 4) ♦ . 5 , XX , YY) 

CALL  GRP ANT (0, XX , -1 , YY , . 1 5 ,0 . ) 

C 

C  SET  TERMINAL  TO  ALPHA  MODE  AND  SCALE  TO  INCHES 
IFI IZETA.EQ.O)  CALL  QRMODE I’ ALPHA’ ,3) 

CALL  QASCL (QS I ZE 1 1 ) , Q 8 1 ZE ( 2 ) .QSIZEI3) .GSIZEI4) ) 
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3S90> 

3600>  C 
3010> 

392Q> 

3630> 

3640> 

3650> 

30S0> 

3970> 

360O> 

3090> 

3700» 

3710> 

3720> 

3730> 

3740* 

3750> 

3760> 

3770> 

3790> 

3790> 

3900>  4942 
391 0> 

3820> 

3930> 

3840> 

39S0> 

3990> 

3370> 

3880> 

3990> 

8900> 

3910> 

3920> 

3980> 

3940> 

8990> 

3990> 

3970> 

3990> 

3990> 

4000> 

401 0> 

4020> 

4030> 

4040> 

4090> 

4090> 

4070> 

4090> 

4090> 

4100>  4199 
41 10» 

4120> 

4130> 

4140*  4997 
4190> 

4190> 

4170>  C 
4190>  C 


X=QSIZE(3) e1 . 05+QSIZE 15) 

Y«QSIZE(4)*.7S 
FREQs . 3/WLEN 

CALL  QRTXT I ’FREQUENCY  ’  ,0) 

CALL  QRPRNT ( 1 ,X .0 , Y , . 12 , 0 . I 
CALL  QRNUMI ’ F  * , - 1 , 3 , FREQ) 

CALL  QRTXT ( ’  QHZ ’ ,0) 

CALL  QRPRNTI-1 ,X,0, Y. .12 ,0. ) 

Y*Y- .  2 

CALL  QRTXT (’POLARIZATION  ’.0) 

CALL  QRPRNTI1 ,X,0,Y, .12,0.1 
CALL  QRTXT (FOLVHI IF IX (FOLARI 1,3) 

CALL  QRPRNTI-1 ,X,0,Y, .12,0.1 
IF IMQRAPH)  QO  TO  4542 
YSY-.2 

PERSH*1 09 . -1 00  .llFLOATINVAL  ID) /AMAXO 1 1 , NTOTAL) 
CALL  QRTXT ( 'X  SHADOWED  ’ ,01 
CALL  QRPRNTM  ,X,0,Y,  .12,0. 1 
CALL  QRNUMI ’ F ’ , -1 , 1 .PER9H) 

CALL  QRPRNT I -1 ,X,0,Y, .12,0.) 

Y*Y- .  2 

CALL  QRTXT I 'PITCH  ’ ,0) 

CALL  QRPRNTM  ,X,0.Y,  .12,0. 1 
CALL  QRNUM I ’ F ’ , - 1 , 3 , PRYTV M 1 1 
CALL  QRPRNT I -1 ,X,0,Y, . 12,0. I 
Y=Y-.2 

CALL  QRTXT ( ’ ROLL  ' ,0) 

CALL  QRPRNTM  ,X,0, Y ,  .  12 ,0 . 1 
CALL  QRNUMI ’ F’ , -1 ,3,PRYTV(2I I 
CALL  QRPRNT I -1 ,X,0, Y , . 12 , 0. I 
YSY-.2 

CALL  QRTXT I ’YAW  ’,0) 

CALL  QRPRNTM, X.O.Y..  12, 0.1 
CALL  QRNUM I ’ F ’ ,*1,3, PRYTV (31) 

CALL  QRPRNTI-1 ,X,0.Y, .12,0.) 

YSY-.2 

CALL  QRTXT (’TILT  ' .0) 

CALL  QRPRNTM  ,X,0,Y,  .12,0. 1 
CALL  QRNUMI ’ F' , -1 , 3 , PRYTV (4) ) 

CALL  QRPRNTI-1 ,X,0,Y, .12,0. 1 

Y»Y-.2 

ISIDEiO 

IFIPHR.EQ.90.)  ISIDEH 
IFIPHR.EQ.-90.)  ISIDE-2 
IFIPHR.EQ.O.)  1 9 1 DE*3 
IFiPHR.EQ.190.)  I9IDE>4 
IFIISIDE.EQ.O)  QO  TO  4559 
CALL  QRTXT ( 9 1 DE 1 1 , 13IDE)  ,ISIDEMSIDE) ) 

CALL  QRPRNT(0,X,0,Y, .12,0.) 

QO  TO  4597 

CALL  QRTXTi’PHIR  ’,0) 

CALL  QRPRNTM, X,0,Y,. 12,0.) 

CALL  QRNUMI’ F' ,-1 ,3,PHR) 

CALL  QRPRNTI-1 ,X,0,Y, .12,0.) 

CONT I NUE 

IF  IMQRAPH)  QO  TO  4559 
IF IDE9UQ.NE.0)  QO  TO  4550 

COMPUTE  AND  QRAPH  LOBE  STATISTICS  IF  REQUIRED 
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4190: 

4200: 

4210: 

4220: 

4230: 

4240: 

4290: 

4200: 

4270: 

4200 

4200 

4300 

4310 

4320 

4330 

4340 

4390 

4300 

4370 

4330 

4390 

4400 

4410 

4420 

4430 

4440 

4490 

4400 

4470 

4490 

4490 

4900 

4910 

4920 

4930 

4940 

4990 

4900 

4970 

4900 

4990 

4000 

4010 

4020 

4030 

4040 

4090 

4000 

4070 

4000 

4090 

4700 

4710 

4720 

4730 

4740 

4790 

4700 

4770 

4790 


C 

IF (KL06E . LE .0)  GO  TO  4991 
CALL  GUP LOTI 11 . ,0. ,0) 
i qhprsi 

WRITE! I OL' II  IGRPR , K  LOBE , GL I  NT , NRCS , NLOBE , THETA .ALPHA . DELTA , 

1  BETA , SCALE , IZETA , I PCODE ,WLEN , L9UG , NB ,BUF , QDAUTO, 

2  FREQ , POLVH 1 1 F I X I POLAR 1 1 , PERSH . P I TCH , ROL  L . YAW . T I LT , PHR , GS I ZE 
CLOSE (UNITs | OL) 

CALL  SPAWNS 

OPEN  I  UN ITs I OL, NAME s' LOBE. TMP’ .TYPEs’OLD’ .RECORDS  I ZE>1 01 , 

2  ACCESS* 'DIRECT' .SHARED) 

C 

4991  IFI ISTATS.EQ.O)  GO  TO  4990 

C 

C  GRAPH  HIGHER  ORDER  BETA  STATISTICS  IF  REQUIRED 

C 

CALL  HOSBD 1 1 , DUMMY , IOE , DUMMY , SLOBE , BUF . NB , 

1  FREQ, POLVHI I F IX IPOLARI > .PERSH.PITCH, ROLL. YAW. TILT. PHR) 

C 

C  CLEANUP  GRAPHS 

C 

4990  IFIGDAUTO. AND. IZETA. EQ.O. AND. DEBUG. EQ.O)  CALL  QRTEKQ 

IFIDEBUG.EQ.1 .AND. COMBIN.EQ.O. AND. .NOT. ATTN)  CALL  GRTEKG 
4999  IF  I IZETA. EQ.O)  CALL  GRPLOTIO. ,0.4,0) 

IFI IZETA. EQ.O)  CALL  GRMODEI ’ALPHA’ ,3) 

IFI IZETA. GT.O. AND. .NOT. MGRAPH)  CALL  QRPLOTi 1 1 . .0 . 5 , -1 1 
IFI IZETA. QT.O. AND. MQRAPH)  CALL  GRUPPN 
IFICMD.EQ.CMDSI9) )  GO  TO  3 
4990  COMB  I  Ns  I ABS I COMB  IN) 

IF (ATTN)  CALL  DETACH! IOR, IER) 

CMDsCLEAR 

CALL  GRPLOTIO., 0.4,0) 

CALL  GRMODEI 'ALPHA' .3) 

GO  TO  70 

4944  READ  1 1 OF, 2)  BUF 
LINEsLINE*1 
GO  TO  3 
C 

C  -PRINT''  -  PRINT  RESULTANT  RCS  VECTOR 

C 

90  CONTINUE 

I F I DEBUG . QE . 1 . AND . CMD . EQ . CMDS 1 10))  GO  TO  4944 
10= IOW 

C  FORCE  OUTPUT  TO  LINE  PRINTER  FOR  DEBUG 
IPCODEsARRM) 

IFIARRID.GE.I .  O.OR.CMD.  NE.  CMDS  1 101)  10*  IOP 

ANQSTHETAI1I 

ANQINCSTHETAI3) 

LP*. FALSE. 

LPFs. FALSE. 

IFIARRI1) .GE. 2.0. AND. DEBUG. EQ.O)  LPs.TRUE. 

IFIARRI1 ) .QE. 3.0. AND. DEBUG. EQ.O)  LPFs. TRUE. 

IF  (LP)  CALL  ASNLUNIIOP.’SY’ ,0) 

IF (LP)  OPEN  I  UN IT> IOP .NAMEs’ LP . LST’ ) 

IF  1 10. EQ. IOP. AND. .NOT. (GOULD. OR. LP))  CALL  ATTACH! 10, IER) 

IF  1 10. EQ. IOP. AND. .NOT. LP)  GOULD*. TRUE. 

IF  (CMD.EQ.CMDSMO)  I  GO  TO  9090 
C 

C  DEBUG  OUTPUT  PRINT  MODE 
MOODBsl ./10««9 
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4790* 

DO  5054  <1*1  ,  NRCS 

V 

4800> 

IF (RCSIJ) . GT.M90DB)  00  TO  9053 

481 0> 

9054 

CONTINUE 

4820> 

WHITE! 10.9099)  CMD.LINE 

4830> 

5099 

FORMAT  CO*  ,  A4 ,  ’  HAS  NO  RCS  CONTRIBUTION  (DATABASE  LINE’ , 13, ’ 1 ’ 1 

4840> 

00  TO  5056 

4850> 

C 

4860> 

9053 

CONTINUE 

4870> 

WRITE! 10.531  CMD.LINE 

m 

4880> 

59 

FORMAT (/’ O’ ,A4,’  COMPONENT  OUTPUT  FROM  DATABASE  LINE ’,131 

4880> 

00  TO  5051 

V 

4800> 

5050 

READ  1 1 OF , 2 , END*9993)  BUF 

481 0> 

LINEUIMEM 

4920* 

CALL  TRIM80INB.BUF) 

/*• 

4930* 

WRITE  1 10,54)  (BUF<l),|s1,NB) 

4840> 

94 

FORMAT  1 / ' 0 ' ,1 20A1 1 

» '  ’ 

4850» 

FREQs. 3/WLEN 

4880> 

WRITEI 10,5057)  FREQ 

4870> 

5057 

FORMAT!'  FREQUENCY IflHZJs* ,F5 . 1 1 

4880> 

IBQsIBUG. AND. 1 

4880» 

IF (CMD.NE .CMOS  1101)  IBQ=1 

9000> 

C 

IFMSTATS.EQ.il  IBQsO 

5010» 

5091 

IF! IRPR.EQ.O)  WRITE  1 10,51 ) 

9020* 

51 

FORMAT  1 ’ 0  RCS  OUTPUT’ /2X, ' ANQLE’ ,5X, ‘RCS’ ,  6X, 

5030* 

9  ’  RCS-S ' , 5X , ’ RCS  +  S ’ /I XI 

Ml 

S040> 

IF(IRPR.EQ.I)  WRITEI 10,571 

5050* 

57 

FORMAT  I ’ 0  RPRCS  OUTPUT  < /2X , ’ ANGLE ' , 2X , ' RCS 1 DBSMI ’ , 2X , 

5080> 

1  ' RCS I  SMI ’ I 

S070> 

C 

.V 

5080» 

C 

9080> 

DO  52  Jsl.NRCS 

9100> 

IF  I IRPR.EQ.O)  00  TO  5058 

9 1 1 0> 

YsRCSI  JIN<l2«RCS2 1 JIMM2 

*5 

91 20> 

YYY310.HAL0010IAMAX1 (CUTOFF, YD 

9 1 3  0  > 

WRITE! 10,5059)  ANO.YYY.Y 

91 40> 

9058 

FORMAT IF7.1.2F1 0.4) 

9190> 

00  TO  59 

91 80> 

C 

:v 

9170> 

5058 

RCSCsAMAXI (CUTOFF , RCS (J ) ) 

91 80> 

SsSQRTIAMAXI 10., RCSIJ )««2-RCS2 IJ ) ) ) 

91 80> 

0BSM*1 0 .ttALQGI 0 IAMAX1 (CUTOFF ,RCSC) ) 

MV 

9200> 

DBSMMS*10.»ALOG10(AMAX1 (CUTOFF, RCSC-SD 

N  • 

9210> 

DB’MPS*1 0 . MAL0Q1 0 1 AMAX1 (CUTOFF ,RCSC*S) ) 

9220> 

WRITEI 10,53)  ANQ , DBSM, DBSMMS , DBSMPS 

9230> 

89 

FORMAT  I F7.1 ,3F10.4,3X,2F10.4) 

'  „ 

9240> 

59 

IFI ISTATS.EQ.O)  00  TO  52 

y 

9290> 

IFISTATS.QE.ANQINC. AND. AMODIANQ, STATS). NE.O)  00  TO  52 

9280> 

C 

9270> 

C 

COMPUTE  HIQHER  ORDER  STATISTICS  FOR  BETA  DISTRIBUTION 

9280> 

c 

hosbd  DOES  NOT  USE  STATS  INTERNALLY  WITH  PRINT 

.  - 

9280> 

c 

9300* 

CALL  HOSBD 1 0 , J , 1 OE , 1 0 , 8L0BE ) 

** 

9310> 

c 

9320> 

52 

ANQ*ANQ*ANOINC 

«,  - 

9330> 

IF (DEBUG. NE.O)  00  TO  5068 

9340» 

c 

o 

9390> 

c 

COMPUTE  AND  PRINT  LOBE  STATISTICS  IF  REQUIRED 

9380> 

c 

5370* 

IFIKLOBE.LE.OI  00  TO  5590 

9380* 

IQRPRsO 

HI 


■  t 
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mo> 

WR 1 TE II OL  *  1 I  1 QRPR , KLOBE .GLINT, NRCS , NLOBE , THETA .ALPHA 

9400> 

1  BE  TA , SC ALE , 1 ZE  TA . 1 PCODE , WL  E  N , L  BUG 

$410> 

CLOSE (UNI T* IOLI 

9420> 

CALL  SPAWNS 

S430> 

OPEN ( UN  1 T* 1 OL , NAME  * ’ LOBE . TMP  * . T YPE  * 1  OLD ' , RECORDS  1 ZE  *1 

9449> 

2  ACCESSs’DIRECT’ .SHARED) 

9490> 

C 

9490> 

9990 

IF (GOULD)  WRITEI 1  OP, $052) 

9470> 

9092 

FORMAT  1 ' 1 • ) 

9490> 

9099 

IF (GOULD. AND. DEBUG. EQ.O)  CALL  DETACH! 10, IER) 

9460> 

1 F IQOULD . AND . DEBUG . EQ . 0)  GOULD* . FALSE . 

9900> 

IF (LP .AND. . NOT.LPF)  CLOSE (UNIT* IOP , D ISPOSE*' PRINT’ ) 

991©* 

IF ILP .AND . LPF)  CLOSE (UNIT* IOP .DISPOSE*’ SAVE ’ ) 

9920> 

IF (LP)  CALL  ASNLUNI IOP, ’QD’ ,0) 

99SQ> 

1 F 1 CMD . EQ . CMOS (10)1  GO  TO  3 

9940> 

CMD*CLEAR 

9990> 

GO  TO  70 

9990> 

C 

9970> 

C 

"PRYT"  -  MASTER  ROTATION  ANGLES 

9990> 

C 

99»0> 

99 

CONTINUE 

S600> 

222 

FORMAT!  'O'/MX.SFS  .31) 

9610> 

C 

COMPUTE  YAW* > RE SUIT 

9<29> 

PRYTV(3)>ARR(3) 

9690> 

C*COS ( ARR ( 3 1 NRAD I 

9640> 

S*SINIARR(3)NRAD) 

9«90> 

CALL  MATPUTIPRVT.1 . .0. ,0. .0. .C.-S.O. ,S,C) 

9990> 

C 

MULTIPLY  BY  R0LL*>ROLLNRESULT*RESULT 

9970> 

PRYTV(2)*ARR(2I 

9660> 

C*COS ( ARR 1 2 ) *  RAD ) 

96I0> 

S>8INIARR(2)ttRAD) 

9700> 

CALL  MATPUT IMAT.C, -8,0. ,S,C,0.,0.,0.,1.) 

9710> 

CALL  MATMUL (MAT, PRYT. PRYT) 

9720> 

c 

MULTIPLY  BY  P ITCH«>P ITCHNRESULTsRESULT 

9790> 

PRYTVI1 )*ARR(1 ) 

9740> 

C*COS(ARR(1 INRAD) 

9790> 

S>SIN(ARRI1)nRAD) 

9790> 

CALL  MATPUT (MAT ,C,0. ,S,0. ,1 . ,0.,-B.O. ,C) 

9770> 

CALL  MATMUL (MAT. PRYT. PRYT) 

97f0> 

c 

MULTIPLY  BY  TILT*>RESULTNTILT*RESULT 

9790> 

PRYTV(4)*ARRI4) 

9IQ0> 

C*COS(ARR(4)NRADI 

9910> 

S>SINIARR(4)NRAD) 

9i20> 

CALL  MATPUT (MAT, C.O. ,-S,0. .1 . ,0. .8,0. , Cl 

9990> 

CALL  MATMUL (PRYT, MAT, PRYT) 

9940> 

c 

MULTIPLY  BY  B*>BNRESULT*RESULT»PRYT 

9«90> 

CALL  MATPUT (MAT, 0. ,-1 .,0..1.,0.,0.,Q.,0..1.) 

9990> 

CALL  MATMUL (MAT .PRYT, PRYT) 

9I70> 

BO  TO  3 

9»I0> 

c 

9S60> 

c 

"PAUSE"  -WAIT  FOR  A  CARRIAGE  RETURN 

9900> 

c 

9910* 

60 

CONTINUE 

S620> 

READ! 1 OR ,61 ,END*6090)  J 

9630> 

61 

FORMAT (A1) 

S640> 

GO  TO  3 

96S0> 

C 

9660> 

C 

"RAMDB”  •  MODIFY  COMPONENT  RCS 

9670> 

C 

9990> 

69 

RAM*10.NN(ARR(1)/10.J 
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9910 

9000 

6010 

6020 

9030 

6040 

6090 

6060 

6070 

6060 

6090 

6100 

6110 

6120 

6130 

6140 

6190 

6160 

6170 

6160 

6190 

6200 

6210 

6220 

6230 

6240 

6290 

6260 

6270 

6260 

9290 

6300 

6310 

6320 

6330 

6940 

6990 

6360 

6370 

6960 

6390 

6400 

6410 

6420 

6490 

6440 

6490 

6460 

6470 

6460 

6490 

6900 

6910 

6920 

6690 

6940 

6990 

6960 

6970 

6960 


00  TO  3 
C 

C  "CLEAR”  -  CLEAR  RCS  VECTOR 

C 

70  CONTINUE 

IPASSsO 

OBJECT:. FALSE. 

NVALIDsO 

NTOTALsO 

DO  71  Jsl.HAXRCS 
RCS2IJI=0. 

RCSUJsO. 

YYVIJ ,1 )®0. 

YYVIJ.2IS0. 

XVIJ)=0. 

C  EXTRA IJ 1=0. 

I F IOPENE I  WRITE* I OE’ J ) DZERO , DZERO , DZERO , DZERO , DZERO , DZERO 

71  CONTINUE 

IFICMD.EQ. CLEAR)  00  TO  3 
C 

C  RESET  SYSTEM  VECTORS 
C 

POLAR=1 . 

C 

C  SET  PYRT  TO  IDENTITY  BY  DEFAULT 
C 

DO  74  J«1,9 

74  PRYT(J)slDENT)J) 

C 

C  RESET  LOBE  PARAMETERS 

C 

KLOBE=0 

LL0BE*0 

MLOBEtO 

NLOBEs-3 

BETA*0. 

C 

C  RESET  RPRCS  PARAMETERS 

C 

IRPRsO 

C 

C  RESET  VALID  VECTOR  AFTER  EACH  OBJECT 

C 

71  VALIDI1 )*-360. 

VALIDI2H360. 

DO  79  J*9,6 
79  VALIDlJItO. 

RAM  si . 

DO  79  J«1 ,9 

79  TRANS IJ)> I  DENT) J) 

DO  76  J«1,3 
79  ORQNI J I *0 . 

IFIIPASS.EQ.il  00  TO  99 
00  TO  9 
C 

C  DEDUO  OUTPUT 

C 

90  IBUQbIF  IXIARRM  1 1 

LBUO* IF IXIARRI2) I 
00  TO  9 


;■«  ^ 


1 


£ 
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6590»  C 
6600>  C 
661 0>  C 
6620>  66 
6690> 

6640> 

6650>  66 
6660> 

6670> 

6660> 

6690> 

6700> 

6710>  C 
6720>  CNN* 
6730>  C 
6740>  C 
6750>  C 
6760>  C«N» 
6770  C 
6760  C 
8760>  C 
6«00>  500 
6610>  C 
66 20>  C 
6630> 

6640> 

6tS0> 

8660 

6670> 

8160 

6800> 

8600 

6610> 

6920  C 
69»0>  C801 
8940  C502 
6950  C 
6960  C 
6970  C 
6960  C 
699  0  610 
7000> 

7010 
7020>  C 
7030  C 
7040  C 
7050  520 
7090> 

7070> 

7060> 

7090>  C 
7100>  C 
7110>  C 
7120>  540 
7190>  541 
7140 
7180 
7160 
7170 
7160 


NEW  WAVE  LENGTH 
IWL= IWL+1 

IFIIWL.EQ.10)  00  TO  9990 
OLDQRFs. FALSE. 

REWIND  I  OF 

IFIKLOBE.QT.O)  CLOSE (UNIT* IOL) 
IF  I GOULD I  CALL  DETACH  I IOF , I ER) 
GOULDs. FALSE. 

L I  HE  =0 
GO  TO  70 


THE  FOLLOWING  COMMANDS  ARE  NEW 
ADDITIONS  TO  RCS  FOR  THE  LOBE 
STATISTICS 

IMHHHtftttltNIHHHHHIftNttNttOMMttttlHtMNNMftllMMtt 
"LOBE”  -  SET  LOBE  SWITCH 
CONTINUE 

IF  I IRPR.EQ. 1 .AND. ARR (1 ) .EQ. 0. )  GO  TO  3 
IF  1 1  STATS.  EQ.  1 .  AND.ARRH  1  .  QT .  0 . )  GO  TO  673 
LL0BE*ARR(1) 

IF  I LLOBE  .GT  .0 .  AND  .KLOBE .  LE  .  01  CALL  NOLOCKIIOL) 

IF ILLOBE .QT . 0 .AND .K LOBE . LE . 01 

1  OPEN  I UN ITs IOL, NAME >’ LOBE. TMF* .TYPE* ’UNKNOWN’ , RECORDS  I ZE* 191 , 

2  ACCESS*' DIRECT’ .SHARED) 

IFILLOBE.GT.O. AND. KLOBE. LE.O)  WRITEIIOL’1)  DUMMY 
IF ILLOBE. QT.OI  MLOBE* LLOBE 
IFILLOBE.GT.O)  KLOBE*LLOBE 
GO  TO  3 

WRITE  1 2, 502 1 

FORMAT  I’  RC8>  —  RPRCS  AND  LOBE  CANNOT  BOTH  BE  ACTIVE’) 

GO  TO  9993 

’’BETA”  -  DEFINE  BISTATIC  ANGLE 
BETASARRI1I 

IFIBETA.LT. S.)  BETASO. 

GO  TO  3 

“ORON"  -  COMPONENT  ORIGIN  VECTOR 

ORGN 1 1 ) »ARR 1 1 1 
ORQN 1 2 ) *ARR 1 2 1 
ORQN 1 3 ) ■ ARR 1 3 ) 

GO  TO  3 

"OELT"  -  FIRST  LOBE  ANGULAR  RANGE 

DO  541  4*1 ,4 
DELTA 14) *ARR IJ) 

DELTA  12 )*AMAX1 (DELTA! 1 ) .DELTAI2) ) 

DELTA  13 )»AMAX1 (1 . , DELTA  1 3 1 1 

IF (DELTA (4) . LT. 4 .0. OR. DELTA  14) .GT. 160. )  DELTA(4)*49. 

GO  TO  3 
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71 90> 

C 

7200> 

C 

7210* 

5S0 

7220> 

7230> 

SSI 

7240> 

C 

7250> 

c«ft»< 

7260> 

c 

7270> 

c 

72B0> 

c 

7290> 

c«««! 

7300> 

c 

7310> 

c 

7320> 

c 

7330> 

560 

7340> 

C 

73S0> 
7300> 
7370> 
73S0> 
7390> 
7400> 
741 0> 
7420> 
7430> 
7440> 

C 

7450> 

C 

7490> 

C 

7470> 

C 

7490> 

570 

7490> 

C 

7S00> 
731 0» 
7520> 
7530> 
7540> 
7SS0> 
7S«0> 
7570> 

C 

7590> 

7590> 

C 

7600> 

7910> 

575 

7920> 

C 

7630> 

CS71 

7640> 

C972 

76S0> 

C 

7B60> 

CC 

7970> 

C97I 

7090> 

CS74 

7990> 

C 

7700> 

C 

7710> 

C 

7720> 

C 

7730> 

7740> 

5S0 

77S0> 

C 

7760> 

C 

7770> 

C 

7790> 

590 

"ALPH"  -  SECOND  LOBE  ANGULAR  RANGE 

DO  SSI  <1*1,4 
ALPHA CJ) SARA <J) 

GO  TO  3 


END  OF  LOBE  COMMANDS 


"RPRCS"  -  RELATIVE  PHASE  RCS 
CONTINUE 

IFIMLOBE.OT.O. AND. ARAM)  .8T.O.)QO  TO  SOI 
IFHSTATS.EQ.1. AND. AAA<1)  .QT.O.)  GO  TO  S71 
I RP R  =M I  NO ( 1 ,MAX1 (0 . , ARR1 1 1 ) I 
IF  I IRPR .EQ.OI  GO  TO  3 
RPRCSG*ARR(2) 

LRPRCS=ARR 13 ) 

IFILRPRCS.LE.il  LRPRCS*5 

LOBE *1 

LLOBEst 

MLOBEsI 

GO  TO  3 

"STATS”  -  HIGHERER  ORDER  STATISTICS 
CONTINUE 

I F 1 1  RPR . EQ . 1 . AND . ARR 1 1 1 . GT . 0 . 1  GO  TO  571 
IFIMLOBE.OT.O. AND. ARRI1I.QT.0.IG0  TO  S73 
I STATSsMI NO (1 ,MAX1 (0. , ARR (till 
STATS*ARR(2) 

IF (STATS . LE .0 . I  STATS*10. 

IFI I  STATS . EQ. 0 . OR .OPENE I  GO  TO  S 

OPEN < UN  I T= I OE, NAME  = ’ACSTMP.DE S’ .TYPE*’ SCRATCH’ , 

1  ACCESS* ’DIRECT' ,REC0RDSIZE*1S,MAXREC*1S1 I 

OPENE*. TRUE. 

DO  STS  J*1  ,111 

«R I TE ( IOE ’ J I DZERO , DZERO, DZERO , DZERO , DZERO, DZERO 
GO  TO  3 

WRITE(2,S72I 

FORMAT  I ’  ACS>  —  RPRCS  AND  STATS  CANNOT  BOTH  BE  ACTIVE’) 
GO  TO  BB93 

VRITE(2,S74) 

FORMAT!'  RCS>  —  LOBE  AND  STATS  CANNOT  BOTH  BE  ACTIVE’) 
GO  TO  9903 

"GLINT"  -  SET  GLINT  SWITCH 

QLINT*MIN0 12 ,MAX1 10. ,ARR(1)I) 

GO  TO  3 

"NEWTH”  -  REPEAT  LOOP  FOR  NEW  THETA  ANGLE  RANGE 
CALL  CNEWTHIKLOBE .ROUTE) 
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7T90> 

7600> 

C 

7610> 

C 

7620> 

C 

7630> 

600 

7B40> 

7650> 

7960> 

601 

7670> 

C 

7690> 

C 

7B90> 

C 

7900> 

7910> 

7920> 

7930> 

90 

7940> 

91 

79S0> 

7960> 

C 

7970* 

7990> 

92 

7990> 

C 

I000> 

c 

•  010» 

c 

9020> 

c 

9030> 

c 

S040> 

c 

6090> 

c 

9060> 

c 

9070> 

100 

90S0> 

8090> 

c 

B100> 
61 10> 
6120> 
61 30> 
6140> 
9150> 
•  160> 
9170> 
61 90> 
6190> 

99 

9200> 

6210> 

9220> 

6230> 

6240> 

c 

•  250» 

c 

6260> 

c 

6270> 

c 

S290> 

c 

9290> 

c 

9300> 

c 

6310> 

c 

6320> 

c 

6330> 

c 

6340> 

c 

6360> 

c 

6360> 

c 

6370> 

9360> 

101 

00  TO  (3,96,9993) .ROUTE 

"GSIZE*  -  SET  GRAPH  SIZE 

DO  601  J*1,5 
QSIZE(J)*ARRIJ) 

CALL  QRPICIQSIZEI1 ) ,GSIZE(2) ,05 1 ZE 131 , GSIZE 1 4)1 
00  TO  3 

DETERMINE  IF  COMMAND  IS  AN  OBJECT 

CONTINUE 

DO  91  J*1,N0BJ 
IOBJ*J 

IF (CMD . EQ.OBJSI J ) )  00  TO  100 
CONTINUE 

WRITE  I IOV, 92)  CMD 

FORMAT ( '  RCS>  —  INVALID  FUNCTION  ,A4, ’“’ ) 

00  TO  9993 


ADD  OBJECT  CONTRIBUTION  TO  RCS  VECTOR 

AIMHIlHHtltHIMMHIItttllltPHIIMNttllMIMtiHHHHHINNIHt 


I F IWLEN . LE . 0 . )  00  TO  9991 

NOTE:  IPASS  IS  ONLY  USED  WITH  RFLAT  IN  CONJUCTION  WITH  LOBE 
IPASS*-1 
IPASS* IPASS+1 
OBJECT*. TRUE. 

ANQ*THETAI1) 

ANQMAXiTHETA (2 14.001 
ANQ I  NOTHETA  (3) 

CPHR*COS I PHR*RAD ) 

SPHR*SIN(PHR*RAD) 

IRCS*1 

CALL  TRANSP ( PRYT , PRYT2 ) 

CALL  MATMUL I PR YT2 , TRANS , COMP ) 

IF  1 1 IBUQ.AND.2) .EQ.O)  00  TO  101 
WRITEI2.222)  TRANS 
WRITEI2.222)  PRYT 
WRITEI2.222)  COMP 

ANQ  IS  RADAR  THETA  ANGLE 
PHR  IS  RADAR  PHI  ANQLE 

IF  END  OF  ANGLE  RANGE  FOR  OBJECT.  OUTPUT  PARTIAL 
LOBE  DATA  IF  REQUIRED,  THEN  SET  NEXT  COMMAND. 

KLOBE»SET  FROM  LOBE  ENTRY  POINT 
MLOBEoLOBE  IS  BET  FOR  RUN 
LLOBE»LOBE  IS  SET  FOR  OBJECT 
LOBE *> LOBE  IB  BET  FOR  CURRENT  THETA 

IF  I IRCS.LE.NRC8)  SO  TO  131 

IF( I RPR. EG. 1 . AND. KLOBE. EQ.O)  00  TO  136 
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1 
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i 

i 

i 

i 


! 

> 

i 

l 

■ 

f 
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9390 

9400 

9410 

9420 

9430 

9440 

9430 

9460 

9470 

9490 

9490 

9900 

9510 

9520 

9530 

9540 

9550 

9560 

6570 

9590 

9590 

9600 

9610 

9620 

9630 

9640 

9050 

9660 

8670 

8690 

8690 

8700 

9710 

8720 

8730 

8740 

8750 

8760 

8770 

8790 

8790 

9900 

8910 

8920 

9930 

9940 

9850 

9960 

9970 

9980 

9990 

9900 

9910 

8920 

9930 

9940 

9990 

9900 

9970 

9990 


i 

i 

I 

i 


IFILL0BE.LE.0)  Q0  TO  72 
C  OUTPUT  LOBE  INFO  TO  TEMPORARY  FILE 
NLOBE*NLOBE*4 
DO  139  J*1 ,4 

139  WR I TE I I0L ’ NLOBE* J )  (SLOBE I J ,1,11,1=1 , NRCS) 

DO  133  I  si , NRCS 

IF  I SLOBE (1 , 1 ,2) .NE .0. )  80  TO  134 

133  CONTINUE 
80  TO  135 

134  NLOBEsNLOBE+4 
DO  137  J*1 ,4 

137  WRITE! I0L ' NLOBE+ J )  (SL0BE1J, I .2) , 1=1 .NRCS) 

C 

C  MAKE  SECOND  PASS  THROU8H  RFLAT  FOR  LOBE 

135  IF (CMD .EQ.OBJS ( 13) .AND. I  PASS . EQ.O)  IPASS=1 
80  TO  72 

C 

C  DETERMINE  IF  VALID  ANGLE  FOR  LOBE  CONTRIBUTION 

C 

131  LOBESLLOBE 

IF (LOBE . LE .0)  80  TO  132 
SLOBE 1 1 , IRCS.1 1*0. 

SLOBE ( 1 , I RC  S , 2 ) =  0 . 

C 

C  SKIP  ANGLE  IF  NOT  VALID  FOR  CURRENT  OBJECT 
C 

132  DO  102  Jsl.5,4 

IF (VALID(J) -VALID! J*1)  I  103,102,104 

103  IFIAN8.LT. VALIDIJ) .OR . ANG .GT . VAL I D I  J*1 ) )  80  TO  102 

IF (VAL IDI J+2 I -VAL ID( J+3) I  105,109,106 

105  IFIPHR.LT.VALI01J+2).OR.PHR.GT.VALIDIJ+3))  GO  TO  102 

80  TO  109 

106  IFCPHR.GT.VALID(J*3J.AND.PHR.LT.VALIDCJ*2JJ  80  TO  102 

80  TO  109 

104  IF(AN0.8T.VALIDIJ«1) .AND.ANO.LT. VALIDIJ))  80  TO  102 

IFIVALID(J*2)-VALIDCJ*3>)  107,109,108 

107  IFCPHR.LT. VALIDCJ+2) .OR . PHR .GT . VAL I D ( J*3 ) )  80  TO  102 

80  TO  109 

108  IF(PHR.GT.VALIDIJ*3) .AND. PHR. LT. VALID! J*2) )  80  TO  102 

80  TO  109 
102  CONTINUE 
C 

C  CURRENT  ANGLE  IS  NOT  VALID.  INCREMENT  ANGLE. 

C 

GO  TO  130 
C 

C  CURRENT  ANGLE  IS  VALID.  FIND  PRIMED  ANGLES  FROM  RADAR  COORDINATES. 
C  NOTE!  THIS  PROGRAM  STORES  MATRICES  BY  ROWS  (NOT  COLUMNS) 

C 

109  CONTINUE 

CANG*COS I ANGuRAD ) 

SANG*  S I N I ANGttRAD I 

CTHP*COMP  (3)«SANG«CPHR*C0MP  <6)»$ANGMSPHR*C0MP  tBUCANG 
THP*ACOSK (CTHP) 

STHP*SIN(THP) 

CPHPS* ICOMP 1 1 ) »SANG«CPHR*C0MP 1 4 ) NSAN8ttSPHR*C0MP 1 7 ) «CANG ) 

SPHPS* ICOMP 1 2 1 NSANGnCPHR+COMP 1 5 ) »SANG«SPHR*COMP 1 8 IMCANG) 

PHP*0. 

IF ( AMS  I SPHPS ) . LT. 1 E-S . AND. ABS (CPHPS) . LT . 1 . E-S )  80  TO  1111 
PHP*ATAN2 I SPHPS, CPHPS) 


38 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 
laurel.  Maryland 


8990> 
9000> 
901 0> 
9020> 

1111 

9030> 

223 

9040> 

C 

90S0> 

C 

9080> 

c 

9070 

c 

9080> 

c 

9090 

c 

9100> 

c 

9110 

c 

9120> 
9130 
91 40> 

c 

91 50> 

c 

9160 

C  F 

9170> 

9180 

9190> 

9200 

9210> 

9220> 

9230 

9240> 

9250> 

C 

9260 

9270> 

9290> 

9290 

9300> 

9310 

1112 

9320> 

C 

9330 

C 

9340> 

C 

9380 

9360> 

110 

9370 

C 

93S0> 

C 

9390> 

C  • 

9400 
941 0> 

200 

9420> 

C 

9430 

C  ” 

9440> 

C 

9480> 

9460 

210 

9470> 

C 

94S0> 

C  ” 

9490 

9S00> 

220 

9810> 

C 

9820> 

C  ” 

9830 

9S40> 

230 

9880 

C 

9S60> 

C  " 

9870> 

C 

9880 

240 

C2PHP*CPHP#CPHP 
SPHP-S I N I PHP ) 
S2PHP3SPHPNSPHP 


AOD  OBJECT  CONTRIBUTION  FOR  CURRENT  ANGLE 


RCSC-0 

IFMOBJ.LE.O.OR.  IOBJ.QE.10J  GO  TO  110 

FIND  UNPRIMED  ANGLES  FOR  OBJECTS  WHICH  USE  POLARIZATION 

CTH&PRYT2 (3)*SANG*CPHR*PRYT2 (6)MSANG<tSPHR+PRYT2 I9)*CANG 
TH-ACOSKICTH) 

Sth=S I N I TH ) 

CPHS= IPRYT2 (1 )*SANGl»CPHR+PRYT2 (4)*SANG*SPHR+PRYT2 (7 )*CANG) 
SPHSs ( PRYT2 1 2  HlSANG»CPHR*PRYT2 1 S )«SANG«SPHR*PRYT2 ( 8  UCANG) 
PH=0. 

IF IABS I SPHS) . LT . 1 .  E-S .AND . ABS (CPHS) . LT . 1 . E-5)  GO  TO  1112 
PHsATAN2(SPHS,CPHS) 

2  CPH=CO$(PH) 

epu*Q I  fel  I PH I 

I F I (I BUG . AND . 4 1 . EQ . 0 )  GO  TO  110 

DTHaTH/RAD 

DPHsPH/RAD 

WRITE(2,223)  DTH.DPH 
GO  TO  SPECIFIED  OBJECT 

GO  TO  1200, 210, 220, 230, 240,250, 280, 270. 280, 290,300, 310, 320, 330. 
S  340,380,380,370),  IOBJ 


OGIVE”  -  OGIVE  A  TRUNCATED  OGIVE 
CALL  OGIVE 
GO  TO  120 


INPUTSsA.B.L 
CALL  CYLIN 
GO  TO  120 


CALL  CONE 
GO  TO  120 

PARAB”  -  PARABOLA 
CALL  PARAB 
GO  TO  120 

SOLID”  -  ELLIPSIOD  A  SPHEROID 
INPUTS*A,B,C 
CALL  SOLID 
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9590* 

ao 

TO  120 

»eoo> 

c 

981 0> 

C  "TORUS" 

•  TORUS 

9820* 

290  CALL 

TORUS 

9630> 

QO 

TO  120 

9940> 

C 

9650> 

C  "WEDGE" 

-  WEDGE 

9660> 

280  CALL 

WE DOE (TRANS, CTH.STH, CRH, SPH) 

9670> 

ao 

TO  120 

9880> 

C 

9690> 

C  "WIRE- 

WIRE 

9700> 

270  CALL 

WIRE (TRANS, CTH.STH, CRH, SPHI 

971 0> 

ao 

TO  120 

9720* 

C 

9730> 

C  "LWIRE” 

-  LOOP  WIRE 

9740> 

280  CALL 

LWIRE (TRANS, CTH.STH, CPH, SPH) 

9750> 

00 

TO  120 

9760> 

C 

9770> 

C  "CAVA" 

CAVITY  AREA 

9780> 

290  CALL 

CAVA 

9790> 

00 

TO  120 

9900> 

C 

981 0> 

C  "CAVD" 

-  CAVITY  DIAMETER 

9820> 

300  CALL 

CAVD 

9890> 

GO 

TO  120 

9840> 

C 

9850> 

C  "CFLAT" 

-  CIRCULAR  FLAT  PLATE 

9860> 

310  CALL 

CFLAT 

9870> 

ao 

TO  120 

9880> 

C 

9890> 

C  "RFLAT" 

-  RECTANGULAR  FLAT  PLATE 

9900> 

320  CALL 

RFLAT ( IPASS) 

991 0> 

IF  1  IPASS.EQ.2)  RCSCsO. 

9920> 

00 

TO  120 

9930* 

C 

9940> 

C  "TWAVE" 

•  TRAVELING  WAVE 

9950> 

330  CALL 

TWAVE 

9960> 

00 

TO  120 

9970> 

C 

9990> 

C  -DIKED" 

•  DIHEDRAL 

9990> 

340  CALL 

DIHED 

10000> 

00 

TO  120 

10010> 

C 

1 0020> 

C  "TCORN" 

•  GENERAL  TRIANGULAR  CORNER 

1 0030> 

380  CALL 

TCORN 

1 0040> 

00 

TO  120 

1 0050> 

C 

iooeo> 

C  "ACORN” 

•  GENERAL  RECTANGULAR  CORNER 

10070* 

380  CALL 

RCORN 

10090* 

00 

TO  120 

10090* 

C 

10100* 

C  "ECORN” 

-  GENERAL  ELLIPTICAL  CORNER 

10110* 

370  CALL 

ECORN 

10120* 

00 

TO  120 

10130* 

C 

10140* 

C 

10180* 

C  ADO  RCS  CONTRIBUTIONS  TO  VECTORS 

10180* 

C 

10170* 

120  CONTINUE 

10190* 

IFIRCSC.LT. 0.1  WRITE! I0W.121 )  CMD.ANQ.RCSC 
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m 


10190 

10200 

10210 

10220 

10230 

10240 

10250 

10200 

10270 

10200 

10290 

10300 

10310 

10320 

10330 

10340 

10350 

10300 

10370 

10300 

10390 

10400 

10410 

10420 

10430 

10440 

10450 

10400 

10470 

10400 

10490 

10500 

10510 

10520 

10530 

10940 

10550 

10500 

10570 

10500 

10590 

10000 

10010 

10020 

10030 

10040 

10050 


FORMAT!’  RCS>  —  ’ ,A4,F7 .1 , E 1 5 .01 

IF (RCSC.GT. CUTOFF )  NVAL ID*NVAL 10*1 

RCSORCSCNRAM 

IF ( IRPR . EQ. 1 )  00  TO  130 

RCS!  IRCSIsRCS! IRCSI+RCSC 

RCS2 1 1 RCS) *RCS2 ( I RCS ) + ACSC*RCSC 

IF ( ISTATS.EQ.O)  00  TO  130 

READ! I 0E ’ IRCS)  T1 ,T2 , T3 ,T4,T5 ,T6 
T1 *T1 +DBLE IRCSC) 

T2*T2+DBLE  (RCSONM2 
T3*T3tDBLE(RCSC)»«3 
T4ST4+DBLEIRCSCMHI4 
SQ=SQRT IRCSC) 

T5ST5+SQ 
T6=DMAX1 (TO, SO) 

WRITE IIOE’ IRCS)  T1 ,T2 ,T3 ,T4 ,T5 ,T0 

INCREMENT  ANGLE 

CONTINUE 

ANQ*ANQ»ANQINC 

i sr  c  s i rcca i 

I F  ( I  PASS  .NE .2 )  NT0TAL=NT0TAL*1 
00  TO  101 


WRITE! IOW, 9992) 

FORMAT!’  RCS>  —  WAVELENGTH  UNDEFINED’) 

WRITE! IOW, 9994)  LINE 

FORMAT!’  RCS>  --  LAST  DATA  BASE  LINE  READ  I  S’. 1 4) 
CONTINUE 

TERMINATE  RCS  PROGRAM 

I F IKLOBE . QT . 0)  CLOSE (UN I T* IOL , D I SPOSE* ’  DELETE ’ ) 

IF (OPENS)  CLOSE ( UN  IT* I OE) 

IF (GOULD)  CALL  DETACH! I0P , IER) 

T IME*SECNDS (T IME) /SO . 

IF! IZETA.EQ.O)  CALL  QRMODE! ’ALPHA’ ,31 
WRITEI2.9999)  TIME 

FORMAT!*  RCS>  --  ELAPSED  TIME  s’.FS.I,’  MINUTES’) 

CALL  EXIT 
END 
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»>  FILE  a  SYO: C 1 1 1 , 1 21RCSCMD . FTN  27-SEP-83  14:34:38  <<< 

10>  C 

20>  SUBROUTINE  CNEWTH (KLOBE , ROUTE  I 

30>  C 

40>  INTEQERM2  ROUTE 

$0>  C 

B0>  INTEQERK2  DEBUG , COMB  I N , QL I  NT 

70>  LOG  I CALM2  GOULD, LP ,OPENE,LPF .OLDGRF , MGRAPH , QDAUTO 

•  0>  COMMON/SWITCH/  ION, THETA  14) .GOULD , LP .OPENE , LPF .OLDGRF .MGRAPH . 

»0>  1  GOAUTO, DEBUG, COMB  IN .GLINT, (BUG, LBUQ, I  STATS, I2ETA 

100>  C 

1 1 0>  ROUTE si 

120>  IF  I  THETA  1 4) .GT. THETA  12 1 1  GO  TO  1 

130>  IF  I .NOT. MGRAPH)  RETURN 

1 40>  IFI IZETA.EQ.O. AND. .NOT. QDAUTO)  RETURN 

1SO>  IFI IZETA.EQ.O)  CALL  GRTEKG 

1S0>  IFI IZETA.NE.O)  CALL  GRPLOTM 1 . ,1 . 5 , -1 ) 

1 70»  RETURN 

180>  C 

190>  1  IFIKLOBE.EQ.1 .OR. GLINT. EQ.1 .OR. ISTATS.EQ.1 .OR. DEBUG. NE.O)  GO  TO  S91 

200>  OLDGRF s. TRUE. 

21 0>  ROUTES2 

220>  RETURN 

230>  C 

240>  S91  WRITE! IOV.S92) 

2S0>  992  FORMAT!’  RCS>  —  MULTIPLE  GRAPH  SEGMENTS  ARE  NOT  ALLOWED  WITH:’/ 
260>  1  LOBE,  STATS,  OR  DEBUG!’) 

270>  ROUTE*3 

290>  RETURN 

290>  END 
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>>>  FILE  ■  SYO : [ 1 1 1 , 1 2 3RCS0BJ . FTN  27-SEP-93  14:34:44  <<< 

0>  C 

5>  C  <»X»X*XI»XI»X*X»XIIX«XNX*XI»X)|X»X»X»X»X»X«> 

10>  SUBROUTINE  OGIVE 

20>  C 

30>  C  WITH  LOBE 
40>  C 

50>  IMPLICIT  REALM  (A-Z) 

60>  INTEGERN2  LOBE 

70>  COMMON /RCSCOM/CUTOFF ,«LEN , P I , P 1 2 , RCSC , THP , CTHP . STHP , PHP , CPHP . SPHP , 

•0>  S  C2PHP.S2PHP, RAD, POLAR. LOBE, 

90>  S  R1,A,B,BETA 

100>  C 

1 1 0>  THPTsTHP 

1 20>  IFITHP.GT.PI2I  THPTsPI-THP 

130>  IF (LOBE .QT.O)  L 2  =  SQRT I R 1 *R 1  - IB  +  R1 -A)  M2) 

1 40>  IFIB.GT.O.)  GO  TO  S 

150>  3  ALPsACOSKII .-A/R1I 

1 80>  RCSCsR1«R1»SIN(ALP)»»2/(4.ttPM 

1 70>  NMA*P 1 2-ALP 

1 30>  IFITHPT.LT. NMA)  GO  TO  2 

1 90>  C 

200>  C  BO-ALP  <  THETA’  <  90 

210>  SIGMAsP IWRINRIN (1 . - (R1 -A) / (R1 hSTHPI I 

220>  IF  (LOBE .  LE  .0)  GO  TO  4 

230>  IFIRCSC.QT. SIGMA)  GO  TO  S 

240>  CALL  TOLOBEI1 .SIGMA, ANCPHP .AttSPHP ,0. ) 

245>  4  RCSC=S I GMA 

2S0>  RETURN 

290>  C 

270>  C  THETA’  *  90-ALP 

290>  6  CALP2-I1 .-A/R1 )M2 

290>  SALPsSINIALP) 

300>  ZALP=R1ttSALPtt (CALP2-SQRT(CALP2ttCALP2*SALPltSALP) I 

31 0>  WALP=ZALP*TAN(NMA) 

320>  CALL  TOLOBE I 1 , RCSC , WALPHCPHP , VALPMSPHP , ZALP ) 

330>  RETURN 

3S0>  C 

370>  C  THETA’  <  90-ALP 

3SO>  2  Z>18.KPI»CO$(THPT)M6«(1.-TANIALP)M2«TANITHPT)M2)M3 

390>  IFILOBE. QT.O. AND. Z.LT. CUTOFF)  GO  TO  S 

400>  IF IZ.LT. CUTOFF)  RETURN 

41 0>  SIQMA*WLENNWLEN«TAN(ALP)M4/Z 

420»  IFILOBE. LE.O)  GO  TO  1 

430>  IFIRCSC.LT. SIGMA)  GO  TO  9 

440>  CALL  TOLOBEI1, SIGMA, 0..0..L2) 

4S0>  1  RCSCMMIN1 1  RCSC,  SIGMA) 

4S0>  RETURN 

470>  C 

490>  C  TRUNCATED  OGIVE  (BETA  IS  ALWAYS  <  ALPHA) 

490>  C 

S00>  S  BET>BETA«RAO 

51 0>  IF (THPT.GT.P 12-BET)  GO  TO  3 

520>  C 

S30>  C  THETA’  s  0 

540>  RCSCsP  MlBHBttTAN IBET)M«2 

SB0>  IFISTHP.LT. CUTOFF. AND. LOBE. QT.O)  CALL  TOLOBE (1 ,RCSC,0. ,0. ,L2) 

5S0>  IFISTHP.LT. CUTOFF)  RETURN 

S70>  C 

590>  C  BETA  <  THETA’  <  90-BETA 
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590>  3 1 GMA=WLEN*B»TAN ( THPT  +  BET ) «M2/ ( 8 .MR IMSTHPI 

•00>  IF(LOBE.GT.O)  CALL  TOLOBE ( 1 , SIGMA .BMCPHP , BMSPHP . 12) 

6 1 0>  I F (BET . LE . THPT . AND . THPT . LE . P 1 2 -BET)  GO  TO  1 

620>  C 

B30>  C  THETA'  <  BETA  OR  THETA'  >  BO-BETA 

S40>  S IQsWLENMBMTAN ( THPT -BET)  <Ht2/ (9 . MP IMSTHP I 

B50>  IF ( LOBE .GT .0)  CALL  TOLOBE 1 2, $  I G, BMCPHP, BMSPHP, -L2) 

6B0»  S I QMA *  $  I  QUA  ♦  3 1 G 

670>  GO  TO  1 

6B0»  END 

690>  C 

700>  C 

705 >  C  <M» <«*><•(><«><<•> <MXMXMXMXM> <MXM> <MXMXM> <MXM> <MXM> 

71 0>  SUBROUTINE  CYLIN 

720>  C 

730>  C  WITH  LOBE 
740>  C 

750>  IMPLICIT  REALMS  (A-Z) 

760»  INTEQERM2  LOBE 

770»  COMMON/RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP . SPHP . 

790*  t  C2PHP.S2PHP, RAD, POLAR, LOBE, 

780>  S  A.B.L 

S00>  C 

«10>  IF (STHP . LE .0. I  RETURN 

S20>  A2=A»A 

830>  B2SBNB 

•  40>  Z= ( A2*CPHP*CPHP*B2»SPHPMSPHP)MMl . 5 

«50>  RCSC  =  2 .  MP I NLMLMA2MB2 / ( WLENMZ  I 

•  60>  C 

S70>  D*8 .MP ImCTHPmCTHPmZ 

SSO>  IFIABS(D) .LT. CUTOFF)  GO  TO  1 

SSO>  C 

900>  C  THETA’  f  90 

910*  SIQMAS2.MWLENNA2MB2MSTHP/D 

920>  RCSC^AMI HI (RCSC .SIGMA) 

930>  IF (LOBE . LE . 0)  RETURN 

940>  CALL  T0L0BEI1 , RCSC/2. , AMCPHP , BMSPHP , 0 . ) 

950>  CALL  TOLOBE (2, RCSC/2. .AMCPHP, BMSPHP, L) 

950>  RETURN 

970>  C 

9«0>  C  THETA’  *  90 

990>  1  IF(LOBE.QT.O)  CALL  TOLOBE (1 , RCSC, AMCPHP .BMSPHP , L/2 . I 

1000>  RETURN 

1010>  END 

1020>  C 
1030>  C 

1  03 5 >  C  <NXMXMXMXMXNXMXMXMXMXMXMXMXMXMXNXNXNXN> 

1040>  SUBROUTINE  CONI 

1 030>  C 

1080>  C  WITH  LOBE 
1070*  C 

10I0>  IMPLICIT  REALMS  (A-Zt 

1090>  INTEQERM2  LOBE 

1100*  COMMON/RCSCOM/CUTOFF ,WLEN,PI,PI2, RCSC , THP , CTHP , STHP , PHP . CPHP . SPHP , 

1110*  9  C2PHP.S2PHP, RAD, POLAR, LOBE, 

1 120>  9  ALP, NU. A, LI  ,L2 

11S0>  C 

1140*  IFICTHP.GE.O. .AND. SIN(ALPMRADI ,GE. STHP)  RETURN 

11S0>  NU2*NUmNU 

1100*  TALP«TAN(ALPMRAD) 
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170*  D*SQRT(S2PHP*NU2MC2PHP) 

1 00*  POLANQsATANI -NU/ (TALPMD) ) 

190>  RCSCs8.NPIM(L2mm1 . 3-L1WH1 . 8 )MN2MTALPNM4 

200>  RCSC=RCSC/ 1 8 .MWLENNNU2MABS (COS (POLANQ)nm3 ) ) 

21 0>  Z=(STHP«TALP»D*NU»CTHP)»*2 

220>  Z*ZN8 .MP IM5THPMD 

230>  IFIZ.QT. CUTOFF)  80  TO  1 

240>  IFILOBE.LE.O)  RETURN 

2SO>  ls(L1«L2)/2. 

280>  CALL  TOLOBE (1 , RCS , LMTALPMCPHP , LMNUMTALPMSPHP , L ) 

270>  RETURN 

280>  C 

290>  C  OFF  -  NORMAL  INCIDENCE 
300>  1  SIGMA? I ( STHP/D) * ICTHPMTALP/NU) ) MM2 

310>  BIOMASS  1 8MAMWLENMNU2NNUNTALP/Z 

S20>  RCSCsAMINI IRCSC,SI8MAMIL1*L2)I 

880>  IFILOBE.LE.O)  RETURN 

S40>  CALL  TOLOBE11 .AMIN1IRCSC/2. .LImSIQMA) .L1MTALPMCPHP, 

850>  1  LImNUNTALPMSPHP , L 1 1 

360>  CALL  TOLOBE(2,AMIN1 IRC8C/2. ,L2mSI8MA) , L2MTALPMCPHP , 

370>  1  L2MNUMTALPMSPHP.L2) 

880>  RETURN 

380>  END 

400>  C 
41 0>  C 

415*  C  «NXNXMXMXMXMXMXMXMXNXmxmxMXMXMXMXNXMXM> 

420>  SUBROUTINE  PARAB 

430>  C 

440>  C  WITH  LOBE 
4S0>  C 

4B0>  IMPLICIT  REALMS  (A-Z) 

470>  INTEQERM2  LOBE 

4S0>  COMMON/RCSCOM/CUTOFF .WLEN.PI , P 1 2 , RCBC , THP , CTHP , STHP , PHP . CPHP . SPHP 

490>  S  C2PHP.S2PHP, RAD, POLAR, LOBE, 

800>  S  P 

S10>  C 

S20>  C  RCS  IS  UNDEFINED  AS  THP  APPROACHES  PI 

S30>  IF  I AB  S I  THP ) .8T.70.MRAD)  RETURN 

S40>  RCSCS4.MPIMPMP/CTHPMM4 

SSO>  IFILOBE.LE.O)  RETURN 

S60>  TTHPsTANITHP) 

870>  CALL  TOLOBE I 1 , RCSC , 2 . MPMTTHPMCPHP , 2 . MPMTTHPmSPHP , -PmTTHPmTTHP ) 

5SO>  RETURN 

S90>  END 

SOO>  C 

S10>  C 

81 8>  C  «MXMXMXMXMXMXMXMXMXMXMXMXMXMXMXMXMXMXM> 

020>  SUBROUTINE  SOLID 

S30>  C 

B40>  C  WITH  LOBE 
SS0>  C 

BB0>  IMPLICIT  REALMS  (A-Z) 

•70>  INTE8ERM2  LOBE 

SS0>  COMMON/RCSCOM/CUTOFF , WLE  N , P I . P 1 2 , RCSC , THP , CTHP , STHP . PHP , CPHP . SPHP 

S90>  8  C2PHP.S2PHP, RAD, POLAR, LOBE, 

700*  8  A,B,C 

710>  C 

720>  RaAMBMC 

730>  D« (AMSTHPMCPHP)MM2 

740>  D*D* (BmSTHPmSPHP)MM2 
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DsD»(C»CTHf>)««2 
RCSC*P InRnR/ (DND) 

IF  I  LOBE .GT .0)  CALL  TOLOIE ( 1 , RCSC , ANSTHPNCPHP , BnSTHPnSPHP , 
1  CMCTHP) 

RETURN 

END 


<«><«><«> <« ><«><«> <«><«> <«><«> <«><«> <«><«> <«><«><«><•><*> 
SUBROUTINE  TORUS 


WITH  LOBE 


IMPLICIT  REALM  (A-Z) 

INTEQERN2  LOBE 

COMMON/RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP , CTHP . STHP , PHP , CPHP . SPHP . 
9  C2PHP.S2PHP, RAD, POLAR .LOBE, 

9  A ,  B 


RCSC=8 . NP INNSmBNANA/WLEN 

IFtSTHP.LT. CUTOFF. AND. LOBE. 8T. 01  CALL  TOLOBE 1 1 , RCSC , 0 . , 0 . , 0 . I 
IFISTHP.LT. CUTOFF)  RETURN 
SIQlap|n(BNA/STHP«B«B) 

3182*0. 

IF (LOBE .QT .0)  CALL  TOLOBE 1 1 ,AM INI (RCSC . SI81 1 .AnCPHP .BNSPHP ,0. 1 
IF (ABS(CTHP) . LE . 6/ (2 .NA) J  80  TO  1 
S I82*P In IBMA/STHP-BMBI 

I F 1  LOBE . QT . 0 1  CALL  TOLOBE 1 2 . AM I N1 1 RCSC . S 1 82 1 . -ANCPHP , -BNSPHP , 0 . 1 
RCSC*AMIN1 (RCSC, SIQ1 *3102) 

RETURN 

END 


4N><N><N><N><N><N><N><N><N><N><NXN>«N><N><N><N><N><N><N> 
SUBROUT  I NE  WEDQE ( COOR . CTH , STH , CPH . SPH I 


WITH  LOBE 


IMPLICIT  REALN4  (A-ZI 
INTE8ERN2  LOBE 

COMMON/RCSCOM/CUTOFF ,WLEN, P I , P 1 2 , RCSC , THP .CTHP . STHP , PHP .CPHP , SPHP , 
»  C2PHP.S2PHP, RAD, POLAR, LOBE. 

8  AN8LE.L.6 
REAL  COOR (9) 


WEOMAX*B/(WLENnWLEN) 

AN8*AN8LENRAD 

M*  IP lNLNL/14 .NANQNAN8) InSINIP INP 1/ (2 .NANO) INN2 
T*COS (P INP 1/12 .NANO)  I 
CaCOSIP InPHP/ANQ) *T 
D«1 , -T 

IFIABS(C). LE. CUTOFF. OR. ABBID) .LE. CUTOFF)  RETURN 


IF IPOLAR.EQ. 2 . )  80  TO  2 

RCSC*2 . N I -COOR 1 3 ) NSPH+COOR 1 6 ) NCPHINN2-1 . 

RCSC«MN|RCSC/C«1 ./DINN2 

80  TO  4 


RCSC82 . N (COOR 1 3 ) NCTHNCPH*COOR ( 6 1 NCTHNSPH-COOR 1 9 ) N 
6  STHINN2-1. 


'A  i'X/uV.  V-  •V'WV 
afNTif>!l  LsiLig  |CL4  ‘aVL&VW  A-S  -  •  V *>N>  •jlV 


fc V5  Wvl  s  >.'J 


OO  MO  M  o  ooo  ooo 
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RCSC*MMIRCSC/C+1 . /D)MM2 

4  DU«AIS I P 1 2 - ACOSK I -COON  1 2 ) MSTHMCPH-COOR 1 6 ) MSTHM 
t  SPH-COOR I9)NCTH) I 

RCSC«AMIN1 INC SC .WEDMAX) 

IF IDUMDU.LT. CUTOFF)  GO  TO  • 

SIGMAaRCSCMIVLEN/ IP INLMDU!  1MM2/4 . 

RCSCaAMINI (RCSC .SIGMA) 

I F (LOSE . LE .0)  RETURN 

CALL  TOLOSE 1 1 , RCSC/2 . , 0 . , 0 . , 0 . ) 

CALL  TOLOBE(2, RCSC/2. .0..0..L) 

RETURN 

C 

5  IF  I  LOBE . LE . 0)  RETURN 
CALL  TOLOBE 1 1 . RCSC , 0 . , 0 . , L/2 . ) 

RETURN 
END 


<M><M><M><N><M><M><M><N><N><N><M><MXN><M><M><N><M><N><M> 

SUBROUT  I NE  VINE ICOOR , CTH , STH . CPH , SPH) 

WITH  LOBE 

IMPLICIT  REALMS  IA-Z) 

INTEQERM2  lobe 

COMMON/RCSCOM/CUTOFF .WLEN.P I ,P 12 . RCSC , THP . CTHP , STHP , PHP , CPHP . SPHP 
S  C2PHP.S2PHP. RAD .POLAR. LOBE, 

S  L.A 

REAL  COORISI.Pm.ZKIS) 

IFISTHP.LE.O. )  RETURN 
GAMMA>1 .7S 

ZK 11 ) *COOR I S ) M I -CTH) -COOR I B ) M I -STHMSPH) 

ZKI2)  ■-  ICOORI7)MI-CTH)  -COORIDMl  -8THMCPH) ) 

ZK I S ) »COOR I 7 ) M I -STHMSPH ) -COOR I S ) M I -STHMCPH ) 

IFIPOLAR.EQ.2. )  GO  TO  2 

P(t)3-SPH 

PI2MCPH 

PIS)aO. 

GO  TO  S 
P(1)*CTHMCPH 
P(2)*CTHMSPH 
PISIM-STH 

SPS>ZK I 1 IMP  M ) *ZK(2)MP (2) *ZK (SIMP (3 ) 

SPSaSPS/SQRT ( ZK I I ) MZK (1 ) ♦ZK I 2 ) MZK I 2 ) *ZK I S ) MZK I S ) ) 

PS*ASINIAMAX1 1-1 . , AMIN1 1 1 . , SPS) ) ) 

RCSC*P INLMLMCOS (PS)MM4/ IP  1 2MP 1 2«ALOG (WLEN/ (GAMMAMP IMA) )MN2) 
IFtSTHP.LT. CUTOFF. AND. LOBE. GT. 01  CALL  TOLOBE I 1 , RCSC , 0 . . 0 . • 0 . ) 

I FISTHP.LT. CUTOFF)  RETURN 
SIGMAi|«LENMTANITHP)MCOS(PS)MM2)MM2 

S I GMA*S I QMA/ 1 1 • . MP I « I P I 2mP 1 2 *ALOG ( WLEN/ (GAMMAMP I MAMSTHP I )MM2 ) I 
RCSCiAM I N1 1 RCSC . 2 . MS  I QMA I 
IF (LOBE . LE .0)  RETURN 
CALL  TOLOBE ( 1 , RCSC/2 ..0..0..L/2.) 

CALL  TOLOBEI2, RCSC/2., 0. ,0. ,-L/2.) 

RETURN 

END 
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2I2S>  C  <•><»> <«><i»> <»><«*> <«>«»><»> <*><»> <»><»><*><«><»><»> <*><»> 

2S30>  SUBROUTINE  LWIRE (COOR ,CTH , STH ,CPH , SPHI 

2B40>  C 

2IS0>  C  WITH  LOBE 
2B(0>  C 

2B70>  IMPLICIT  REALM  IA-ZI 

2BB0»  INTEGERS  LOBE 

2BB0>  COMMON/RCBCOM/CUTOFF .WLEH.PI , P 1 2 , RCSC , THP , CTKP . STHP , PHP . CPHP , 8PHP , 

2900>  $  C2PHP.S2PHP, RAD, POLAR. LOBE, 

2t!0>  B  A 

2B20>  REAL  COOR (B I 

2B30>  C 

2940>  I F (POLAR . EQ . 2 . )  BO  TO  2 

2950>  QAMMA*ACOSK ( -COOR ( 2 1 «SPH*COOR ( 3 1 «CPH I 

2BB0>  00  TO  3 

2B70>  C 

2990*  2  QAIMA*AC08K ( COOR 1 2 ) NCTHNCPH+COOR 1 9 ) NCTHNSPH-COOR 1 8 ) »STH I 

2990>  C 

3000>  3  KA2* (4.*P INA/VLEN1MSTHP 

3010>  C02*C0SI0AMMA)««2 

9020>  80208 1 NIQAMMA  !  NN2NCTHPNCTHP 

3030>  CALL  BESJ (KA2 , 0, JO. . 1 , IER) 

3040>  IF ( IER.NE .0)  00  TO  9 

3080>  CALL  BESJ (KA2 , 2 , J2 , . 1 . IER) 

3030>  IF! IER. ME. 01  00  TO  9 

3070>  RCSC: | SQ2C«CQ2 ) »J0* ( S02C-CQ2 ) N J  2 

3080>  RCSC*P I MANAnRCSCMRCSC 

3090>  IF (LOBE . LE . 0)  RETURN 

3tOO>  I F ( STHP . LE . CUTOFF )  CALL  T0L0BEI1 .RC5C.0. ,0. ,0. ) 

31 1 0>  IF (STHP. OT. CUTOFF)  CALL  TOLOBE ( 1 , RCSC , AttCPHP , AN8PHP , 0 . I 

3120>  RETURN 

3130>  B  RCS*0. 

3140»  RETURN 

3150>  END 

31B0>  C 
91 70>  C 

31 79>  C  «■> <«><«> <«><«><*> <*><«><«> <«><*> <«><*> <«><«> <»><«> <«><«> 

3 1 80>  SUBROUTINE  CAVA 

3190>  C 

3200>  C  WITH  LOBE 
3210>  C 

3220>  IMPLICIT  REALN4  (A-ZI 

3230>  INTEQERN2  LOBE 

3240>  COMMON/RCSCOM/CUTOFF ,f LEN . P I , P 12 . RCSC , THP , CTHP , STHP , PHP , CPHP , SPHP , 

32S0>  9  C2PHP.S2PHP, RAD, POLAR, LOBE, 

3290>  9  A 

3270>  C 

3280>  IF(ABS(THPI .QT.70.NRAD)  RETURN 

3290>  T*(4.MPI*A/(WLENnWLEN) )N*1 .5 

S300>  RCSC* . 4MTNVLENNWLEN 

3310>  TTHPsTANITHP) 

3320>  IFITTHP.LT. CUTOFF. AND. LOBE. QT.O)  CALL  T0L0BEI1 .RCSC.O. .0. ,0. ) 

3330>  IFITTHP. LT. CUTOFF  I  RETURN 

3940>  S I OMA* . BttP I NA/TTHP 

3350»  RCSC*AMIN1 (RCSC, SI OMA I 

33B0>  IF (LOBE . LE . 0)  RETURN 

3370>  D*SQRT (4 .MA/P 1 1 

3390>  DS2*D«STHP/2. 

3390>  CALL  T0L0BEI1 ,RCSC,DS2NSTHPttCPHP,D32NSTHP«SPHP,DS2«CTHP) 

3400>  RETURN 
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<»><«> <»><*><«><*><«> <«> <»X«X«X«> <M>  <«>  <«><»> <«X«XN> 

SUBROUTINE  CAVD 
WITH  LOBE 

IMPLICIT  REALM  (A-Z) 

INTEGERS  LOBE 

COMMON/RCSCOM/CUTOFF , WLEN , P I ,P 1 2 , RCSC , THP ,CTHP , STHP , PHP .CPHP , BPHP , 
S  C2PHP.S2PHP, RAD, POLAR, LOBE, 

$  D 

IFIABS(THP) .GT.70.dRAD)  RETURN 
RCSC*  .05d(2.dPI  dD/WLEN )  ddSdWLENdWLEN 
TTHPsTANITHP) 

IFCTTHP.LT. CUTOFF. AND. LOBE. QT.O)  CALL  TOLOBEI1.RCSC.O..O..O.) 
IFCTTHP.LT. CUTOFF)  RETURN 
SIGMA& . 05d(2 .dP ldD)dd2/TTHP 
RCSCsAM INI! RCSC , S I QMA ) 

IF  1  LOBE . LE . 0)  RETURN 
DS2>DdSTHP/2. 

CALL  TOLOBE  ( 1 . RCSC , DS2dSTHPdCPHP , DS2MTHPMPHP , DS2NCTHP ) 

RETURN 

END 

<*XdXdXdXdXdXdXdXdXdXdXdXdXdXdXdXdXdXd> 

SUBROUTINE  CFLAT 
WITH  LOBE 

IMPLICIT  REALM  (A-Z) 

COMMON/RCSCOM/CUTOFF , WLEN , P I , P I 2 , RCSC , THP , CTHP , STHP , PHP . CPHP , SPHP . 
S  C2PHP.S2PHP, RAD .POLAR .LOBE, 
t  A 

IF  CABS (THP) . 6E . P 1 2)  RETURN 
RCSCM.dP  ld(P  IdAdR)  dd2/  (WLENdWLEN) 

D*S .dP  I dSTHPdTAN  I  THP )  ««2 

IFID.LT. CUTOFF. AND. LOBE. QT.O)  CALL  TOLOBE (1 .RCSC. 0. .0. ,0. ) 

IF  ID. LT .CUTOFF )  RETURN 

SIQMA*AdWLEN/D 

RCSC  s  AM  INK  RCSC ,  2 .  «S  I  QMA ) 

IF (LOBE . LE .0)  RETURN 

CALL  TOLOBE ( 1 , RC8C/2 . , AdCPHP , AdSPHP , 0 . 1 
CALL  TOLOBE  12  .RCSC/2  . ,  -AdCPHP  ,  -AMPHP  ,0 . ) 

RETURN 

END 

<*XdXdXdXdXdXdXdXdXdXdXdXdXdXdXdXdXdXd> 

SUBROUTINE  RFLATI IPASS) 

WITH  LOBE 

IMPLICIT  REALM  (A-Z) 

INTEQERdZ  LOBE, IPASS 


/vyviviv.vy.v^.i .v .y  .yo 
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3890  COMMON/RCSCOM/CUTOFF , WLEN.PI , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP . SPHP , 

4000>  9  C2PHP.S2PHP, RAD, POLAR, LOSE, 

401 0>  9  A , B 

4020>  C 

4030>  IF IABSITHP) ,0E , P 12)  RETURN 

4040>  RCSC  =  P I M 1 8 . «A*B/WLEN)*tt2 

4050>  S2THP* STHP* STHP 

4060>  S4THP*S2THP*S2THP 

40/0>  I F I S4THP . LT . CUTOFF )  QO  TO  4 

40S0>  IFIS2PHP.LT. CUTOFF)  GO  TO  2 

40#0>  IFIC2PHP.LT. CUTOFF)  00  TO  3 

4100>  SIQMA-VLENNWLEN/  I64.NP  IHM3*S4THiP*S2PHP*C2PHP) 

4 1 1 0  >  S IQMA a 2 . NS  I  QUA 

4120>1  RCSC* AM INI I  RCSC,  2.  IIS  IQMA) 

41 30>  IF (LOBE . LE . 0)  RETURN 

4140  C 

41 50*  RCSC4*RCSC/4 . 

4180  IF ( IPASS.NE  .2)  CALL  TOLOBEH  , RCSC4 , A , B , 0 . ) 

4170>  IF  I IPASS.NE .2)  CALL  TOLOBE ( 2 , RCSC4 , -A , B , 0 . ) 

4180>  IF  I IPASS.EQ.2)  CALL  TOLOBE M .RCSC4 , A , -B,0 . ) 

4180  IF! IPASS.EQ.2)  CALL  TOLOBE (2 .RCSC4 , -A, -B,0. I 

4200>  RETURN 

4210>  C 

4220  C  PHI’  *  0  OR  180 

4230  2  IFI  IPASS.EQ.2)  RETURN 

4240>  SIQMA*B«B/I4.MP1HS2THP) 

4250>  IF (LOBE . LE . 0)  00  TO  1 

4260  SIQSAMINKRCSC/2.  .SIQMAI 

4270*  CALL  TOLOBE 11 ,SIQ. A, 0. ,0. ) 

4280>  CALL  TOLOBE 12 , S 1 0 , -A , 0. , 0 . ) 

4280  S  RCSCsAMI  N1 I RCSC , 2 . *S  IQMA) 

4300>  RETURN 

4310>  C 

4320  C  PHI  *  a  90  OR  -80 

4330  3  IFI  IPASS.EQ.2)  RETURN 

4340>  S IQMA = A* A/ (4 .*P IHS2THP ) 

4380  IF  I  LOBE .  LE  .  0)  00  TO  1 

4380>  SI0*AMIN1 (RCSC/2 . , SIQMAI 

4370»  CALL  TOLOBE II , SIQ.O . ,B, 0 . ) 

43S0>  CALL  TOLOBE I2,SIS,0.,-B,0.) 

4380  00  TO  S 

4400>  C 

441 0>  C  THETA’  a  0 

4420>  4  IFI IPASS.EQ.2)  RETURN 

4430  IF  I  LOBE . LE . 0)  RETURN 

4440>  CALL  TOLOBE II ,RCSC , 0. ,0. ,0. I 

44SO  RETURN 

4400  END 

4470>  C 

4480  C 

44S8>  C  <*> «»><»> <«><«> <*><*> <«><«><«> <*><*> <«><*><«> <*><«> <«><*> 

4480  SUBROUTINE  TWAVE 

4482>  C 

4484>  C  WITH  LOBE 
4S00>  C 

4810  IMPLICIT  REALM  IA-Z) 

4S20»  INTE0ERM2  LOBE 

4830  COMMON/RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP , SPHP , 

4840  S  C2PHP.S2PHP, RAD, POLAR, LOBE, 

4880  8  P ,  L ,  QAMMA ,  A 
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oo>  c 

70  IF ( POLAR. EQ. 1 . )  RETURN 

«0>  KL*2 .MR ImL/WLEN 

«0>  I F (KL . LE .4. )  RETURN 

OO  Q*2.35ttAL0810(2.«KL)*.5S 

10>  T«1.-P«CTHP 

20>  IFIT.LE. CUTOFF)  RETURN 

30>  F*STHP*SI N (KL*T/ 12 .NP) 1 /T 

40  RCSC* (GAMMA*F*F«WLEN/Q)Mft2/P  I 

45>  IF(LOBE.QT.O)  CALL  TOLOBEM , RCSC , AftCPHP , ANSPHP , -L/2 . ) 

30>  RETURN 

00>  END 

70>  C 
I0>  C 

iB5>  C  «»>  <»>«*><*><»><»><*>  <»><*><)»> <»><«><»><«><«><«> <*><»><*> 

»0>  SUBROUTINE  DIHED 

OO  C 

10>  IMPLICIT  REALM  IA-Z) 

20  INTEQERH2  LOBE 

30>  COMMON / RCSCOM/CUTOF  F , WLEN  p  P I . P 1 2 , RCSC , THP . CTHP . 8THP . PHP . CPHP . SPHP . 

40>  S  C2PHP.S2PHP, RAD .POLAR, LOBE, 

SO>  S  A,B,C 

BO>  C 

70>  IFIPHP.LE.O. .OR.PHP.GE.PI2)  RETURN 

’75  >  PI4SPI/4. 

BO>  QAMMA=ATAN2IB,A) 

90>  DU*P I 2-THP 

95>  C  COMPUTE  NORMAL  COMPONENT 

IOO>  I F ( PHP . 8T . GAMMA I  80  TO  1 

II  0>  RCSCM .  «P  I N I AOB/VLEN I  «N2«4 .  NS2PHP 

I20>  80  TO  2 

I30>  1  RCSCs4.*PHllCM/WLEN)M2«4,»C2PHP 

40>  2  I FIDUM2.8T. CUTOFF)  80  TO  4 

41 >  IF(LOBE.LE.O)  RETURN 

IS0>  4  IFISTHP.LT. 0.5)  80  TO  3 

I80>  S I GMAsVLENtfWLE  NNRC5C/ 1 8 . * ( P I MB*DU )  <H>2 1 

1 70>  RCSC* AM INI (RCSC, SIGMA) 

IB0>  RETURN 

IB0>  3  RCSC*0. 

»00>  RETURN 

H 0>  END 

I20>  C 
130  C 

I35>  C  <N> <«><»> <«><n> <«><«><«> <»><«> (»><li><l><» <«> <«><«> <«><«> 

I40>  SUBROUTINE  TCORN 

ISO  C 

ISO>  C  WITH  LOBE 
170  C 

ISO>  IMPLICIT  REALM  IA-Z) 

IBO>  INTE8ERN2  LOBE 

I00>  COMMON/ RCSCOM/CUTOF F,«LEN,PI,PI2, RCSC , THP , CTHP , STHP . PHP , CPHP . SPHP , 

HO  I  C2PHP.S2PHP, RAD, POLAR, LOBE, 

>20  S  A.B.C 

130  C 

•40  I F I  THP .  LE .  0 . 0 .  OR .  THP .  GE .  P 1 2 )  RETURN 

ISO*  IFIPHP.LE, 0.0. 0R.PHP.8E. PIS)  RETURN 

ISO  CALL  ORDER(STHPNCPHP,STHP«SPHP,CTHP,L,M,N) 

•70  CALL  ORDER IA, B, C , T1  ,T2,T3) 

ISO  LAU/T3 

ISO*  MS«M/T2 
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S100>  NC*N/T1 

511 0>  I F ( LA+MB .GT . NC)  QO  TO  1 

5120*  AREAM  .*L*M*T1 / I LA+MB+NC) 

51 30>  QO  TO  2 

51 40>  1  AREAs|lA»LA*MB»MB*NC«NC)/<LA*MB*NC) 

51 50>  AREA>AMB»C» I LA*MB«NC-2 . mAREA ) 

51 60>  2  RCSC s 4 . HP  I  MARE AM ARE A/ (VLENoWLEN) 

5170>  IF (LOBE .LE .0)  RETURN 

51B0>  D«1  ./SQRTM  ./AM2*1 ./BM2  +  1 . /CM2) 

51»0>  CALL  TOLOBE 1 1 , RCSC , DNSTHPHCPHP , DttSTHPttSPHP , DttCTHP) 

S200>  RETURN 

S210>  END 

S220>  C 
S250>  C 

5235>  C  <»> <»><«> <»><«>  <<»><«> <«><«><«><*><*><*> <*><«> <»><»> <*><»> 

5240>  SUBROUTINE  RCORN 

5250>  C 

S2S0>  C  WITH  LOBE 
5270>  C 

52S0>  IMPLICIT  REALM  (A-ZI 

5290>  INTEQERM2  LOBE 

S300>  COMMON/RCSCOM/CUTOFF ,WLEN , P I , P 1 2 . RCSC , THP , CTHP . STHP . PHP , CPHP . SPHP . 

531 0>  S  C2PHP.S2PHP, RAD .POLAR, LOBE, 

S320>  5  A , B ,C 

5330>  C 

5340>  I F I THP . LE . 0 . 0 . OR . THP . GE . P 1 2 1  RETURN 

5350>  I F I PHP . LE . 0 . 0 . OR . PHP . QE . P 1 2 1  RETURN 

S360>  CALL  ORDER  I STHPMCPHP , STHPMSPHP , CTHP , L ,M, N) 

5370>  CALL  ORDER (A , B ,C . T1 , T2 , T3 ) 

5SS0>  LASL/T3 

53»0>  MBSM/T2 

5400>  NC=N/T1 

541 0>  I F (MB . LT.NC/2 . I  QO  TO  1 

5420>  AREA*  L*T2* ( 4 . «T1 -N/MB I 

S430>  QO  TO  2 

S440>  1  AREAM. NL«M*T1«T1/N 

S4S0>  2  RCSCM.«PI«tAREA/VLENIM2 

5460>  I F I  LOBE . LE . 0)  RETURN 

S470>  D>1  ,/SQRTM  ./A««2*1  ,/B«N2«1  ./CM2) 

54S0>  CALL  TOLOBE I 1 , RCSC , DttSTHPttCPHP , DNSTHPHSPHP , DMCTHP I 

S4S0>  RETURN 

SSOO>  END 

5510>  C 

5520>  C 

5525>  C  mxnxmxnxnxmxm* <«><«> <«x«x«xnx»x«xi» <«><«><«> 

SS30>  SUBROUTINE  ECORN 

SS40>  C 

5550>  C  WITH  LOBE 
S5S0>  C 

5570>  IMPLICIT  REALM  (A-Z) 

SSS0>  INTEQERM2  LOBE 

S5tO>  COMMON/RCSCOM/CUTOF F , WLE N , P I , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP , SPHP , 

5B00>  S  C2PHP.S2PHP, RAD, POLAR, LOBE, 

5510*  5  A.B.C 

5020*  C 

5030*  I F ( THP . LE .0.0 .OR. THP ,QE . P 1 2 )  RETURN 

5040*  IF (PHP . LE . 0. 0 .OR .PHP .QE ,P 12)  RETURN 

SOSO*  CALL  ORDER (STHPNCPHP , STHPOSPHP ,CTHP , L ,M,N) 

SOSO*  CALL  ORDER(A,B,C,T1,T2,T3) 

5070*  LA2>(L/T3)M2 


1 

* 


a 


N 


52 
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5810 

5890 

5700 

5710 

5720 

5750 

5740 

5750 

5780 

5770 

5790 

5790 

5900 

5910 

5920 

5950 

5940 


MB2s (M/T2INN2 
NC2s (N/T1 I»»2 
i F INC2 . GE . LA2+MB2 )  00  TO  1 

AREAS  I N/T1 )  NATAN  HI  LA2  «M92  HH>2 'NC2MNC2 )  / 1 4 .  MHC2«  I L/TS  HMM/T2  m 
AREAS  ARE  A*  (II/T2 I  NATAN  (H  LA2  *NC2 )  NN2 -MB2NMB2 1  / 

8  (4.NMB2NIL/TS)N(N/T1) )) 

ARE A s ARE A + ( L/T3 ) NATAN  I ( IMB2  +  NC2 INN2-LA2NLA2 ) / 

8  (4.NLA2N(M/T2)N(N/Tim 
QO  TO  2 

AREAS (M/T2 1  NATAN ( 2 . N ( L/T5 ) N ( N/T1 ) / IMB2+NC2 - LA2 ) ) 

AREAsAREA* IL/T3) NATAN  1 2 . N (M/T2 InIN/TI I/I LA2-IIB2  +  NC2 ) ) 
AREAsANBnCnAREA 
RCSC>4 . nP In ( AREA/WLEN) NN2 
IF (LOBE . LE . 0)  RETURN 
Dsl ,/SQRTM ./ANN2+1 ./BNN2+1 ./CNN2) 

CALL  TOLOBEM , RCSC , DnSTHPnCPHP .ONSTHPNSPHP .DNCTHP) 

RETURN 
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»>> 


i 


c 

c 

c 


c 

c 


FILE 
10>  C 
20> 

30* 

40* 

50* 

SI* 

S2> 

S3* 

54* 

55* 

60* 

70* 

75* 

30* 

90* 
100* 
110* 
120* 
130* 
170* 
110* 
190* 
220* 
222*  . 
224*  Cl 
220*  C 
230* 
240* 
290* 
260* 
270* 
210* 
290* 
300*  C 
310* 
320* 
330* 
340* 
3S0*  C 
360*  1 
370*  C 
390* 
390* 
400* 
410* 
420* 
430* 
490* 
470* 
490* 
492* 
494*  9 
490* 
900*  C 
810*  3 
820* 
830* 
840* 
880*  C 


■  8Y0:  HU  ,  121RCSSUB .FTN 
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4 

C 

8 


THIS  ROUTINE  ASSUMES  MATRICES  ARE  ROW  STRUCTURED 

SUBROUTINE  TOLOBEI IC.RCSL ,X,Y,Z) 

INTEGERS  DEBUG, COMB  IN, QL I  NT  „ 

LOGICALN2  GOULD, LP .OPENE , LPF .OLDQRF .MGRAPH.GDAUTO 
COMMON/SWITCH/  IOW.THETAI4) .GOULD, LP .OPENE , LPF , OLDQRF .MGRAPH, 

1  QDAUTO, DEBUG, COMBIN, GLINT, IBUG.LBUG. ISTATS, IZETA 

REALM  SLOBE  14,191 ,21 

REALM  ORQN f 3) . COMP  1 9 ) , MTRAN 13,3), TMP 1 3 ) 

COMMON /RC SLOB/ I OL , ORQN , COMP , SLOBE . BETA , 

9  THR . CTHR , STHR , PHR , CPHR , SPHR , KLOBE 

COMMON/RCSCOM/CUTOFF , WLEN , P I , P 1 2 , RCSC , THP , CTHP , STHP , PHP , CPHP , SPHP , 

9  C2PHP.S2PHP, RAD, POLAR .LOBE, 

9  ARR 

REALM  RCSI191)  .RCS2M91) 

COMMON  /RPRCS/  I  RPR, IRCS, RAM.RCS, RCS2 

IF (LBUG.NE.O)  WRITEI2.7)  IC.RCSL ,X,Y,Z 
FORMAT  I’  TOLOBEs ' , II , 4E1 1 .  3 ) 

TMP(1J»X 

TMP(2)sY 

TMPI3XZ 

CALL  MATVEC (COMP, TMP, TMP) 

TMP 1 1 1 >TMP 1 1 ) ♦ORQN 1 1 ) 

TMP(2)*7MP  (2)  *0RQM2) 

TMP « 3 )  ■  TMP 1 3 ) ♦ORGN 1 3 ) 

CALL  MATPUT IMTRAN , CTHARCPHA  .CTHRMPHR , -STHR . 

1  SPHR, -CPHR, 0. , 

2  -STHRtlCPHR ,  -STHRMPHR ,  -CTHR) 

IF (BETA. LE. 8.)  GO  TO  1 

CALL  MATVEC IMTRAN, TMP, TMP) 

I F 1 1 LBUQ . AND .2) . EQ . 0 )  GO  TO  9 

WRITE  12, 4)  IC,THR,PHR,X,Y,Z, ORGN. MTRAN 

FORMAT I1X, 1 2 , 9E1 0 . 3/3E1 1 .3/3E11 .3/3E11 .3) 

IF  I IRPR.EQ.O)  GO  TO  3 
ARQM .  oP  I NTMP 1 3 )  /WLEN 
SQiSQRT (RCSLtfRAM) 

RCS II RCS ) *RCS 1 1 RCS ) ♦ SQnCOS I ARG ) 

RCS2  I  IRCS)  *RCS2 1 1 RCS) *SO*S IN(ARQ)  . .  . . . 

IF  I (LBUQ. AND. 1 ) .NE .0) WRITE  12 ,9 1 IRCS.RCSL .RAM.TMP ,RCSI IRCS) , RCS2I IRCS) 
FORMAT  I ’  RCSL* ’ , 1 3 ,7E1 1.3) 

IF (KLOBE . LE .0)  RETURN 

SLOBEM  ,  IRCS,  Id iRCSLNRAM 
SL0BEI2,  IRCS,  ICUTMPM) 

SLOBEM,  IRCS,  1C) *TMP (2 1 
9LOBEI4,  IRCS,  ICHTMPI3) 
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PV 


c; 


A 


3 
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i 
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N 


K 

I 


s 


S60> 
s7o>  2 
50Q>  c 
S90> 
0OO> 

01O>  c 

02O>  C 
830>  C 
04O> 
650> 
00O> 
870> 
880> 
09O> 
700> 
710> 
720*  C 
730>  C 
740>  C 
750>  C 
700> 
770> 
70O> 
790> 
800> 
«10> 
020> 
930> 
840>  3 
330>  2 
00O>  1 
970* 
980>  4 
09O> 
900> 

91 0>  C 
920>  C 
930>  C 
940>  C 
9S0> 
90O> 
970> 
930> 
990» 
1000> 
1010>  3 
1020>  1 
I030> 
1040>  4 
1090> 
1O0O> 
1070>  C 
1090>  C 
1090> 
1100> 
1110> 

1 1 20> 
1130> 
1140> 
11S0> 


IF ( ILBUQ. AND.  1 ) .NE .0)  WRITE  12, 2)  1C,  IRCS,  ( S LOBE  I  I ,  IRCS,  10,1=1,4) 
FORMAT ( '  TOLOBE*  * , 212 , El  1 . 3 , 3X , 3F0 .3) 

RETURN 

END 


SUBROUTINE  ORDER IA1 ,A2,A3,B1 .B2.B3I 
T1=AMIN1 (At ,A2,A3I 
T3  *AMAX 1 1 A  t , A2 , A3 ) 

B2»A1*A2*A3-T1-T3 

B1=T1 

B3  =  T3 

RETURN 

END 


THIS  ROUTINE  ASSUMES  MATRICES  ARE  ROW  STRUCTURED 

SUBROUTINE  MATMULIA.B.C) 

REALM  AI9)  ,B(9)  ,T(9)  ,C(9) 

DO  1  Js1,3 
ji3=i  j-i  m3 
DO  2  K=1 ,3 
JK=J13*K 
TIK*J13I«0. 

DO  3  Lai, 3 

TIK*J13)aTIK*J13JtA(J13*LU»BIK*IL-1J»3) 

CONTINUE 

CONTINUE 

DO  4  J*1,9 

CIJ)sT(J) 

RETURN 

END 


THIS  ROUTINE  ASSUMES  MATRICES  ARE  ROW  STRUCTURED 

SUBROUTINE  MATVEC (A , B ,C) 

REALM  AI9) , BIS), CIS) , TIM) 

DO  1  Jal.S 
J13s ( J-1 )n3 
TU)aO. 

DO  3  Lai ,3 

TU)aTU)+A(J13*L)«B(L) 

CONTINUE 
DO  4  J*1 ,3 
CIJ)*TU) 

RETURN 

END 


SUBROUTINE  MATPUT1A ,A1 1 , A1 2 , A1 3 , A21 ,A22,A23,A31 .A32.A33) 
REALM  A (9) 

All )>A1 1 
AI2XA12 
A 1 3 )* A 1 3 
A(4)*A21 
AI5XA22 


1  .'."T  r  .  p.  r:  r.'1'.  r:  *7  r- r -  r-  *v  r  r; »~.  »■.  n 
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I 

*  I 

-£  , 


.  S  , 


l< 


i> 

y- 


180> 

AI6)*A23 

170> 

AI7)*A31 

130> 

A(3)*A32 

180> 

A(8)*A33 

200> 

RETURN 

210> 

END 

220> 

C 

230> 

c 

240> 

SUBROUTINE  TRANSPIC.E) 

280> 

REAL  C (91 ,0(91 .E (9t 

280> 

Dill *C (11 

270> 

DI2)*d4> 

2I0> 

DI3IXCI7I 

1 2  9  0  > 

D  141*021 

9  00  > 

D(SlsCIS) 

1 9 1 0  > 

DlllsGIII 

U20> 

0(7)  *03) 

1 9  9  0  > 

DI8)*08) 

I940> 

0(8)  *08) 

1 9  5  0> 

DO  1  J*1,3 

1 3B0> 

1 

E(J)*DIJ) 

1 370» 

RETURN 

3«0> 

END 

I300> 

c 

I400> 

c 

1 4 1 0  > 

c 

I420> 

SUBROUTINE  BESJ IX , N, 8J , D, IER) 

1 43  0> 

c 

I440> 

c 

X*AR8UHENT  OF  THE  J  BESSEL  FUNCTION  DESIRED 

1 430> 

c 

NsORDER  OF  THE  J  BESSEL  FUNCTION 

I460> 

c 

BJsRESULTANT  J  BESSEL  FUNCTION 

I470> 

c 

DsREQUIRED  ACCURACY 

I480> 

c 

lERsRESULTANT  ERROR  CODE 

490> 

c 

0*NO  ERROR 

I800> 

c 

1*N  IS  NEOATIVE 

I810> 

c 

2*X  IS  NEOATIVE 

520» 

c 

3*REQUIRED  ACCURACY  NOT  OBTAINED 

1830* 

c 

4*RAN0E  OF  N  COMPARED  TO  X  NOT  CORRECT 

S40> 

c 

I880> 

BJ*0. 

I860> 

IF  IN) 10 , 20,20 

IS70> 

10 

IER«1 

I8I0> 

RETURN 

880> 

20 

IF (X)33 . 30,31 

I800> 

•0 

IF IN.8T.0)  RETURN 

161 0> 

BJ*1 . 

1620* 

RETURN 

ieio> 

38 

IER*2 

I840> 

RETURN 

I880> 

31 

IF IX*18. 132,32,34 

ie«o> 

31 

NTE8T*20.*10.«X-X»»2/3 

I870> 

00  TO  30 

I880> 

34 

NTEST*80.«X/2. 

I8I0> 

30 

IF(N-NTE8T)40,3B ,31 

1 7  00> 

38 

IER*4 

mo> 

RETURN 

1 7310* 

40 

IER*0 

I720> 

N1*N«1 

I740» 

BPREV*0. 

1790* 

C 

l  \ 
V 


i 


I 

J 


ll 


-5 


mV 


V"T 


V?J 


% 


v: 


3 


A 

V. 
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treo> 

C 

COMPUTE  STARTING  VALUE  OF  M 

1770> 

C 

1790> 

IFIX-9. 190,90,60 

1790> 

90 

MA*X*6. 

i«00> 

GO  TO  70 

1910> 

90 

MAsI .4»X*60./X 

1«20> 

70 

MB=N+ IF  IX IX) /4*2 

1 930> 

MZERO*MAXO IMA ,MG) 

t!40> 

C 

1«90> 

C 

SET  UPPER  LIMIT  OF  M 

1900> 

C 

1970> 

MMAX*NTEST 

1M0> 

100 

DO  190  M*MZERO,MMAX,3 

1990> 

C 

1900> 

C 

SET  FIM).F(M-I) 

1 91 0> 

C 

1920> 

FM1 *1 . E - 1 5 

1 930> 

FM*0. 

1940> 

ALPHA *0. 

19S0> 

IF (M- IM/2IW2) 120. 1 10,120 

1 990> 

110 

JTi-1 

1970> 

GO  TO  130 

1990> 

120 

JTal 

1990> 

130 

M2«M-2 

2000> 

DO  160  K31.M2 

2010> 

MKsM-K 

2020> 

RMK32 . NFLOAT IMK ) NFM1 /X-FM 

2030> 

FM*FM1 

2040> 

FM13BMK 

2090> 

IF (MK-N-1 )1 90 , 140 , 1 90 

2090> 

140 

BJsBMK 

2070> 

190 

JT3-JT 

2090» 

S*1 *JT 

2090» 

190 

ALPHA3ALPHA«BMKNS 

2100> 

BMK32.PFM1/X-FM 

21 1 0> 

IF (Ml ISO , 170, ISO 

2120> 

170 

RJsBMK 

2130> 

190 

ALPHA* ALPHA+BMK 

2140> 

B4 sB J /ALPHA 

21 90> 

1 F 1  ADS (B J -BPREVI -ADS ( D«B4 11200,200,190 

2190> 

190 

BPREVsDJ 

2170> 

IER33 

2190> 

200 

RETURN 

2190> 

END 

2200> 

C 

2210> 

C 

2220> 

SUBROUTINE  TRIMSOINB.B) 

2230> 

LOG  1  CALM 1  8(901 . BL 

2240> 

DATA  BL/'  '/ 

2290> 

NB*90 

2290> 

DO  1  Jsl.90 

2270> 

IFIBINB) .NE.BLI  RETURN 

2290> 

1 

NBsNB-1 

2290> 

RETURN 

2300> 

END 

2310» 

C 

2320> 

c 

2330> 

FUNCTION  ACOSKIAI 

2340> 

AC0SK*AC09(AMAX1 C-1 . , AM INI (1 . , A) ) l 

2390> 

RETURN 

2390> 

END 

>3 


'.-'A 
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.-1 

10>  CMMI 

20>  C 

H  H 

000 

ssss 

BBBB 

DDDD 

30>  C 

H  H 

0  0 

s 

B  B 

D  D 

40>  C 

HHHHH 

0  0 

sss 

BBBB 

D  D 

s' 

90>  C 

H  H 

0  0 

s 

B  B 

D  D 

60>  C 

H  H 

000 

ssss 

BBBB 

DDDD 

7 

90 

90 

100 

110 

120 

130 

140 

ISO 

160 

170 

180 

190 

200 

210 

220 

230 

240 

290 

260 

270 

280 

290 

300 

310 

320 

322 

324 

330 

340 

390 

360 

370 

390 

390 

400 

410 

412 

414 

416 

416 

420 

430 

440 

490 

460 

470 

490 

490 

900 

910 

920 

930 


C 

C 

C 


11 


4 

12 


C 

20 

C 


SUBROUTINE  H0SBDM0UT ,  LOOP ,  I0E ,  I  OH  .SLOBE  ,BUF ,  NB, 

1  FREQ, POL. PERSH, PITCH, ROLL, YAW. TILT.PHR) 

REAL*8  T1 , T2 , T3 , T4, 81 .92 , S3 ,84. Til, 822, VI .U2.U3.U4, 01 ,02 
REALMS  B ,N , P , ACSMAX , RC8MI N , TS , T6 
REALM4  SLOBE (6, 181 1 .SAVE  130) 

BYTE  BUFI1) 


A 

"*j 


INTEQERM2  DEBUG,  COMB  IN,  BUNT 

LOQICALM2  GOULD , LP , OPENE , LPF , OLDGRF .MGRAPH , GDAUTO 

COMMON /8WI TCH/  I OW . THETA  1 4 1 , GOULD , LP , OPENE . LPF , OLDGRF .MGRAPH . 

1  GDAUTO , DEBUG , COMB  I N , QL I  NT , I  BUG . LBUQ . I  STATS . 1 2ETA , STATS . NRCS . I OR 


.N 

_N 


LPsLOOP 

IF ( IOUT.EQ.O)  GO  TO  20 


COMPUTE  REQUIRED  ANGLES  FOR  GRAPH 


SLOPE*FLOAT (NRCS-1 I /AMAX1 1.001, THETA ( 2 ) -THETA ( 1 ) ) 
CONSTsI . -SLOPEMTHETA 1 1 ) 

ANQaSTATS 

NS«0 

IF (ANQ.QE .THETA (1 ) )  GO  TO  12 

ANG*ANQ*STATS 

GO  TO  11 

ANQ*ANQ»STATS 

NSaNS-1 

IF  IANQ.GT.THETAI2) I  GO  TO  13 
LP* I F I X I  SLOPEMANQ*CONST* . 001 ) 

NS*NS«1 


S 


READ  I IOE’ LP)  T1,T2,T3,T4,TS,T6 


S1*T1 

T12«T1MT1 

S2*IT12-T2)/2.D0 


3 


IFIS1MS2.EQ.0. DO. AND. IOUT.EQ.O)  RETURN 
IFIS1MS2.EQ.0.D0)  GO  TO  4 


Ci 


S3* (T12MT1 -3 .D0MT1MT2*2 . D0mT3)/6 .00 

84* IT1 2MT1 2-6 . D0MT1 2MT2+3 . D0mT2mT2*6 . D0MT1MT3-6 . D0«T4) 724 . DO 


V1*S1 

U2*2 . D0MS2 

U3*12.D0NS3 

S22*S2m82 

U4*8 . D0MS22+12 . 00MS1 MS3«1 32 . D0MS4 


Q1 *6  .  D0MS3/ ( SI MS2 ) 

Q2*-1 . 9DO+3 .D0mS1mS3/S22*33 .D0MS4/S22 


.N 

V- 


58 
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CALL  GiNOM! ’ F ’ , -1 , 2 , PE ASH) 

CALL  GRPRNTI-1 , X , 0 , Y , . 01 , 0 . ) 

Y»Y-DY 

CALL  GRTXT! ’PITCH  ’ ,20) 

CALL  GRNUM ( ’ F ' , - 1 , 2 , P I TCH ) 

CALL  GRPRNTI-1. X,0,Y, .08,0.) 

Y*Y-DY 

CALL  GRTXT! ’ROLL  ’ ,20) 

CALL  GRNUM ( ' F ’ , -1 ,2, ROLL) 

CALL  GRPRNTI-1 ,X,0 , Y , .08 , 0 . ) 

Y*Y-DY 

CALL  GRTXK ’YAW  ’,20) 

CALL  GRNUM! ’F’ ,-1.2, YAW) 

CALL  GRPRNT(-1 ,X,0,Y, .08,0. ) 

Y*Y-DY 

CALL  GRTXT (’TILT  ’ ,20) 

CALL  GRNUM! ’ F ’ , -1 ,2, TILT) 

CALL  GRPRNT !-1,X,0,Y,.08,0.) 

Y=Y-DY 

CALL  GRTXTI’PHIR  ’ ,20) 

CALL  GRNUM! ' F ' , -1 ,2,PHR) 

CALL  GRPRNTI-1 ,X,0, Y ,  .08,0.) 

CALL  GRMOVE 1 0 . , 0 . ) 

DO  10  4*1 , 90 

T*FLOAT(J)/2.8 

W*T»(T*2 . 1 / 1 T ♦ 1 .  MM»2 

G*2 .NTlHT+3  . )  / 1  (T*1 .  IN (T+2 . ) ) 

CALL  GRDRAWIW.G) 

CALL  QRDRAWI1 . ,2. ) 

CALL  GRDRAW! .8,0.) 

DO  14  4*1 , NS 

CALL  QRSYM I SLONE ( 2 , J ) , 3 LOBE  1 3, J) , .1,0. ,42) 

CALL  GRTXT I’  ’ ,1) 

CALL  GRNUM! ’ I ’ ,-1 , IF IXISLOBE (1 , J) ) .DUMMY) 

CALL  GRPRNTI-1 ,SLOIE(2,4) , 0, SLONE  18, J) , .08,0. ) 

PREPARE  KURT08IS  GRAPH 

IF! . NOT .GDAUTO)  CALL  GRMODEI ’ALPHA’ ,0) 

IF IQDAUTO. AND  IZETA.EQ.O)  CALL  GRTEKG 
IF (IZETA.EQ.O. AND. .NOT. GDAUTO)  READII0R.18)  IJK 
FORMAT (A2) 

IF! IZETA.EQ.O)  CALL  GRMODE ( ’ ERASE’ ,8) 

IF  I IZETA.GT.O)  CALL  GRPLOT (8 . 8 ,0. , -1 ) 

CALL  QRAX 18(0,1. 8, 1.8, 8. ,0. ,0. ,1 . ,8,2, ’WIDTH’ .8,2) 

CALL  GRAX 18(8, 1. 8, 1. 8, T. 8, 80., -2. ,8. ,8, 2,’ KURTOS 1 8 ’ ,8,2) 
CALL  QR8CL (0. , -2 . , 1 . ,8 . I 

CALL  GRSXYM. , 8.8.X, Y) 

CALL  GRTXTIIUF.NR) 

CALL  GRPRNT (0,X,0,Y,.18,0.) 

CALL  GR8XY(8.8,8.8,X,Y) 

CALL  GR8XY 18 . 8 ,8 . 3 ,X ,DY) 

DY«Y-DY 

CALL  GRTXT (’FREQUENCY  ’,30) 

CALL  GRNUM! ’F’ ,-1,2, FREQ) 
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mo> 

1890> 

)700> 

171 0> 
1720> 
1780> 
1740> 
1750> 
1790> 
1770> 
1780> 
1780> 
1800> 
1910> 

1  V20> 
1880> 

1 140> 
1880> 
190O> 
1870> 
19«0> 
1I90> 
1900> 

1 110> 
1I2«> 
1990> 
1940> 
1950> 
199«> 
1970> 
1990>  C 
1990> 
2000> 

201 0> 
2020> 
20J0> 
2040>  18 
2080> 
2090> 
2070> 
2090> 
2090> 
2100> 

2 1 1 0> 
2120> 
2190>  19 
2140> 

21 80>  C 

2190> 

1170> 

2 1 90> 
2190> 
2200>  17 
220S> 
2210> 
2220> 
2290> 


>»  FILE  • 
10* 

20> 

*0> 

40> 

80> 

80> 

70> 


CALL  GRTXTI ’  GHZ’ ,4) 

CALL  GRPRNTI-1 , X.O.Y, .08,0.) 

YaY-DY 

CALL  GRTXT I *  POLAR  I ZAT ION  ’,80) 

CALL  GRTXT (POL, S) 

CALL  GRPRNTI-1 .X.O.Y, .09,0.) 

Y«Y-0Y 

CALL  GRTXT I '%  SHADOWED  *,20) 

CALL  GRNUMI’F’ ,-1  .2.PERSH) 

CALL  GRPRNTI-1 ,X,0,Y, .08,0.) 

Y»Y-DY 

CALL  GRTXTI 'PITCH  ' ,20) 

CALL  GRNUMI'F’ ,-1 ,2, PITCH) 

CALL  GRPRNTI-1. X,0,Y. .08,0.) 

YaY-DY 

CALL  GRTXT I ‘ROLL  ' ,20) 

CALL  GRNUMI'F' ,-1 ,2, ROLL) 

CALL  QRPRNTI-1 ,X,0,Y, .08,0.) 

Y«Y-DY 

CALL  GRTXT I 'YAW  ’ ,20) 

CALL  GRNUMI’F’ ,-1 ,2, YAW) 

CALL  GRPRNTI-1. X,0,Y, .08,0.) 

YaY-DY 

CALL  GRTXTI 'TILT  ’ ,20) 

CALL  GRNUMI’F’ ,-1,2, TILT) 

CALL  GRPRNTI-1, X,0,Y, .09,0.) 

YaY-DY 

CALL  GRTXTI’PHIR  ‘ ,20) 

CALL  GRNUMI’F’ ,-1,2,PHR) 

CALL  GRPRNTI-1, X.O.Y,. 09,0.) 

CALL  QRMOVE I 0 . , 0 . ) 

DO  18  Ja|, SO 
T IF LOATIJ 1/2.5 
W=T*(T*2.)/(T*1 .1**2 
Q:0 .NT* IT«4 . )/(T*2. ) MM2 
CALL  GRDRAWIW.G) 

CALL  GRDRAWI1 . ,8. ) 

W=1 . 

DO  10  Jal.SO 
WaW-.OI 
Z»4.HW-2 . 

ZN* (Z*2 . ) / 12 .ttW-Z) 

ZN2*ZN«ZN 

G«8 .«ZN2«IW»1 . / IZN29W) -8 . /ZN-1 . /ZN2) / 1 (ZN«1 . )«MZN*2 . ) ) 
CALL  GRDRAWIW.G) 

CALL  GRDRAWIO. ,-1.8) 

DO  17  Jal.NS 

CALL  GR8YM(8LO0E(2,J) , 8LO0EI4, J) ,.1,0. ,42) 

CALL  GRTXTI’  ’,1) 

CALL  GRNUMI’ I ’ ,-1 , IFIXI8L0BEI1 ,J)) , DUMMY) 

CALL  GRPRNTI-1 .8L0BE 12 , J) ,0, SLOBE (4, J) ,.09,0.) 

CALL  GRRE8ISAVE) 

CALL  GRHV('H’) 

RETURN 

END 


SYO : [ 1 1 1 , 1 2 1 SPAWNS . FTN  27-SEP-98  14:89:42  <<< 

SUBROUTINE  SPAWNS 
INTEQERN2  EXSTATI9) 

DATA  TASK/9RMCR. . ./ 

CALL  SPAWN  I  TASK,,,  1 1 ,  .EXSTAT, ,  ’  RUN  LOSE/TASMLOSE’  ,19, , 
CALL  ST0PFRI11) 

RETURN 

END 


IDSW) 
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>»  FILE  >  SYO:  Mil,  1 21LOBE .FTN  27-SEP-83  14:38:4*  <<< 


30>  C 

L 

000 

8BSB 

EEEEE 

40>  C 

L 

0  0 

B  B 

E 

S0>  C 

L 

0  0 

BBBB 

EEEE 

•  0>  C 

L 

0  0 

8  B 

E 

70>  C 

LLLLL 

000 

BBBB 

EEEEE 

COMPUTE  RCS  LOBE  STATISTICS 


CMMNNMNNNNNNlHINIHHINNMNNMItNNNNMMNttMMNMMItNNNMNNMttNNNNIHt 

c 

REALM  SLOBE44, 1S1  >  ,*10(4,180)  .NF.NDMSI  .31  .SAVE  1301 
REALM  THETAI4I .ALPHA  14) .DELTAI4) .SCALE (•) ,3(2,2) ,CI2.2) 
REALM  DI2.2.2.2) , AMN 1 2 1 , DRS 1 2 ) , TMP ARR (2,2) .QSIZE(S) 
INTE8ERM  OUNT.RR.SS 
LOO  I  CALM  ODAUTO 
BYTE  BUFISO) 


DIDIX) >44877. KCOSM . S707983-ATANIX/6SS .31 1 


IOW*  2 

CALL  ASNLUNI IOW, 'Tl ' ,0) 
I0L*1 

CALL  ASNLUNI lOL.'SY' ,0) 
I0T>9 

CALL  ASNLUNI lOT.'SY' .0) 


OPEN  I UN  IT* I OL, NAME*' LOBE. TMP’ ,TYPE>'0LD' .RECORDS  I ZED  31 . 
1  ACCESS* 'DIRECT' .READONLY, SHARED) 


READ (I0L ' 1 )  I  OUT .KLOBE ,OL I  NT . NRCS , NLOBE , THETA .ALPHA , DELTA , 

1  BETA , SCALE , I ZETA . I PCODE.WLEN , LBUQ.NB, BUF . ODAUTO , 

2  FREQ. POL . PERSH . P I TCN . ROLL , YAW, T I L T , PHR , OS  I ZE 
LBUQ4*LBUO . AND . 4 

LBUQ8*LBUQ . AND . 8 


IF  I  DELTA  13) .EQ.O. )  WRITE(IOW.S) 

FORMAT  I '  LOBE>  —  DELTA  HAS  NOT  BEEN  DEFINED’) 
IFIDELTAI3) .EQ.O. )  00  TO  33 


I F I LBUQ4 . NE . 0 )  WR I TE 1 2 , 227 )  I OUT .KLOBE , QL I  NT , NRCS , NLOBE , I ZETA , I PCODE . 
1  THETA, DELTA, ALPHA, SCALE, BETA, WLEN 
FORMAT(1X,7l4/!X,SFS.2/1X,SFS.2/!X,SFt.2l 


IO>IOW 

IFI IOUT.EQ. 1 I  00  TO  4 


INITIALIZE  LOBE  FOR  PRINTED  OUTPUT 


IFI IPCODE. EQ.O)  00  TO  7 
10*3 

CALL  ASNLUNI 10, 'OD' ,0) 
00  TO  7 


INITIALIZE  LOBE  FOR  QRAPHICAL  OUTPUT 
IOQM 
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S90> 

sso> 

970* 

5S0> 

S90> 

C 

000> 

C 

B10> 

C 

620k 

C 

0S0> 

C 

640» 

C 

660k 

C 

660k 

C 

670k 

680k 

7 

660k 

222 

700k 

C 

710k 

C 

720k 

730k 

740k 

760k 

760k 

770k 

7S0k 

760k 

C 

800k 

c 

810k 

c 

620k 

c 

•30k 

c 

840k 

c 

880k 

•60k 

c 

•70k 

660k 

660k 

600k 

610k 

21 

620k 

c 

630k 

640k 

690k 

660k 

22 

670k 

228 

660k 

660k 

C 

1000k 

C 

1010* 

1020k 

1030k 

c 

1040* 

c 

1060* 

c 

1060* 

1070* 

1080* 

1060* 

1100* 

c 

1110* 

c 

1120* 

1130* 

1140* 

1 

CALL  8R I NIT  10, IZETA, 100, 1 ER) 

CALL  8RHVCH’) 

CALL  ORPICIOSIZEM)  ,aSIZE(2) .OSIZE(S)  .8SIZEU)) 
IFI IZETA. EQ. 01  CALL  8RDEV I’ TEK ’ , ’ HARD ’ , I ER ) 


LBUO  DEFINITIONS: 

LRU8a1s>T0L0BE  TRANSFORMATIONS 
LBU8b2s>T0L0BE  OUTPUT 
LBU0b4s> INPUT  A  LOBE  OUTPUT 
LBU8iS»0LINT  OUTPUT 

NANOa (DELTA  1 2 ) -DELTA (1 ) ) /DELTA  1 3  >  >1  .0001 

IFI IOUT. EO.O. OR. IZETA. NE.O)  WRITE! 10,222)  NAN8 

FORMAT  ( ’  0>>  COMPUTE  RCS  LOBE  STATISTICS  FORMS.’  THETA  ANOLES  «’) 

COMPUTE  CONSTANT  TERMS 

SNBETbSINIBETA) 

CSBET1 aCOSIBETA) *1 . 

SOL> .2BBBES 
P|aS.f4tSS28 
TWOP I *P 1*2 . 

RADsPI/110. 

IF ( IOUT .EQ. 1 >  CALL  8RSCL (SCALE) 1 ) .SCALE  12) .SCALE IS) .SCALE (4) ) 


PERFORM  ENTIRE  SET  OF  LOBE  STATISTIC  COMPUTATIONS 
(READ  ENTIRE  TEMPORARY  FILE  EACH  TIME) 

FOR  EACH  REOUIRED  DELTA  ANOLE 

AN8*DELTAI1) 

IF IAN8.8E . THETA  1 1 ) )  80  TO  22 

ANaaAN8sDELTAISI 

NAN8SNAN8-1 

IF INAN8. LE . 0)  STOP  ’LOBE  ERROR’ 
ao  TO  21 

F IRSTA>ANa 

SLOPEaFLOAT (NRCS-1 I /AMAX1 ( . 01 . THETA (2) -THETA (1)1 
CONST>1 . -SLOPEKTHETA I 1 ) 

IFILBUa4.NE.0l  WRITE I2.22M)  SLOPE .CONST. ANR 
FORMATI’  SLOPE, CONST. ANaa > ,3F1 0.3) 

DO  100  LL>1 , NAN8 

DETERMINE  NEXT  INDEX  INTO  LOBE  ARRAY 
IF  I AND. 8T. THETA  12 ) )  80  TO  101 
La  I F I X I  SLOPE BANS ‘CONST* .001) 

COMPUTE  BUM  AND  STANDARD  DEVIATION 

NlSa-3 

SUMaO. 

STDaO. 


NLB*NLB*4 

IFINLB.8T.NL0BE)  DO  TO  I 

READI IOL’NLB*1 1  ISLOBE 1 1 ,J) ,Jal , LI 


THE  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

laurel  Maryland 


1150>  I F ( SLOBE II, LI. EQ. 0.1  QO  TO  1 

11«0>  DO  19  K*2.4 

1170>  19  READ! I0L'NLB«K)  (SLOBEIK, 9) , 9»1 ,L) 

1 190>  C 

1 190>  IF ILBUQ4.NE . 0)  «H I TE 1 2 , 22 1 )  I5L09E 1 1 , LI , I *1 , 4 J 

1200>  211  FORMAT  I’  RCS.POSN*’ .4E11 .3) 

1 2 1 0>  C 

1220>  N*N*1 

1290>  DO  2  J>1 ,4 

1240>  2  918(9 ,N) sSLOBE (J . LI 

12S0>  C 

129©>  SUM*SUM*SI8I1 ,N) 

1270»  STD*STD*SI8( 1  , N)**2 

1280>  IF (N. LT. 1901  80  TO  1 

1290»  IF  I99.EO.  1901  1RITEI2.9) 

1300>  9  FORMAT ( ’  RCS>  —  LOBE  WILL  USE  ONLY  THE  FIRST  180  COMPONENTS’) 

1310»  C 
1 320>  C 

1330  3  IFIN.LE.il  80  TO  100 

1340>  IFMIOUT. EQ.O. OR.  I2ETA. HE. OJ  .AND.LL.EQ.il  WRITEI2.24)  N 

1350  24  FORMAT  11  OX , ‘ USING’ ,14,’  SUBCOMPONENTS  <<’ I 

1390  STD*SQRTISUM«SUM>STO) 

1370  IFI ILBU841 .NE .01  WRITEE2 ,2231  N, SUM. STD 

1990>  223  FORMAT ( ’  N.SUM.STD*' , IS.2E11 .3) 

1390  C 
1400  C 

1410  C  COMPUTE  R  C  S  LOBE  I 

1420  C 
1430  C 

1440  8*0. 

1450  N13N-1 

1490  C 

1470  DO  5  I31.N1 

14SO  M3|*1 

1490  DO  5  J3|1  ,N 

1500  5  a*G+$l8(1 , I IttSIQM ,9)N( (810(4, 1 1-518(4, 9)>«CS9ET1« 

1510  9  (SIS  12, 1 1-SI8I2,  J))MSNBET)Mtt2 

1820  NF*SORT(Q)/(SOLMSTD) 

1530  IFI  ILBU84)  .NE  .0)  MTR I  TE  ( 2 . 2241  NF 

1540  224  FORMAT!'  NF>',4E11.3) 

1550  NF3NFP1.E9 

1590  C 
1570  C 

1390  C  COMPUTE  R  C  S  LOBE  II 

1S90>  C 
1900  C 
1910>  C 

1920  C  COMPUTE  A  MATRIX  (A(1 ,2)  *A(2 ,1 1) 

1930  C 

1940  A1130. 

1950  A12*0. 

1990  A22>0. 

1970  C 

1990  DO  9  1*1  ,N1 

1990  1 1  b  |  *1 

1700  DO  9  J*I1  ,N 

1710  C 

1720>  71*818(1, ll«SI8(1, 91 

1730>  C 

1740>  T2> ISI8 12, I ) -918(2,9) MCSBET1 -151814,11-510(4,9) ImSNBET 
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1790>  i 

1790> 

1770> 

17*0> 

1790> 

1300> 

1I10> 

1I20> 

1930> 

1M0> 

mo> 
i»eo> 
1«70> 
1990> 
1990> 
1900» 
191 0> 
1920> 
1930> 
1940> 
1990> 
1990> 
1970> 
1900* 
1990> 
2000> 
201 0> 
2020> 
2090> 
2040> 
2090> 
2090> 
2070> 
2090> 
2090> 
2100> 
211 0> 
2120> 
21«0> 
2140> 
21»0> 
2190> 
2170> 
2190> 
2190> 
2200> 
2210> 
2220> 
2290> 
2240> 
22S0> 
2290> 
2270> 
2290> 
2290> 
2900> 
2I10> 
2220> 
2«30> 
2940> 


T3* IS  IQ 1 3, 1 1 -S IG 13 , J I )*CSBET1 

A1 1 *A1 1 ♦T1PT24T2 
A12*A12*TH*T2»T3 
A22*A22*T1»T3«T3 

CONTINUE 

A21SA12 


COMPUTE  FI DELTA I  VS  DELTA 

DELT3360./DELTAI4) 

DELT*2. 

DEL*0. 

IDEL*1 

IFILBUO4.NE.0)  WH I TE (2,225)  DELT.A1 1 ,A12 ,A21 . A22 
FORMAT  I ’  DELT** ,E11 .3/*  A*’ ,2E11 .3/3X.2E11 .31 

IF (DEL .OT .339 .9)  80  TO  11 
NDI IDEL , 1 1 *DEL 
D1*SINIDEL«RAD) 

D2  * -COS  I DE  LPRAD ) 

T1*A11«01«A12ND2 

T2*A21«01*A22«D2 

A*D1»T1*D2«T2 

ND I IDEL , 2 ) (SORT (A) / I WLENNSTD) 

DEL*DEL*DELT 
IDEL*IDEL«1 
80  TO  10 

NDI IDEL,1 1*360. 

NDI I  DEL, 2) *ND II ,21 

DO  19  J*1 , IDEL 

ND(J,1)«AMOD(NDU,1)«90.,9S0.l 
IFIIOUT.EQ.il  80  TO  15 


WRITE! 10,121  AN8.NF 

FORMAT  I ’ORCS  LOBE  OUTPUT’ ,5X, ’ THETA* ’ ,FS . 1 , ’ 
9  SX, ’ DELTA’ .IX, ’ND’ /I 

DO  13  J*1 , IDEL 

WRITE! 10,14)  ND I J , 1 1 , NO (J , 2 1 
FORMAT I1X,F9.2,F10.9) 

80  TO  100 

8RAPH  NF  AND  1/ND 

CONTINUE 

FMAX*0. 

FM IN* 1000. 

DO  IS  J*1 , IDEL 

NDU,  21*1  ,/INOIJ,  2INRAD) 

FMAX*AMAX1 IFMAX.NDIJ ,2) J 


NFI8HZI* ' ,FS .3/ 


*,  *  R.  ^  *  ' 
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2990 

2360 

2970 

2310 

2390 

2400 

2410 

2420 

2430 

2440 

2490 

2490 

2470 

2490 

2490 

2900 

2910 

2920 

2930 

2940 

2990 

2900 

2970 

2990 

2990 

2000 

2010 

2020 

2030 

2040 

2090 

2000 

2070 

2700 

2710 

2720 

2730 

2740 

2790 

2700 

2770 

2700 

2790 

2000 

2010 


10 

C 


C 

17 


C 

100 

101 

C 


C 

c« 

c 

c 

c 

c 

c 


2030 

2040 

2000 

2000 

2000 

2010 

2920 

2940 


IF (FMAX .EQ.MD I J ,2)1  FANQ*ND ( J , 1 1 
FM IN: AM  INI IFMIN.NDI J ,211 
MDIJ.2lcNDIJ.2i 
CONTINUE 

IF IFANG.QE . 190 . )  FANG=FANG-1 90 . 

I F (FANG . QT . 90 . )  FAN0«FANG-1 90 . 

CALL  ORRXV (ANCI.O.  ,X, Y) 

Y*8SIZE (2) - .99 
CALL  SRSAV I  SAVE ) 

CALL  8RINITI1 ,0.0. lERi 
IFILL.NE.il  00  TO  17 

CALL  QRTXTI ’NF I0HZ) * ,0) 

CALL  ORPRNT I *1 ,2 . 1 , -1 , Y* . 93, . 1 , 0. i 
CALL  ORTXTI ’MAJOR’. 0) 

CALL  ORPRNT 1*1 ,2.1 ,-1 ,Y*.39, .1 ,0. ) 

CALL  ORTXTI  MINOR’ ,0) 

CALL  ORPRNT(«1,2.1.-1 ,Y*.23. .1 ,0.1 
CALL  ORTXTI • ANGLE ’ .0) 

CALL  ORPRNT («1,2.1,-1,Y«.09,.1,0.) 

CALL  ORNUMI’F’ .-1.3.NFI 

CALL  ORPRNT I 0,X,-1 ,Y+.93,.10,0.) 

CALL  GRNUMI’F’ .-1,2, FMAX) 

CALL  ORPRNT 10 ,X,-1 , Y« .39, . 1 ,0. ) 

CALL  ORNUMI’F’ ,-1,2, FMINI 
CALL  QRPRNTIO.X, -1 , Y« .23, .1 ,0. ) 

CALL  QRNUMI’F' .-1,2, FANG) 

CALL  ORPRNT I 0,X,-1 ,  Y+ . 03 , . 1 , 0. ) 

CALL  ORCPI ’POLAR’ ) 

Y*03IZE(2i 

CALL  ORP ICIX.Y ,X*1 . ,Y+1 . ) 

CALL  ORSCL (0. ,0. , (SCALE  13) -SCALE 1 1 ) )/6 . 9,300. ) 

CALL  QRPICIO. ,Y-.29,11.,Y«.9) 

CALL  QRVI ND I - 1 DE L , ND 1 1 ,2) , ND 11,1), IDUM.O.DUM.DUM, IDUM) 
CALL  ORUPPN 
CALL  GRRESISAVE) 

ANQ* AN8+DE  LTA 1 3 ) 

CONTINUE 

IFIOLINT.EQ.O)  CALL  EXIT 


COMPUTE  OLINT  LORE  III 


OPEN I UN  I T« I OT , NAME ■’ LORE . SCR ’, TYPE ■ ’ SCRATCH ’ ) 

DO  30  IALP>1 ,91,90 
PLOT 

IPI IOUT.EQ.O)  00  TO  90 

IF  I .NOT .ODAUTO)  CALL  QRMODE I ’ ALPHA' ,0) 

IF! IZETA.EQ.O.AND.ODAUTO)  CALL  GRTEKG 

IFIIZETA.EQ.O. AND.. NOT. ODAUTO)  READ! ION, 471  UK  > 


66 


1 
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944» 

47 

FORMAT IA2) 

950> 

IF ( IZETA.8T.0I  CALL  QRPLOTC11 . ,9 .5 , -1 ) 

960> 

CALL  QRINITM  ,0,0,  IERI 

•70> 

C 

900> 

50 

ALP : FLOAT ( 1  ALP-1 ) 

990> 

REWIND  IOT 

000> 

C 

010> 

C 

LOOP  ON  TNETA  INCREMENTS 

020> 

ANQ*THETA(1) 

090> 

DO  95  Lal.NRCS 

040 > 

IF (ALP .EQ.O. )  QO  TO  29 

0S0> 

c 

000> 

READI IOT, 201  N , XBAR , YBAR , SUM , B 

070> 

IFIQL INT.EQ. 1 .OR.N.LE. 1 )  00  TO  97 

0«0> 

00  TO  99 

090> 

c 

100> 

29 

NLBa-9 

tio> 

XBAR*0. 

120> 

TBAR=0. 

190> 

SUM>0. 

140> 

STD»0. 

150> 

N=0 

100> 

c 

170> 

91 

NLB*NLB*4 

1f0> 

IF INLB.QT.NLOBE)  00  TO  90 

100> 

READ!  I0L’NLB*1 )  ISLOBEM  ,4)  ,J*1  ,LI 

200» 

IFISL0BEI1  , L)  .EQ.O)  00  TO  91 

210> 

DO  92  Ks2,4 

220> 

92 

READI IOL ' NLB+K)  ISLOBEIK, J) ,  J*1 , L 9 

290> 

C 

240> 

N*N*1 

290> 

DO  99  Ja1,4 

260> 

99 

8IQIJ,N)sSL0BEU,L) 

270> 

C 

COMPUTE  ECaifaa] 

290> 

SUMsSUMtSIQM  , N) 

290> 

C 

COMPUTE  SUMIaifaaNgratk-XI 

»00> 

XBAR>XBAR«SI0I1 ,N)n9IQI2,N) 

310> 

c 

COMPUTE  SUMIai«na«9raak-Y) 

*20> 

YBAR*YBAR*SIQ(1 ,N)ftSI0(9,N) 

990> 

00  TO  91 

940> 

c 

9S0> 

90 

IF  IN. EQ.O)  WRITEI IOT, 20)  N 

«00> 

IF  IN. EQ.O)  00  TO  97 

970> 

C 

COMPUTE  fraak-Xbar 

990> 

XBAR «XBAR/ SUM 

990> 

C 

COMPUTE  fraak-Ybar 

400> 

YBARaYBAR/SUM 

410> 

C 

420> 

STDaO. 

4«0> 

IF IOL INT.EQ. 1 .0R.N.EQ.1 )  WRITEI IOT, 20) 

440* 

IFIQIINT.EQ. 1 .OR.N.EQ. 1 )  00  TO  97 

4S0> 

C 

490> 

BI1 ,1 )*0. 

470> 

Mil ,2) *0. 

490» 

9(2,1)10. 

490> 

B(2,2)a0. 

900> 

C 

«10> 

N1«N-1 

920> 

DO  94  lal.NI 

»90> 

11*1*1 

, XBAR, YBAft 
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3340s 

3550> 

C 

3960k 

3370> 

33B0> 

3990> 

C 

3600> 
301 0> 
3620> 
3630> 
3640k 

34 

3990> 

3660> 

20 

3070> 

C 

3680» 

39 

3690> 

3700> 

C 

371 0> 

C 

3720k 

c 

3730> 

3740> 

3790> 

3760> 

226 

3770 

3760k 

3790> 

C 

3600k 

37 

3610k 

3620k 

3630k 

C 

3640k 

3690k 

3660k 

C 

3670k 

36 

3660k 

C 

3990k 

39 

3900k 

C 

3910k 

3920> 

C 

3930k 

C 

3940k 

C 

3990k 

3960k 

93 

3970k 

3990k 

3990k 

4000k 

4010k 

94 

4020k 

99 

4030k 

4040k 

C 

4090k 

C 

4060k 

c 

4070k 

91 

4060k 

C 

4090k 

4100k 

4109k 

4110k 

4120k 

DO  34  Jail  ,M 

T 1 =  S I Q 1 1 , 1 1  MS  I Q 1 1 ,  J) 

T  2  =  S I Q ( 2 , l)+$IQ(2,J)*2. ttXBAR 
T3*S IQ  13 , I )*SIG(3,J)-2. NYBAR 

B(1 ,1)=BI1 ,1)*T1»T2*T2 
B(1 , 2 1 *B ( 1 , 2) ♦T1HT2NT3 
8(2,1)88(1,21 
B(2,2)sB(2,2)*T1«T3»T3 
CONTINUE 

WRITE! IOT.20)  N, XBAR.YBAR, SUM, B 
FORMAT 1A2.7A4) 

AlsCOS(ALP) 

A2=SIN(ALP) 

COMPUTE  GLINT  STANDARD  DEVIATION 

BBaB( 1 , 1 ) MAI  MAI *B (2 , 1 ) H A  2  *»  A 1 *B ( 1 ,2)»A1*A2*B(2 ,2)«A2*A2 
STD=SQRT IBB) / (2 .ttSUM) 

IF (LBUGI . HE .0)  WR I TE ( 2 , 226 )  N , SUM , STD , XBAR , YBAR , ANG , B 
FORMAT! ’ON, SUM, STDs' , I3.2E11 .3/ 

1  ’  XBAR , YBAR , ANQs ’ , 3E1 1 . 3/ 

1  ’  B= ' , 4E1 1  .31 

NOIL, 1 ) >ANG 
ND( L , 3 ) sSTD 

IFIALP.EQ.BO. )  GO  TO  36 

ND(L ,2) sXBAR 
GO  TO  39 

ND(L,2)3YBAR 

ANQ= ANG+THETA ( 3 ) 

I F ( IOUT.EQ.1 )  GO  TO  $1 

PRINT  GLINT  DISPLACEMENT  VECTOR 

WRITE! 10,33)  ALP 

FORMAT  I 'ORCS  GLINT  OUTPUT  ALPHAS’ ,F6.1/ 

1  ’  ANGLE  DISPLACEMENT  DIS-S  DISCS’ /I 
DO  94  LPH.NRCS 
T1 *ND I LP , 2 ) -ND I LP , 3 ) 

T2*ND(LP,2)*ND(LP,3) 

WRITE  1 10,93)  NDILP.1 ) ,ND(LP,2) ,T1 ,T2 
FORMAT IF6.1 . F1 1 .3 , F 1 0 . 3 , FI  0 . 3 ) 

GO  TO  92 

GRAPH  GLINT  DISPLACEMENT  VECTOR 

I F ( IZETA.EQ.O)  CALL  GRMODE I ’ ERASE ' ,3) 

CALL  QRPL0TIQSIZEI3) , QS I ZE (21 , 0 ) 

CALL  GRPL0TIQSIZEI1 ) .QSIZEI2)  ,1 ) 

DAXsQS I ZE ( 3 ) -GS I ZE  1 1 ) 

CALL  GRAX ISIO.GSIZE 1 1 ) ,GSIZE (2) - . 99 , DAX , 0. , SCALE  1 1 ) .SCALE (3) , 
1  I F I X ( SCALE ( 5 ) ) , I F I X ( SCALE ( 6 ) ) , ’OBSERVATION  ANGLE’, SB. 1) 
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4ieo>  c 
«m> 

4170> 

4  190> 

4190>  C 

4200> 

4210> 

4220> 

4230> 

4240> 

4250> 

4200> 

4270>  C 

42I0>  C 

4290>  C 

4S00> 

4310> 

4320> 

4330> 

4340> 

4990> 

4360> 

4370> 

4390> 

4390> 

4400> 

4410> 

4420> 

4430> 

4440> 

4450> 

4490> 

4470> 

4490> 

4490> 

4S00> 

4910> 

4S20> 

4S30> 

4S40> 

4»50> 

4S60> 

4570> 

4»90> 

4990> 

4900> 

461 0> 
4620> 
4930> 
4640» 
4990>  C 
4660>  C 
4670> 
46I0> 
4990>  C 
4700>  C 
471 0> 
4720> 
4730»  70 
4740>  C 


DAYsQS I ZE 1 4 ) -QS 1 ZE 1 2 ) 

CALL  GRAX I S (3 ,GS I ZE 1 1 ) ,GS I ZE (2) , DAY ,90. , -360. ,360. ,5.2. 
1  'DISPLACEMENT  (DEG)’, 99, 21 

IF  I IZETA.EQ.O)  CALL  GRMODE I 'ALPHA* ,2) 

CALL  GRTXTIBUF.NB) 

CALL  GAPRNTI0,5.5,-1 ,7. , .15,0.) 

CALL  GRTXTI’ MONOSTATIC  GLINT  FOR  ALPHAS’ ,01 
CALL  QRNUMI *F’ ,-1 ,1 ,ALP) 

CALL  QRPRNTIO, S .5,-1 ,6 .75 , . 1 5 ,0. ) 

IF! IZETA.EQ.O)  CALL  GRMODE! 'ALPHA' ,3J 

DISPLAY  RADAR  PARAMETERS 


XS9.0 
Y=7 . 

DVs.2 

CALL  GRTXTI 'FREQUENCY  ’,30) 

CALL  QRNUMI ’ F ’ ,-1 ,2, FREQ) 

CALL  GRTXTI’  GHZ’ ,4) 

CALL  GRPRNTI-1 ,X,0,Y, .12,0.) 
YsY-DY 

CALL  GRTXTI  POLARIZATION  ’,50) 
CALL  GRTXTIP0L.3) 

CALL  QRPRNTt-1 ,X,0,Y, .12,0. ) 
YsY-DY 

CALL  GRTXTI  X  SHADOWED  *,20) 
CALL  QRNUMI ’F’ ,-1,2,PER8H) 

CALL  GRPRNTI-1 ,X,0,Y, .12,0. ) 
YsY-DY 

CALL  GRTXTI 'PITCH  ’ ,20) 

CALL  GRNIIMI ’ F ’  ,-1,2, PITCH) 

CALL  GRPRNTI-1 ,X,0,Y, .12,0.) 
YSY-DY 

CALL  GRTXTI’ ROLL  ’,20) 

CALL  GRNUMI ’ F ’ , -1 ,2, ROLL) 

CALL  GRPRNTI-1, X,0,Y, .12,0.) 
YSY-DY 

CALL  GRTXTI ’YAW  ’,20) 

CALL  QRNUMI ’F’ ,-1,2, YAW) 

CALL  QRPRNT 1-1, X,0,Y,. 12,0.) 
YSY-DY 

CALL  GRTXTI ’TILT  ’,20) 

CALL  GRNUMI ’F' ,-1,2, TILT) 

CALL  GRPRNTI-1 ,X,0,Y, .12,0.) 
YsY-DY 

CALL  QRTXTI’PHIR  ’  ,20) 

CALL  QRNUMI ’ F ’ ,-1 ,2 ,PHR) 

CALL  GRPRNTI-1. X,0,Y, .12,0.1 


CALL  QRMOVE INDI1 , 1 ) .DIDINDI 1 ,2) ) ) 

DO  70  JRCSS2.NRCS 

CALL  QRDRAWINDI JRCS , 1 ) , D IDIND I JRCS , 2) ) I 
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4790> 

4760> 

C 

4770» 

C 

4790 

4790 

4900> 

49 

4910 

4920> 

4930 

99 

4940> 

4990 

4990 

C 

4970 

C 

4990> 

c 

4990> 

c 

4900 

c 

4910> 

c 

4920> 

c 

4930 

c 

4940 

92 

4990 

4990 

c 

4970> 

c 

4990 

4990 

9000> 

c 

9010> 

9020 

9030> 

9040 

9090> 

9090> 

9070> 

c 

9090 

41 

9090> 
9100 
S110> 
9 1 20> 
9130 

42 

9140> 

C 

9 1 9  0  > 
9190> 
9170> 

49 

9130> 

C 

9190> 

9200> 

C 

82 1 0  > 
9220> 

c 

9280 

9240> 

9290 

c 

9290 

40 

S270> 

C 

9280> 

S290> 

C 

9100 

8S10> 

C 

9320 

S990> 

9340 

IFtQLINT.EQ.il  QO  TO  S2 

GRAPH  DISPLACEMENT  ♦  /-  STD 

CALL  QRM0VEINDI1 .1) .DIDINDI1 ,2)*NDI1 .Sill 
DO  49  J-2'NACS 

CALL  QRDRAWINDIJ ,1 1 , DIDlNDtJ , 2 ) ♦ ND  <  J  ,S)  1 1 
CALL  QRMOVE INDt 1 , 1 1 , DID IND( 1 , 2) -ND ( 1 ,3)1) 
DO  99  J*2 , NRCS 

CALL  QRDRAWtNDtJ.il .DIOtNDIJ, 2) -NDtJ. 31)1 
CALL  QRUPPN 
CALL  QRRES(SAVE) 


COMPUTE  QLINT  LOBE  IV 


REPEAT  STD  CALCULATIONS  FOR  ELLIPSES 

LOOP  ON  OELTA  INCREMENTS 
ANQ=FIRSTA 
DO  49  LL&1.NANQ 

DETERMINE  NEXT  INDEX  INTO  LOBE  ARRAY 
IF(ANQ.QT.THETA(2) I  QO  TO  39 
L= I F I X ISLOPEHANQ+CONST* . 001 1 

NLBa-3 

XBAR30. 

YBAR*0. 

SUMsO. 

STD»0. 

N*0 

NLBsNLB«4 

IFtNLB.QT.NLOBEI  QO  TO  40 
READ! IOL’NLB+1 1  ( 8  LOBE ( 1 ,J) ,J  =  1 , LI 
I F ( S LOBE ( 1 , L I . EQ.O)  QO  TO  41 
DO  42  K*2 ,4 

READ! IOL ’ NLB+K)  (SLOBEtK. Jl .4*1 .LI 
N«N*1 

DO  43  J31.4 
8 1 Q  ( J , N I *SLOBE ( J , L I 
COMPUTE  Elaioma] 

SUM*SUM*S IQt 1 . N I 
COMPUTE  SUM  I  •  i  flma»9  r  ok  -X I 

XBARaXBARtSIQM ,N)«SIQ(2.N) 

COMPUTE  8UMIai9aia«9raak>Y) 

YBAR«YBARpSIQ(1  .NINSIQO.NI 
QO  TO  41 

IFtN.LE .1 1  QO  TO  49 
COMPUTE  9r««k*Xbar 
XBARiXBAR/SUM 
COMPUTE  sraak-Ybar 
YBAR8YBAR/SUM 


DO  49  MM*1 ,2 
DO  49  NNsi  ,2 
BtMM.NNHO. 
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5350 

5380 

5370 

5310 

5380 

5400 

5410 

5420 

5430 

5440 

5450 

5480 

5470 

5480 

5480 

5500 

5510 

5520 

0530 

5540 

5550 

5580 

5570 

5580 

5580 

5600 

5010 

5620 

5030 

5040 

5050 

5000 

5070 

5000 

5060 

5700 

5710 

5720 

5730 

5740 

5750 

5700 

5770 

5700 

5700 

5000 

5010 

5020 

5030 

5040 

5050 

5000 

5070 

5800 

5000 

5000 

5010 

5020 

5000 

5040 


CIMM,0N)*0. 

DO  48  RR*1  ,2 
DO  48  SS*1.2 

46  DIMM, 00, RR, 551*0. 

C 

01*11-1 

DO  44  1*1 .01 
11*1*1 

DO  44  J*I1 ,0 
C 

T1 *S IQ (1 . I »»S IQ ( 1 . J ) 

T2 *51012, I )*SIGI2, JI-2.NXBAR 
T3*8IQ(3, I ) *810(3, J) -2. NVBAR 

T4*( 1510(4, I) -51014, J)  UCSBET1* <51012, I ) -01012, JI)»*0IETJ»»2 
T5* 1 5 10 12, 1 ) -5 10 (2 . J ) IMCSBET1 - f 510(4 , I )-SIQI4,JI IO50BET 
T6* 1510(3 , 1 1 -S 10(3 , J I IttCSBETI 
C 

TMPARR (1,11 *T2nT2 
TMPARR 11,2) *T2MT3 
TMPARR 12,11 *TMPARR (1,21 
TMPARR (2 ,21*13013 
DO  44  MM*1 ,2 
DO  44  00*1,2 
TEMP*TlOTMPARRIMM,00) 

0 (MM, 00) *0 IMM , 00) *TEMP 
C (MM.OO) *C(MM,00) *TEMPoT4 

D (Ml , 00 , 1 , 1 ) *D (MM , 00 . 1 , 1 ) *TEMP0T50TS 
D (MM , 00 , 1 , 2 ) >D (MM , 00 . 1 , 2 ) «  TEMPOT5OT0 
D(MM,00 ,2 , 1 ) *D(MM,00, 1 ,2) 

D (MM , 00 , 2 , 2 ) *D IMM , 00 , 2 , 2 1 *TEMPOT0OT6 
44  COOT  I OUE 

C 

AlsCOSIALP) 

A2*SI0IALP) 

C 

C  COMPUTE  QLIOT  STAODARO  DEVIATION 

BB*B(1 , 1 INA10A1 *0(2,1 10A2NA1 *011 ,2)OA1oA2*RI2 ,2)0A20A2 
STD*SQRT ( BB I / 1 2 . NSUM) 

CC*C 1 1 , 1 )MA1*A1 *C (2 , 1 )*A2*A1 *C 1 1 , 2)NA1MA2*C (2 , 2)NA2«A2 
C 

C  COMPUTE  OlOtlta.M 
C 

BF*SQRT(CC)/(2.ftSTD«SOL) 

0F*0F«1 .EO 
C 

IF (LBUQ8 ,0E .0)  VRITEI2.220)  0,5UM,8TD.XBAR,TBAR,B,C,D 
220  FORMAT (’00, SUM, STD*’ ,13, 2E 11.3/ 

1  ’  XBAR.YBAR*’ , 2E 1 1 .3/ 

1  ’  B*’ , 4E 1 1 .3/ 

1  '  C*’ ,4E11 .3/ 

1  41 ’  D«’ . 4E1 1 .3/)) 

C 

C 

C  COMPUTE  SHOT  LOBE  V 

C 

C 

AMNMIsAI 

AM0(2)*A2 

DELT*360./DELTA(4) 
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V 


DEL*0. 
f  DEL  =  1 


I F (DEL . GT .339 . 9 1  GO  TO  61 

MOIIDEL.DsDEL 

DRS ( 1 ) :S I N IDELttRAD) 

DRS (2 ) =-COS IDELmRAO) 


DDsO. 

DO  69  MM=1 ,2 
DO  69  NNB1.2 
DO  $9  RR=1,2 
DO  69  SSsI .2 

DD*DD+D  IMM,  NN ,  RR ,  SSIftAMN  IMMINJ 


IS ( SSI 


NDI IDEL ,2) *SQRT (DO) / (2 .MWLENftSTD) 


BELsDEL*DELT 
IDEL*IDEL+1 
GO  TO  60 


NDI IDEL.1 (=360. 

NDI IDEL,2)sND(1 ,21 


DO  68  J=1,IDEL 

NDI J , 1 ) sAMOD (NDI J . 1 ) *90. , 660. I 


IFIIOUT.EQ.il  GO  TO  6S 


PRINT  NF  AND  ND 


WRITE! 10,62)  ANQ.NF.ALP 

FORMAT  I ’ORCS  GLINT  OUTPUT ’ , 5A , ’ THETA* ’ , F6 . 1 , 1 
5X , ’ DELTA’ ,5X, ’ND’ ,5X. ’ALPHAS’ .F6.1/I 


NF I GHZ I > ' , F8 . 3/ 


DO  63  Jsl , IDEL 
TEMPH  .  / 1  ND  I J  , 21NRAD) 

WRITE  1 10,641  NDU.1)  .NDIJ.2)  .TEMP 
FORMAT (1X.F9.2.2F10.3) 

GO  TO  43 


GRAPH  NF  AND  1/ND 


CONTINUE 

FMAXsO. 

FMINslOOO. 

DO  66  J*1,I0EL 
ND  I J ,2) *1 . / INDl J , 2INRAD) 
FMAXsAMAX 1 1 FMAX . ND U , 2 1) 
IFIFMAX.E0.NDIJ.2n  FANQsNDU.1l 
FMIN3AMIN1  (FMIN.NDU  ,2) ) 

NDIJ , 2) sNDI J ,2)9.1 
CONTINUE 

I P (FANG. GE. 130. I  FANQsFANG-190. 
IF  IFANQ.QT.90. 1  FANQsFANG-190. 


CALL  GRRXYIANQ.O. ,X,Y) 
YsQIIZE 12) - .98 
CALL  QRINITI1 .0,0, IER) 
IFILL.NE.il  GO  TO  67 
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tsso>  c 
osoo> 

0970> 
9380> 
9990> 
6600> 
9010> 
6620> 
9930> 
9940>  C 
9090>  97 
6660> 
9679> 
6690> 
9990> 
9700> 
9710> 
9720> 
8730>  C 
9740> 
9790> 
6790> 
9770> 
67»0> 
9790> 
6900> 
0«10> 
5«20>  C 
9830>  C 
9I40>  49 
9I90>  39 
9990> 
9970> 
9990> 


CALL  QRTXT! ’NF(QHZ) ’ ,01 

CALL  QRPRNTI+1 .2 . 1 , -1 , Y« . 93 , . 1 ,0. ) 

CALL  QRTXT ( ’MAJOR’ ,01 

CALL  ORPRNT (♦! ,2.1 ,-1 , Y* .33 , . 1 ,0. ) 

CALL  8RTXTI ’MINOR’ ,01 

CALL  QRPRNT (+1,2.1,-1,Y*.23,.1,0.) 

CALL  QRTXT I ’ANQLE' .01 

CALL  QRPRNT (*1 ,2.1,-1,Y*.09,.1,0.) 

CALL  QRNUMI ’ F ’ , -1 ,3 ,NF) 

CALL  QRPRNT(0,X,-1 ,V*. 93.. 10, 0.1 
CALL  QRNUMI * F ’ . -1 , 2 , FMAX ) 

CALL  QRPRNT (0,X, *1 , Y* . 38  , .  1 ,0.) 

CALL  QRNUMI ’ F’ , -1 , 2.FMINI 
CALL  QRPRNT 10 ,X, -1 , Y* . 23 , . 1 ,0 . I 
CALL  QRNUMI ’ F ' ,-1,2,FANQ) 

CALL  QRPRNT 1 0, X, -1 , Y+ . 08 , . 1 , 0. ) 

CALL  QRCP I ’POLAR’ ) 

Y3QSIZEI2) 

CALL  QRPICIX, Y.X+1 . , Y*1 . ) 

CALL  QR8CLI0. .0. . (SCALE (3) -SCALE II ))/9 .9 .390. ) 

CALL  QRPICIX-1 . .Y-.28.X+1 . , Y*.9I 

CALL  QRVINDI-IDEL.NOII ,2) , ND 1 1 , 1 ) , IDUM.O.DUM.DUM, IDUM) 
CALL  QRUPPN 
CALL  QRRESISAVE) 


AN9=ANQt DELTA  13) 
CONTINUE 
CLOSE  I UN IT> I OL) 
CLOSE  I UN  ITS  I OT) 
CALL  EXIT 
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APPENDIX 

A  PLOTTING  COORDINATE  SYSTEM  FOR  THE  APL  RCS/STATISTICS  CODE 


Background 

The  following  remarks  are  applicable  to  the  APL 
RCS/Statistics  Code  described  in  this  document.  In 
consequence  of  the  code’s  sequential  development 
over  a  rather  long  time,  several  peculiarities  exist.  Be¬ 
cause  they  are  not  explicitly  pointed  out  in  the  body 
of  the  text,  their  significance  may  be  overlooked, 
particularly  in  regard  to  recent  demands  on  the  code 
that  were  not  anticipated  during  its  development. 

The  pecularities  have  mainly  to  do  with  coordinate 
systems.  In  early  work,  a  master  coordinate  system 
(using  unprimed  variables)  was  chosen  arbitrarily 
with  the  z  axis  pointing  forward,  the  x  axis  to  port, 
and  the  y  axis  up.  That  master  system  was  connected 
to  the  component  coordinate  systems  by  transforma¬ 
tions  called  TRANS  in  the  code. 

Since  a  comparison  of  theoretical  calculations  with 
experimental  data  is  always  desirable,  and  since  the 
RATSCAT  facility  is  perhaps  the  chief  supplier  of 
such  data,  it  became  advisable  to  introduce  the  RAT- 
SCAT  standard  into  our  master  system.  RATSCAT 
uses  the  convention  z  axis  forward,  x  axis  up,  and  y 
axis  pointing  to  starboard.  Rather  than  change  the 
code  structure,  the  RATSCAT  master  system  was 
connected  to  our  master  system  by  a  simple  linear 
transformation.  The  shadowing  of  one  component 
by  another  was  now  defined  in  the  RATSCAT  system 
(the  instructions  are  called  VALID  in  the  code).  Since 
the  target  aspect  modifications  of  pitch,  roll,  yaw, 
and  tilt  are  also  available  in  the  RATSCAT  system  it 
was  decided  to  introduce  these  transformations  into 
the  code. 

Statistical  information  about  the  lobing  structure 
of  an  RCS  pattern  can  be  regained  provided  that  the 
component  scatterer  locations  can  be  specified. 
When  this  development  was  completed,  the  locations 
were  specified  by  giving  the  position  of  the  origin  of 
each  component  system  in  the  RATSCAT  system 
(this  instruction  is  called  ORGN  in  the  code). 

Now  we  have  the  requirement  to  describe  RCS 
over  an  arbitrary  planar  cut  or  over  a  conical  cut 
about  an  arbitrary  cone  axis. 


Modifications 

Because  shadowing  and  component  origins  are 
specified  in  the  RATSCAT,  or  radar-fixed,  coor¬ 


dinate  system,  the  pitch,  roll,  yaw,  and  tilt  transfor¬ 
mations  are  disallowed  for  other  than  small  angular 
excursions.  But  they  are  exactly  the  transformations 
one  would  like  to  use  to  satisfy  the  requirement  for 
arbitrary  planar  and  conical  cuts.  The  best  way  to 
satisfy  the  requirement  is  to  eliminate  the  original 
pitch,  roll,  yaw,  and  tilt  capability  and  replace  it  with 
a  new  system  called  PLOTTING.  The  target-fixed 
system  stands  fast  while  the  plotting  system  axes 
(originally  coincident  with  the  radar-fixed  axes)  are 
“pitched,”  “rolled,”  or  “yawed”  to  a  new  posi¬ 
tion.  Thus,  the  original  terminology  is  retained  but 
with  the  meaning  just  defined.  (Actually,  one  can 
think  of  the  target  itself  being  pitched,  rolled,  or 
yawed,  but  in  an  angular  direction  opposite  to  that  of 
the  axes  movement.) 

In  addition  to  this  change,  the  original  unprimed 
coordinate  system  is  eliminated.  Henceforth, 
TRANS  expresses  a  transformation  from  the  primed 
to  the  RATSCAT  system  (in  a  manner  exactly 
analogous  to  the  way  TRANS  previously  expressed  a 
transformation  from  the  primed  to  the  unprimed 
system). 


Transformations 

There  are,  therefore,  three  systems  to  deal  with: 
the  plotting  or  display  system,  D;  the  target-fixed  or 
RATSCAT  system,  R ;  and  the  primed  system  of  a 
component,  P. 

D  is  connected  with  R  by  a  linear  transformation, 

C: 


xD  =  c „xR  -I-  ct2yR  +  c„zR 


yp  =  c2|X*  +  c22yR  +  c2}zR 


xk  +  ci2yR  +  c»zR 


where  the  caret  denotes  unit  vector. 

The  user  specifies  a  region  given  by  angular  inter¬ 
vals  for  the  display  polar  and  azimuth  angles  0D  and 
4>D,  respectively.  For  them,  he  must  calculate  the 
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corresponding  angles  in  the  R  system  (where 
shadowing  is  defined)  as  follows: 

cos  0R  -  c,3  sin  0D  cos  4>D  +  c2J  sin  0D  sin  4>d 

+  Cjj  cos  0D 

cost" 

_  cn  sin  0D  cos  4>o  +  c21  sin  0D  sin  4>D  +  c3l  cos  0D 
sin  0R 

(note  the  subscript  R  in  the  denominator  above) 
sin  <t>K 

_  C|2  sin  0D  cos  <j>D  +  c22  sin  0D  sin  <j>D  +  cK  cos  0D 
sin  0R 

(and  here  too). 

If 

cn  sin  dD  cos  <t>D  +  c2l  sin  0n  sin  4>D 
+  cJt  cos  0D  —  0, 

then  =  ±  x/2  when 

c,2  sin  6d  cos0o  +  c22  sin  0D  sin  4>D  +  c32  cos  0D 

is  qreater  than  or  less  than  zero,  respectively. 

If 

c,2  sin  0D  cos^0  +  c22  sin  0P  sin  <t>D 
+  cJ2  cos  6n  =  0, 


and  U  systems  (in  the  same  order  as  R  and  D  are 
connected  above).  It  is  now  identical  with  the  process 
that,  henceforth,  connects  the  P  and  R  systems. 
Furthermore,  since  the  U  system  is  now  eliminated, 
the  transformation. 


-1  0  0 

L  0  0  ij 


must  also  be  eliminated. 

To  recapitulate,  in  the  revised  arrangement  there 
are  only  the  three  coordinate  systems:  the  display 
system,  D ;  the  target-fixed,  or  RATSCAT,  system, 
R;  and  the  primed,  or  component-fixed,  system,  P. 
Henceforth,  the  three  instructions  TRANS,  VALID, 
and  ORGN  will  be  referred  to  the  R  system. 

Pitch,  Roll,  Yaw 

The  linear  transformation  C  can  be  expressed  as 
any  one  of,  or  a  combination  of,  pitch,  roll,  or  yaw, 
denoted  by  matrixes  P,  R,  and  Y,  respectively. 
Taking  the  original  position  of  D  to  be  coincident 
with  R  and  moving  D  while  keeping  R  fixed,  the 
angles  are  defined  as  follows: 

Pitch  through  an  angle  p  is  a  rotation  of  D  about 
y„,  positive  from  zH  toward  xK.  C  takes  the  special 
form  P,  where 


cos  p  0  -sinpl 
P=  0  1  0 

sinp  0  cos p  J 


Roll  through  an  angle  r  is  a  rotation  of  D  about  zR, 
positive  from  xR  toward  yR.  C  takes  the  special  form 
R  where 


then  <(>R  sb  0  or  t  when 


c„  sin  0D  cos  <>0  +  c2,  sin  0D  sin  <t>D 
+  c31  cos  0D 


is  greater  than  or  less  than  zero,  respectively. 

If  0R  «  0,  then  <t>R  is  indeterminate  because  the 
direction  of  interest  is  identical  with  z„ . 

As  was  indicated  above,  these  formulas  represent 
the  same  process  used  previously  to  connect  the  P 


cos  r  sin  r  0 

R  =  -  sin  r  cos  r  0 

0  0  1 


Yaw  through  an  angle  y  is  a  rotation  of  D  about 
xR,  positive  from  yR  toward  zR.  C  takes  the  special 
form  Y,  where 


1  0  0 
Y  -  0  cos.y  sin^  . 

.0  -  sin  y  cos  y 
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When  applied  successively,  P,  R,  and  Y  each 
represents,  in  some  order  (these  transformations  are 
not  commutative),  a  rotation  of  the  coordinate 
system  from  its  last  position  to  a  new  position.  The 
complete  transformation  C,  equal  to  a  product  of  P, 
R,  and  Y  in  some  order,  represents  a  transformation 
from  the  R  system  to  the  D  system. 


Note  that 


PR  = 


0  0-1 
0-1  0 
-10  0 


*  RP. 


Examples 

Consider  a  90*  pitch  of  axes,  followed  by  a  180* 
roll,  with  zero  yaw.  Then 


-1 

0 

0 

0 

0 

- 1 

0 

0 

1 

C  =  RP  = 

0 

-1 

0 

0 

1 

0 

= 

0 

-1 

Q 

0 

0 

1 

1 

0 

0 

1 

0 

0 

For  example,  with  the  conditions  0  =  60*,  0*  s  <f> 
<  180*,  the  transformation  RP  produces  a  portside 
viewing,  from  front  to  rear,  over  a  conical  cut  at  a 
30*  elevation  above  the  horizontal  plane.  In  contrast, 
the  transformation  PR  produces  a  portside  viewing, 
from  rear  to  front,  over  a  conical  cut  at  a  30* 
declination  below  the  horizontal  plane. 
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